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FOREWORD 


This  study  involves  the  palladium-boron  alloy  materials  prepared  at  the  Naval  Research 
Laboratory  (NRL)  by  Dr.  M.  Ashraf  Imam  (see  NRL/MR/6 1 70-96-7803 ,  January  9, 
1996).  This  new  material  was  developed  as  part  of  a  collaborative  program  with  NRL 
and  the  Naval  Air  Warfare  Center  Weapons  Division  (NAWCWD),  China  Lake,  that  was 
funded  by  the  Office  of  Naval  Research  (ONR).  Studies  at  NAWCWD  showed  that  the 
best  reproducibility  for  excess  power  was  obtained  using  the  palladium-boron  materials 
supplied  by  NRL  (see  NAWCWPNS  TP  8302,  September  1996).  The  new  experimental 
studies  described  in  this  report  were  conducted  by  Dr.  Melvin  H.  Miles  at  the  New 
Hydrogen  Energy  (NHE)  laboratory  in  Sapporo,  Japan.  Dr.  Melvin  H.  Miles  received  a 
six  month  appointment  as  a  Guest  Researcher  sponsored  by  the  New  Energy 
Development  Organization  (NEDO)  of  Japan.  Dr.  Melvin  H.  Miles  expresses  his 
appreciation  to  Dr.  N.  Asami  and  Mr.  K  Matsui  for  providing  him  with  this  research 
opportunity.  This  experiment  was  conducted  in  a  special  Dewar-type  colorimetric  cell 
silvered  at  the  top  that  was  developed  by  Drs.  Martin  Fleischmann  and  Stanley  Pons. 
The  detailed  analysis  of  the  experimental  data  presented  in  this  report  was  conducted  by 
Dr.  Martin  Fleischmann.  An  independent  method  of  data  analysis  developed  by  Dr. 
Melvin  H.  Miles  while  he  was  in  Japan  was  presented  in  his  NEDO  Final  Report  and 
shows  similar  trends  for  the  excess  heat  effect. 


LaAv  YPficsWh  iUim  '  Dr  Martin  Fleischmann 


DRS.  FLEISCHMANN,  MILES  AND  IMAM 


Executive  Summary 


The  cathode  for  the  heavy  water  electrolysis  experiment  was  synthesized  and  processed 
using  elemental  palladium  and  boron  (see  report  #  NRL/MR/  6170-7803)  at  the  Naval 
Research  Laboratory.  Washington  D.C.,  USA.  The  electrolysis  experiment,  as  a  part  of 
continuation  of  the  work  at  Naval  Air  Warfare  Center  Weapons  Division  (see  report  # 
NAWCWPNS  TP  8302).  was  carried  out  by  Dr.  Melvin  H.  Miles  during  his  stay  at  the 
New  Hydrogen  Energy  (NHE)  Laboratory.  Sapporo,  Japan.  In  this  experiment  an 
ICARUS-2  electrochemical  unit  and  data  acquisition  system  were  used  similar  to  those 
utilized  in  the  earlier  work  in  the  IMRA-Europe  Laboratories. 

The  experimental  data  and  the  preliminary  evaluation  [carried  out  by  members  of  the 
NHE  group]  is  contained  in  a  series  of  (kR')n  spreadsheets  where  (k|<')n  is  the  pseudo 
radiative  lower  bound  heat  transfer  coefficient  based  on  the  assumption  of  a  zero  rate  of 
excess  enthalpy  generation.  The  NHE  team  also  quoted  a  single  value  of  the  true  heat 
transfer  coefficient  determined  during  the  third  measurement  cycle  [Day  3  of  the 
experiment].  The  methodology  used  by  the  NHE  group  was  non-standard  [in  the  sense 
that  the  method  used  was  not  that  specified  for  the  ICARUS-system]:  it  appears  that  this 
coefficient  is  related  to  the  “true  integral  heat  transfer  coefficient,  (k|<'j362  based  on  an 
extrapolation  procedure  and  using  the  “forward  integration”  of  the  experimental  data. 

It  is  shown  that  the  value  of  (k|<')362  determined  is  incorrect:  it  is  evidently  too  small 
because  it  predicts  negative  rates  of  excess  enthalpy  generation  even  for  the  third 
measurement  cycle  [where  (kiff^  was  determined]  as  well  as  for  other  measurement 
cycles  in  the  experimental  sequence.  A  major  point  of  difficulty  was  evidently  the  very 
early  development  of  "positive  feedback”  [an  increase  of  the  rate  of  excess  enthalpy 
generation  with  temperature]  for  this  Pd-B  cathode.  In  the  general  case,  the  development 
of  such  “positive  feedback”  makes  it  impossible  to  evaluate  the  "true  heat  transfer 
coefficients”.  However,  in  this  experiment  the  effects  were  relatively  small,  and  it  is 
shown  that  such  effects  can  be  taken  into  account  in  the  determination  of  the  “true  heat 
transfer  coefficient,  (k|<')262,”  based  on  the  backw  ard  integration  of  the  experimental  data. 
The  evaluation  of  (k«')262  was  the  procedure  recommended  for  use  with  the  ICARUS- 
systems.  A  number  of  further  errors  in  the  evaluation  given  by  the  NHE  group  have  also 
been  corrected. 

It  is  also  shown  that  the  value  of  (kiffi^  agrees  with  the  maximum  value  of  (kR')n 
determined  in  the  early  parts  of  the  experiment  (here  on  Day  2]  before  the  onset  of  excess 
enthalpy  generation  but  after  the  completion  of  the  exothermic  charging  of  the  electrode. 
It  also  agrees  with  the  values  of  (kR')n  determined  on  Day  61  of  the  measurement 
sequence  during  which  there  was  a  near  zero  rate  of  excess  enthalpy  generation. 

The  use  of  the  correct  value  of  the  "true  heat  transfer  coefficient,  (k|<')2<,2  ’  show's  that 
there  were  low  levels  of  excess  enthalpy  generation  throughout  the  experimental 
sequence  [except  Day  61  ].  The  cell  was  evaporated  to  dryness  on  Day  68  of  the  sequence 
and,  during  this  day,  the  specific  rate  of  excess  enthalpy  generation  rose  to  ~27  W  cm'3. 
Measurements  on  Day  68  are  interpreted  in  detail,  and  it  is  shown  that  the  specific  rate  of 
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excess  enthalpy  generation  before  evaporation  to  dryness  agrees  with  that  immediately 
following  the  completion  of  this  step  which  is  an  example  of  “  Heal-after-Dealh”.  A 
further  example  of  such  “  Heat-afler-Dealh"  was  observed  on  the  (Inal  day  of  the 
experiment  [  Day  69]  as  well  as  on  the  lowering  of  the  current  dcnsitv  on  Day  26  of  the 
sequence. 

Compatisons  are  made  with  other  experiments  principally  in  a  set  of  footnotes. 
Suggestions  made  for  further  evaluations  could  form  addenda  to  this  report. 


Section  A 


A.l  Introduction 

This  experiment,  designated  as  FP2-97120402-M7c2  by  the  New  Hydrogen  Energy 
(NHE)  group  was  carried  out  by  M.H.  Miles  during  his  stay  in  1997-98  in  the  Sapporo 
Japan  Laboratory  of  the  NHE  group.  The  experiment  was  started  at  10:00  a.m.  on 
December  5,  1997  and  terminated  at  10:00  a.m.  on  February  12,  1998.  This  implies  that 
the  experiment  was  terminated  on  the  70th  day  following  the  start  of  the  experiment. 
1  he  electrode  used  had  the  form  of  a  cylindrical  rod  of  dimensions  4.71  mm  in  diameter 
and  20.1  mm  in  length.  This  yields: 

Volume  of  electrode  =  0.350  cm3  ^  1) 

Area  of  electrode  =  3.15  cm2  (A  2) 

This  electi ode  was  prepared  at  the  Naval  Research  Laboratory.  Washington,  D.C.  with  a 
composition  of  99.5%  Pd+0.5%  B  (wt%). 

Palladium-  Boron  Cathode 

Preparations  of  the  palladium  alloy  with  low  boron  contents  were  attempted  to  keep 
within  the  miscibility  gap.  Three  compositions  of  the  palladium-boron  alloy  were 
prepared  at  the  Naval  Research  Laboratory  using  an  arc  melter  with  w  ater-cooled  copper 
hearth.  Palladium  sponge  with  five-nine  purity  and  high  purity  boron  were  used  to 
produce  these  alloys.  The  three  compositions  had  nominal  boron  concentrations  of  0.75, 
0.50,  and  0.25  weight  percent  boron.  These  alloys  were  identified  by  their  nominal 
compositions.  The  glow-discharge  mass  spectroscopy  (GDMS)  analvses  of  the  three  as- 
prepared  alloys  showed  boron  contents  of  0.62.  0.38  and  0.18  weight  percent.  Processing 
ol  the  electrodes  introduced  copper  [<26  parts  per  million  by  weight],  tungsten  [<2.2 
parts  per  million  by  weight],  and  platinum  [<47  parts  per  million  by  weight]  The  lower 
boron  content  alloy  had  the  lowest  pick  up  of  these  elements.  The  cast  alloys  consisting 
of  in egulat  lod  shapes  were  swaged  and  rolled  to  sizes  depending  on  the  requirements  of 
the  experiment.  The  samples  were  annealed  for  two  hours  at" 650°C  resulting  in  an 
average  grain  size  of  90  pm. 

X-ray  diffraction  studies  were  carried  out  to  characterize  the  three  alloys  using  a  Phillips 
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diffractometer  with  a  generator  setting  of  50  kV,  30  mA  and  a  copper  target.  Two  distinct 
phases  of  the  same  cubic  structure  with  different  lattice  parameters  were  found  in  all 
three  alloys.  The  lattice  parameter  in  one  phase  remained  constant  with  changes  in  boron 
content  whereas  the  lattice  parameter  of  the  other  phase  increased  with  the  increase  in 
boron  content.  As  the  boron  content  increased,  the  fraction  of  one  phase  decreased  at  the 
expense  of  the  other  phase,  as  expected.  The  change  in  the  lattice  parameters  with  boron 
content  occurred  in  the  phase  where  the  starting  lattice  parameters  are  the  same  as  pure 
palladium. 

Experimental  Details 

The  Pd-0.5  wt%  B  electrode  was  cleaned  and  polished  using  the  standard  diamond  paste 
method  at  NHE.  The  specimen  was  observed  using  an  optical  microscope,  and  the 
surface  looked  clean  and  highly  polished  with  some  fine  circular  polishing  lines.  The 
electrode  was  spot  welded  on  the  side  to  a  platinum  lead  wire.  Quick  setting  epoxy  was 
used  to  cover  the  spot  welded  area,  the  top  of  the  electrode  as  well  as  the  end  of  the  glass 
tubing  containing  the  platinum  lead  wire.  The  cell  was  assembled,  placed  in  the  water 
bath,  and  connected  to  the  ICARUS  2  system.  The  cell  was  then  filled  with  90  cm  of  0.1 
M  LiOD  by  using  a  large  syringe  (50  cm3).  It  was  determined  that  82  cm  of  the  solution 
filled  the  cell  to  the  bottom  edge  of  the  silvered  portion.  This  mark  was  frequently  used  to 
determine  the  amount  of  D20  required  in  re-filling  the  cell.  The  D20  used  throughout  this 
experiment  was  from  Isotec,  Inc.  [99.9  atom%  D],  The  cell  dimensions  were  25.0  cm  in 
height  with  the  top  8.0  cm  silvered.  The  outside  diameter  of  the  cell  was  4.2  cm  with  an 
inner  diameter  of  about  2.5  cm. 

The  cell  design  [the  ICARUS- 1  Type]  is  illustrated  in  Fig.  A.l.  The  cell  number  was 
given  as  38  and  the  experiment  was  carried  out  in  position  A.2  [i.e.  in  position  2  of 
thermostat  tank  A]  using  an  ICARUS-2  Type  electrochemical  polarization,  control  and 
data  acquisition  system.  During  the  same  time  period,  the  Pd-B-Ce  alloy  prepared  by 
Naval  Research  Laboratory  was  run  in  position  A.  1  of  the  thermostat  tank  and  Pd-Ce  was 
run  in  position  A.3  in  two  other  Dewar  cells.  These  experiments  are  discussed  in  another 
report  [see  ICCF-8  Proceedings].  The  electrochemical  system  consisted  of  a  separate  Hi- 
Tek  DT2101  potentiostat  wired  up  as  a  galvanostat  for  each  cell.  These  potentiostats/ 
galvanostats  are  capable  of  delivering  currents  of  ±1A  at  output  voltages  up  to  -  +100V. 
A  separate  potentiostat/  galvanostat  was  also  used  to  deliver  constant  currents  to  the 
resistive  heater  used  to  calibrate  the  cell.  The  system  was  controlled  by  a  486  data 
acquisition  computer  which  also  controlled  an  Hewlett  Packard  44705A  multiplexer  and 
data  acquisition  system.  This  data  acquisition  system  was  on  an  IEEE-GPIB  bus  so  that  it 
would  be  anticipated  that  there  would  not  have  been  any  timing  errors  introduced  into  the 
measurements  [see  Footnotes  C.l  and  D.l],  The  system  is  fully  described  in  the 
associated  Handbook,  Document  Version  (2.0)  (February  1995)'  '[see  Footnotes  C.2  and 
D.2], 

The  three  ICARUS- 1  type  cells  were  filled  with  90.0  cm3  of  D2O+0.1M  LiOD  and 
inserted  into  a  large  water  thermostat  whose  temperature  was  independently  controlled  by 
a  Techne  TE-8A  stirrer/heater/regulator  unit  [see  Footnote  C.3  and  Figure  A. 5].  The 
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water  thermostat  was  in  turn  maintained  in  a  room  whose  temperature  was  generally  0  to 
2°C  higher  than  that  of  the  water  thermostat.  This  study  will  focus  only  on  the  Pd-B 
experiment  positioned  in  the  center  of  the  bath  [position  A. 2],  It  should  be  noted  that  the 
LiOD  concentration  for  this  cell  was  increased  to  0.2  M  on  Day  53  of  this  study. 

At  the  7th  International  Conference  on  Cold  Fusion  19th-24th  April  1998  we  [i.e.  M.H. 
Miles  and  M.  Fleischmann]  decided  that  a  report  on  M.H.  Miles'  visit  to  the  NHE 
Laboratories  should  include  an  assessment  of  the  Pd-B  experiment  using  as  far  as  is 
possible  the  ICARUS  Methodology.  We  decided  that  this  assessment  should  include: 

(i)  comments  on  differences  between  the  execution  of  the  experiments  carried  out 
in  the  NHE  Laboratory  and  the  procedures  recommended  for  use  with  the 
ICARUS-1  and  ICARUS-2  Systems  (l)(2); 

(ii)  a  commentary  on  the  calibration  of  the  ICARUS-1  Type  cell  used  in  the 
Pd-B  experiment; 

(iii)  an  assessment  of  the  presence  [or  absence]  of  excess  enthalpy  generation  in 
the  Pd-B  experiment; 

(iv)  an  investigation  of  the  presence  [or  absence]  of  the  effects  of  "positive 
feedback";  [e.g.  see  (6)>(7)(8)]; 

(v)  An  investigation  of  the  presence  [or  absence]  of  the  effects  of  "Heat  after 
Death";  [e.g.  see  (9)’(10)]; 

(vi)  an  assessment  of  the  implications  of  the  evaluations  of  the  Pd-B  experiment 
for  further  research  in  this  field. 


There  are  also  some  secondary  reasons  for  choosing  the  Pd-B  experiment  for  further 
detailed  evaluations  [see  Footnote  C.4], 


This  Report  is  divided  into  the  following  Sections: 


Section  A: 
Section  B: 

Section  C: 
Section  D: 

Section  E: 


the  main  text. 

a  description  of  the  methods  of  data  evaluation  as  characterized  by  the 
relevant  heat  transfer  coefficients. 

a  set  of  Footnotes  and  Comments  relating  to  Sections  A  and  B. 
a  further  set  of  Footnotes  and  Comments,  which  perhaps  should  have  a 
restricted  circulation. 

A  set  of  footnotes  and  comments  whose  circulation  should  be  restricted. 


A.2  Some  Important  Preliminary  Definitions,  Descriptions  and  Evaluations 

The  various  radiative  heat  transfer  coefficients  used  in  this  report  are  defined  as 


(kR')ijj 

where  R  is  radiative  heat  transfer,  i  =  1,2,3  denotes  "differential",  "backward  integration" 
and  "forward  integration"  respectively,  1=1,2  denotes  “lower  bound”  and  “true” 
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respectively  [see  Footnote  C.12]  and  j,  if  used,  denotes  the  time  period  of  the 
measurement  cycle  as  given  below 

j=5,  times  somewhat  above  the  origin 

j=6,  times  somewhat  above  ti  [application  of  calibration  pulse] 

j=7,  times  somewhat  above  t2  [cessation  of  calibration  pulse] 

j=8,  combination  of  times  for  j=6  and  j=7 

These  values  of  “j”  cover  time  periods  where  the  cell  temperature  shows  large  changes. 
In  addition,  the  radiative  heat  transfer  coefficient  (tcRo),^  is  also  frequently  used  and  is 
obtained  from  extrapolations  of  the  straight  line  forms  of  the  calorimetric  equations  [see 
Equation  B.  1 1  and  B.  1 2], 

The  protocol  used  for  the  experiment  was  as  follows: 

(i)  the  electrode  was  first  of  all  polarized  for  two  days  without  any 
application  of  calibration  pulses; 

(ii)  on  the  third  day  [and  on  all  subsequent  days  including  days  68  and  69 
when  the  cell  had  reached  dryness]  calibration  pulses  were  applied; 

(iii)  changes  of  current  density  were  made  frequently  at  times  close  to  10:00 
am  of  the  measurement  cycles;  these  changes  of  current  density  are 
shown  in  Table  A.l  and  in  Fig.  A.2B  found  at  the  back  of  the  text; 

(iv)  the  cell  was  "topped  up"  with  D20  whenever  this  was  judged  to  be 
necessary  at  10:00  am,  the  start  time  of  all  the  measurement  cycles;  the 
cell  was  then  left  to  equilibrate  for  9  hours  followed  by  the  application  of 
calibration  pulses  of  6  hour  duration  [the  start  time  of  all  the  calibration 
pulses  was  7:00  pm];  the  cell  was  then  again  left  to  equilibrate  for  a 
further  9  hours  before  reaching  10:00  am  of  the  next  day  of  the 
experimental  sequence; 

(v)  as  is  evident  from  (iv),  the  duration  of  the  measurement  cycles  was  24 
hours;  the  addition  of  D20  and  the  changes  in  the  current  density  were 
always  made  by  M.  H.  Miles  at  his  discretion.  Everything  else  was 
controlled  by  the  NHE  computer  program. 

This  protocol  differs  substantially  from  that  specified  for  the  operation  of  the  ICARUS-1 
and  ICARUS  -2  Systems,  which  was  as  follows (1)(2>; 

(ia)  the  electrodes  were  to  be  polarized  for  4  days  [i.e.  two  measurement 
cycles,  see  (va)  below]  without  any  application  of  calibration  pulses; 

(iia)  on  the  5th  day  [i.e.  for  the  third  measurement  cycle]  and  for  nine  further 
measurement  cycles  calibration  pulses  were  to  be  applied  as  specified  in 
(iva)  below;  this  was  to  be  followed  by  two  further  measurement  cycles 
without  the  application  of  calibration  pulses  and,  in  turn,  by  ten  further 
cycles  with  calibration  pulses.  A  total  experiment  duration  of  48  days  was 
therefore  specified  for  the  initial  phase  of  the  work. 
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(iiia)  the  initial  experiments  were  to  be  carried  out  at  a  single,  low  current 
density,  typically  <  250  mA  cm'2  or  788  mA  for  this  Pd-B  electrode;  in 
later  experiments  a  single,  low  current  density  was  to  be  applied  for 
various  initial  durations  followed  by  a  raising  of  the  current  density  to 
values  typically  >  1  A  cm'2;  this  protocol  was  in  broad  accord  with  that 
used  in  previous  investigations  (b)'(l6);  changes  of  current  density  were  to 
be  made  at  the  beginning  of  the  measurement  cycles.  The  large  size  of  the 
Pd-B  electrode  [A=3.5  cm2]  limited  the  current  density  to  about  320 
mA/cm2. 

(iva)  cells  were  to  be  "topped  up"  at  the  start  of  each  measurement  cycle;  the 
cells  were  then  to  be  left  to  equilibrate  for  12  hours  and  calibration  pulses 
of  12  hour  duration  were  then  to  be  applied;  the  cells  were  then  again  to 
be  left  to  equilibrate  for  a  further  24  hours  so  as  to  reach  the  start  of  the 
next  measurement  cycle. 

(va)  as  is  evident  from  (iva),  the  duration  of  the  measurement  cycles  was  to  be 
48  hours. 

Further  important  aspects  of  the  protocols  specified  for  the  ICARUS-1  System  are  given 
in  Footnote  C.5.  These  protocols  were  designed  with  two  major  aims:  firstly,  to  verify 
the  modeling  of  the  cells  which  is  the  basis  of  the  methods  of  data  evaluation;  secondly, 
to  facilitate  these  methods  of  data  evaluation.  Both  of  these  aspects  are  outlined  in 
Sections  B.2  and  B.3  together  with  the  associated  Comments  and  Footnotes. 

Needless  to  say,  with  the  progression  of  the  work,  parts  of  these  protocols  were  expected 
to  become  redundant.  Furthermore,  the  protocol  (ia)  -  (va)  had  been  designed  for  the 
execution  and  evaluation  of  appropriate  "blank  experiments,"  and  it  was  unlikely  that  this 
protocol  would  be  entirely  suitable  for  experiments  on  the  Pd-D20  type  systems 
[especially  experiments  in  which  there  were  changes  in  the  current  density].  We  also 
have  to  take  special  note  of  an  unsatisfactory  situation.  With  the  exception  of  some 
fragmentary  analyses  of  data  sets  collected  during  the  start-up  of  the  NHE  project (3)  we 
have,  as  of  now,  no  comprehensive  data  sets  and  evaluations  of  "blank  experiments"  [see 
Footnote  C.27],  We  therefore  have  to  rely  on  experiments  carried  out  at  IMRA  Europe 
and  on  the  analyses  of  simulations  e.g.  see  Reference  13  for  this  aspect  of  the  work;  i.e. 
we  are  unable  to  investigate/validate  the  operation  of  the  ICARUS-1  and  -2  Systems 
installed  in  Sapporo. 

As  far  as  the  second  aspect  is  concerned,  we  have  to  take  note  of  the  fact  that  in  view  of 
the  changes  in  the  protocols  made  by  both  M.  H.  Miles  and  the  group  at  NHE  [(iv)  and 
(v)  as  against  the  specified  (iva)  and  (va)],  we  are  unable  to  apply  the  simplest  method  of 
calibration  leading  to  (kR’)2  [see  Section  B.2]  in  view  of  the  inadequate  relaxation  of  the 
temperature  following  the  application  of  the  heater  pulses  [see  Footnotes  C.28  and  C.29]. 

We  now  consider  further  the  major  differences  between  the  operation  of  the  Pd-B 
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experiment  and  the  conditions  used  in  previously  reported  investigations,  e.g.  see 
References  12,  15  and  16.  Apart  from  the  frequent  changes  of  current  density  [see 
Footnote  C.30,’  Fig.  A.2B],  we  can  see  that  these  current  densities  were  mostly  in  the 
vicinity  of  the  threshold  value  required  for  the  onset  of  the  phenomenon  of  excess 
enthalpy  generation  (12)  [see  Footnote  C.31].  Furthermore,  the  cell  temperatures  were 
mostly  below  the  level  required  for  the  onset  of  "positive  feedback",  Fig.  2.A  (  )(  ’  and 
which  leads  to  a  marked  increase  in  the  rates  of  excess  enthalpy  generation  [see 
Footnotes  C.32  and  C.33].  The  conditions  in  the  cell  therefore  remained  in  the  vicinity  of 
the  region  of  onset  of  "positive  feedback"  and,  under  these  conditions  we  would  not 
expect  to  see  a  marked  build  up  in  the  rate  of  excess  enthalpy  generation  [see  Footnote 
C.34], 

Consideration  of  Fig.  A.2B  also  allows  us  to  decide  on  the  measurement  cycles  likely  to 
provide  examples  of  "Heat-after-Death"  [objective  (v)  of  this  investigation]  as  was 
pointed  out  in  the  original  investigation  <9)'(lll).  It  would  be  expected  that  this  phenomenon 
would  be  observable  under  several  distinct  conditions,  which  include  [using  the  original 
classification  numbers]: 

Cell  full:  cell  operated  at  intermediate  temperatures;  cell  current  then 

reduced  in  stages;  [original  number  1] 

Cell  empty:  cell  allowed  to  boil  dry;  cell  then  maintained  at  the  rail  voltage  of 
the  galvanostat;  [original  number  5] 

Cell  empty:  cell  allowed  to  boil  dry;  cell  disconnected  from  the  galvanostat. 
[original  number  6] 

Consideration  of  the  "hard  copy"  of  the  data  sets  shows  that  condition  #5  applies  to  part 
of  Day  68  of  the  sequence  measurement  cycles  [see  Section  A. 8]  while  condition  #6 
applies  to  part  of  Day  69  of  this  sequence  [see  Section  A.9].  Consideration  of  Fig.  A.2B 
shows  that  condition  #1  is  likely  to  apply  to  several  of  the  measurement  cycles.  The 
effects  would  be  expected  to  be  most  marked  for  Day  26  [reduction  of  the  cell  current 
from  above  to  below  the  threshold  for  excess  enthalpy  generation,  reduction  in  cell 
temperature  from  above  the  level  for  the  onset  of  "positive  feedback  to  below  this  level]. 
Attention  is  confined  in  this  Report  to  this  particular  day  [see  Section  A.  10]  although  it  is 
evident  that  there  are  several  further  regions  of  time  which  might  well  give  examples  of 
"Heat-after-Death"  following  Scenario  1. 

In  this  Section  we  should  also  consider  a  further  difference  between  the  protocols  for  Pd- 
B  experiment  and  those  used  in  earlier  studies;  namely,  the  schedules  of  addition  of  D2O 
to  make  up  the  losses  due  to  electrolysis.  Table  A.l  and  Fig.  A. 3  illustrate  the  effect  on 
the  change  in  volume  of  the  electrolyte  in  a  cell  in  an  hypothetical  experiment  carried  out 
first  at  a  cell  current  of  200  mA  for  29  days  followed  by  a  cell  current  of  500  mA  and 
with  a  daily  schedule  of  D20  additions.  It  can  be  seen  that  the  mean  volume  of  the 
electrolyte  falls  by  some  1 .21  cm3  between  the  two  time  regions.  We  can  estimate  that 
this  would  cause  a  decrease  of  the  mean  value  of  (kR’)i  2  by  0.15%  or  of  (kR’)22  by  - 
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0.075%.  Such  small  changes  are  close  to  the  error  limits  quoted  for  the  instrumentation 
and  can  normally  be  neglected.  However,  the  magnitude  of  the  changes  are  above  the 
error  limits  which  can  actually  be  achieved  [e.g.  see  (13)  and  Section  A. 4]  and  should  be 
taken  into  account  in  evaluations  carried  out  at  the  maximum  achievable  precision  and 
accuracy. 

Fig.  A.4  shows  the  effects  of  D2O  additions  as  actually  used  in  this  Pd-B  experiment. 
This  figure  assumes  that  the  D20  is  lost  only  by  electrolysis.  Experimental 
measurements,  however,  show  an  additional  4%  of  the  D20  was  carried  out  of  the  cell  by 
the  gas  stream.  The  exact  volumes  of  D20  were  recorded  throughout  this  experiment  [see 
Footnote  C.35],  It  can  be  seen  that  the  expected  changes  in  (kR')i2  due  to  the  electrolyte 
volume  now  lie  between  -0.3  and  +0.6%  and  of  (kR')22  between  -0.15  and  +0.3%,  changes 
which  should  certainly  again  be  taken  into  account.  These  numbers  will  change 
somewhat  due  to  the  additional  loss  of  D20  by  the  gas  stream  [see  Footnote  C.36]. 

The  schedule  of  additions  leads  to  an  important  conclusion.  We  find  that  by  Day  67  the 
total  volume  of  D20  added  was  262.5  cm3  whereas  the  total  volume  electrolyzed  was 
253.3  cm3.  Earlier  in  this  experiment  on  Day  39,  the  total  volume  of  D20  added  was 
1 12.5  cm3  whereas  the  total  volume  electrolyzed  was  108.2  cm3.  It  is  evident  that  the 
volume  of  D20  added  is  consistently  about  4%  larger  than  the  volume  electrolyzed  and, 
therefore,  there  could  not  have  been  any  recombination  of  the  deuterium  and  oxygen 
produced  by  electrolysis.  This  is  in  agreement  with  earlier  measurements  (12)  [see 
Footnote  C.37]  and  numerous  measurements  by  other  authors  [see  Footnote  C.38], 

The  horizontal  lines  in  Figs.  A.3  and  A.4,  delineate  the  volumes  of  D20  below  and  above 
which  we  would  expect  the  electrolyte  level  to  fall  below  the  base  of  the  silvering  in  the 
upper  part  of  the  cell,  Fig.  A.l,  or  to  approach  the  base  of  the  Kel  F  plug  at  the  top  of  the 
cell.  The  results  of  Figure  A.4  are  calculated  theoretically  by  assuming  that  D20  is 
removed  only  by  electrolysis.  Experimentally,  the  electrolyte  was  observed  to  fall  below 
the  base  of  the  silver  line  prior  to  the  D20  addition  only  on  Day  25  (-1  cm3).  Day  41  (-1 
cm3).  Day  48  (-2  cm3),  Day  50  (-1.5  cm').  Day  53  (-0.5  cm3).  Day  55(-2  cm3),  and  Day 
67  (-6  cm3).  It  can  be  seen  that  the  electrolyte  level  remained  within  the  space  defined  by 
this  silvered  portion  throughout  most  of  the  measurement  cycles  [see  Footnote  C.39]. 
However,  we  can  see  that  the  electrolyte  level  may  have  approached  the  base  of  the  Kel  F 
plug  at  the  start  of  several  of  the  measurement  cycles  following  the  "topping  up"  of  the 
cells.  In  work  carried  out  at  IMRA  Europe,  it  has  been  established  that  such  "overfilling" 
of  the  cells  leads  to  an  anomalous  increase  of  the  pseudo  radiative  heat  transfer 
coefficients  by  4  -  5%  of  the  values  which  apply  at  the  mean.  This  increase  in  (kR‘)  is 
almost  certainly  due  to  an  increase  in  the  conductive  contribution  through  the  Kel  F  plug 
to  the  overall  heat  transfer  from  the  cell  [c.f.  equation  (B.2)].  Experimentally  the 
"overfilling"  of  the  cell  in  this  Pd-B  experiment  occurred  on  Day  26  (92.0  cm3).  Day  43 
(94.0  cm3).  Day  50  (93.5  cm3),  and  on  Day  57  (94.0  cm3).  This  was  done  intentionally 
because  of  long  weekends  when  it  was  not  possible  to  enter  the  laboratory. 

This  type  of  "overfilling"  behavior  may  apply  also  to  Day  61  which  is  a  measurement 
cycle  for  which  we  can  get  important  confirmatory  evidence  of  the  "true  heat  transfer 


8 


coefficient"  which  applies  to  the  operation  of  the  cell  [see  Fig.  A.20,  Section  A.5]. 
Fig.A.20  shows  the  expected  increase  in  (kR')ii  at  times  close  to  the  "topping  up"  of  the 
cell.  Another  explanation  is  that  the  large  addition  of  D20  (9.0  cm3)  at  room  temperature 
[23°C]  on  Day  61  suddenly  cooled  the  cell  contents,  thus  (kR')n  is  too  large  during  the 
cell’s  recovery  from  this  cooling  shock. 


Finally,  we  also  make  a  number  of  preliminary  assessments  of  the  form  of  the 
temperature-time  and  cell  potential-time  series  for  Day  3,  i.e.  the  third  measurement 
cycle  of  this  experiment,  Figs.  A.5  -  A.7.  The  data  for  this  day  are  of  special  importance 
because  the  group  at  NHE  has  quoted  a  value  of  the  "true  heat  transfer  coefficient"  as 
given  by  their  method  of  evaluation  for  this  day.  This  value  of  the  "true  heat  transfer 
coefficient"  was  then  used  in  the  evaluation  of  all  the  measurement  cycles.  The 
evaluation  by  NHE  is  considered  further  in  Section  A. 3  while  Section  A. 4  gives  the 
application  of  the  ICARUS  Methodology  to  this  particular  data  set. 

Columns  1  -  6  of  Table  A.2  give  the  "raw  data"  for  the  third  measurement  cycle  [the 
remaining  columns  are  discussed  in  Section  A. 3].  We  can  immediately  draw  a  number  of 
important  conclusions.  Thus  Fig.  A.5  gives  a  plot  of  the  temperature  of  the  water  bath 
versus  time  for  the  first  32,400  s  of  the  measurement  cycle  [the  period  0  <  t  <  ti 
preceding  the  application  of  the  heater  calibration  pulse]  while  Fig.  A. 6  gives  plots  of  the 
cell  temperature  versus  time  for  the  same  period  and  for  both  positions  in  the  cell  where 
the  temperature  was  measured,  see  Fig.  A.  1.  It  is  evident  that  the  "noise'  level  of  the 
measurements  in  the  water  bath  [o  =  0.0088  K,  mean  =  295.198  K]  is  much  higher  than 
that  of  the  measurements  of  the  cell  temperature,  Fig.  A.6.  This  difference  [a  =  0.0027 
K]  is  to  be  expected  because  the  water  bath  is  controlled  by  a  single  thermal  impedance 
whereas  the  cell  is  controlled  by  two  impedances  in  series.  At  the  same  time,  the  "noise" 
in  the  temperature  of  the  water  bath  is  much  higher  than  that  in  the  original 
measurements  with  the  ICARUS-Sy stems  [a  =  0.003  K] (l2)  and,  in  our  experience,  such 
an  increase  is  due  to  inadequate  control  of  the  room  temperature.  The  room  temperature 
was  usually  0-2  K  above  the  bath  temperature,  i.e.  22  to  24  °C.  Typical  variations  in  the 
temperature  of  this  room  are  shown  elsewhere  [Journal  of  Electroanal.  Chem.,  Vol.  482, 
pp.  55-65,2000], 

It  will  be  evident  that  the  "noise"  in  the  measurements  of  the  temperature  of  the  water 
bath  is  one  factor  which  will  limit  the  precision  of  the  "lower  bound  heat  transfer 
coefficients,  (kR')i  i  via  its  effect  on  the  temperature  function,  fi(0),  see  Section  B.3.  The 
value  a  -  0.0088  K  is  outside  the  range  specified  for  the  ICARUS- 1  system  if 
measurements  are  made  at  low  cell  temperatures.  Table  A.2.  By  contrast,  the  "true  heat 
transfer  coefficients"  are  not  affected  by  such  fluctuations  because  the  temperature 
function  f2(0)  [see  Section  B.3]  is  determined  by  the  cell  temperature  alone. 

It  can  be  seen  that  the  variation  with  time  of  the  cell  temperature  measured  at  the  two 
positions  in  the  cell,  is  systematic,  Fig.  A. 6.  Moreover,  it  is  clear  that  there  is  a 
systematic  difference  in  temperature  between  the  two  positions,  which  must  be  due  to 
either  one  or  two  errors  in  the  calibration  [see  Footnote  C.40],  For  these  measurements 
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we  obtain  mean  [0si,Ort  thermistor  -  0|ong  thermistor]  =  0.0045  K  and  Cohort  thermistor  "  0long  thermistor] 
=  0  0027  K.  The  mean  gives  an  indication  of  the  accuracy  in  (kR')n  which  we  can  expect 
to  achieve  The  error  ~  0.05%  is  somewhat  above  the  target  for  the  precision  of  the 
measurements,  errors  <  0.01%,  which  is  hardly  surprising.  The  standard  deviation  gives 
double  the  value  of  the  expected  standard  deviation  for  the  measurements  with  one 
thermistor.  We  can  see  that  this  value,  ~  0.00135  K,  will  not  affect  the  accuracy  ol  the 
determination  of  any  version  of  the  "true  heat  transfer  coefficient  .  However,  we  shoul 
note  that  it  is  evidently  desirable  to  calibrate  the  thermistors  so  that  we  can  make  the 
temperature  measurements  to  within  +  0.001  K  [see  Footnote  C.41]. 

Differences  in  temperature  between  those  given  by  the  "short  thermistor"  and  "long 
thermistor"  will  be  considered  further  in  Section  A.7  dealing  with  Day  68  as  the  cell  is 
being  driven  to  dryness  and  in  Section  A.9  dealing  with  "Heat-after-Death"  on  Day  69. 

Finally,  we  consider  the  plots  of  the  "raw  data",  for  Day  3,  Fig.  A.  7.  We  can  see 
immediately  the  inadequacy  of  restricting  the  calibration  pulse  to  6  hours  because  the 
temperature  has  not  relaxed  to  equilibrium  in  this  time  period  [see  Footnote  C.42]. 
However,  in  this  particular  case  there  is  an  evident  complication  because  of  the  very  early 
establishment  of  "positive  feedback".  This  effect  can  be  seen  most  directly  from  the 
delayed  relaxation  of  the  temperature  to  the  base  line  following  the  cessation  of  the  heater 
calibration  pulse  [the  base  line  is  given  by  the  extrapolation  of  the  0  - 1  series  observed 
before  the  application  of  the  calibration  pulse].  Evidently,  the  raising  of  the  cell 
temperature  by  the  calibration  pulse  has  led  to  an  increase  in  the  thermal  output  from  the 
cell  which  persists  following  the  termination  of  the  calibration,  i.e.  a  form  of  positive 
feedback"  [see  Footnote  C.43].  The  calibration  of  such  a  system  can  obviously  only  be 
achieved  with  many  restrictions  and  with  great  difficulty  [see  Footnote  C.44]. 


The  interpretation  of  Fig.  A.7  will  be  considered  further  in  Sections  A. 3  and  A. 4. 

A.3  The  NHE  Interpretation  of  This  Pd-B  Experiment  (M7c2) 

As  has  already  been  pointed  out,  the  NHE  interpretation  of  the  Pd-B  experiment  rests  on 
the  determination  of  the  "true  heat  transfer  coefficient"  on  Day  3  of  the  measurement 
cycles.  Apart  from  the  citation  of  the  value  of  this  coefficient  (0.793504  x  1  O'  WK'  )  in 
the  header  for  the  spreadsheet  for  Day  1,  the  information  given  by  NHE  is  contained  in  a 
set  of  spreadsheets,  which  appear  to  be  related  to  the  (kR')n  spreadsheets  of  the  ICARUS 
Methodology  for  analyzing  the  data.  The  description  of  these  spreadsheets  is  contained 
in  Footnote  C.45.  In  addition  to  the  description  contained  in  this  Footnote,  we  have  to 
take  note  of  the  following  observations: 

a)  it  is  not  clear  how  the  value  of  the  "true  heat  transfer  coefficient"  was 
determined  nor  which  of  the  definitions  of  the  heat  transfer  coefficients 
may  have  been  used  [see  Section  B.3].  The  (kn')- 1  -  spreadsheets  [or 
(kR’)3i  spreadsheets]  (l>,(2)  which  would  have  allowed  an  investigation  of 
this  value  [and  which  were  part  of  the  ICARUS  data  processing  system] 
were  not  given  [see  also  Section  A.4].  However,  it  is  likely  that  this  was 
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the  coefficient  (kR')32  and  it  will  be  assumed  here  that  this  was  the  case, 
i.e.  we  will  assume  that  the  values  of  the  excess  enthalpies  given  in 
columns  (12)  and  (18)  of  the  spreadsheets  were  based  on  calculations 
using  the  single  value 

(kR')32  =  0.793504  x  10'9  WK4 

b)  it  is  also  not  clear  to  what  extent  the  values  of  the  "true  heat  transfer 
coefficient"  and  of  the  excess  enthalpies  may  have  been  affected  by  the 
value 

CpM  =  490  JK'1 

used  in  the  calculations.  Values  as  high  as  this  applied  to  cells  used  prior 
to  1992,  and  the  Handbooks  for  the  ICARUS  Systems  contained 
instructions  for  changing  this  [and  other]  parameter(s)  depending  on  the 
value  found  using  the  methods  of  evaluation  outlined  in  the  Handbooks 
(1),(2)  [see  also  Sections  A. 4  and  B.3]. 

It  should  be  noted  that  the  "guesstimate"  of  the  water  equivalent  of  the  cell 

is: 

CM-  contribution  of  D20  in  the  electrolyte  +  contribution  of  the 
glass  in  the  inner  cell  wall  =  (419  +  31)  JK'1  =  450  JK  1 

The  remaining  components  of  the  cell  [LiOD,  metals,  glass  framing, 
heater,  thermistor,  a  proportion  of  the  Kel-F  plug]  will  contribute  only  a 
small  additional  term  to  CpM.  It  follows,  therefore,  that  observations  of 
CpM  far  above  or  below  450  JK'1  indicate  malfunctions  of  the  methods  of 
data  evaluation. 

(c)  as  has  been  noted  elsewhere  [see  Section  B.2],  the  values  of  the  rates  of 
evaporative  cooling,  columns  (9)  and  (16),  given  in  the  spreadsheets 
cannot  be  calculated  using  the  instructions  given  in  the  Handbooks  for  the 
ICARUS- 1  and  -2  Systems (l)  (2).  The  differences  are  not  important  at  low 
temperatures  [such  as  those  which  apply  to  Day  3  of  the  measurement 
cycles.  Table  A.2]  but  become  significant  at  temperatures  close  to  the 
boiling  point.  However,  at  such  elevated  temperatures  other  factors 
neglected  in  the  calculations  carried  out  by  NHE  become  even  more 
important  [see  Sections  B.2  and  A. 7]. 

(d)  it  is  apparent  that  the  enthalpy  inputs  given  in  column  8  of  the  NHE 
spreadsheets  have  been  calculated  using  1.54  V  as  the  thermoneutral 
potential,  E°H,  whereas  most  other  authors  have  used  the  value  1.527  V. 
The  circumstances  leading  to  our  choice  of  the  value  1.54  V  have  been 
described  elsewhere  [see  Footnote  C.46],  The  value  1.527  V  has  been 
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used  in  deriving  Tables  A. 3  -  A. 9). 

While  considering  the  values  of  the  enthalpy  inputs,  column  (8)  of  the 
NHE  spreadsheets,  it  should  be  noted  that  these  have  been  given  to  4 
significant  figures  whereas  they  should  have  been  given  to  5  significant 
figures.  It  is  not  clear  at  the  present  time  whether  the  values  of  (k.R')n 
listed  in  the  spreadsheets  were  calculated  using  the  enthalpy  inputs 
rounded  to  4  figures  or  whether  a  higher  precision  was  maintained  in  the 
calculations. 

(e)  the  most  serious  shortcoming  of  the  NHE  calculations  is  that  the  input  due 
to  the  calibration  heater  has  been  entered  as  zero  rather  than  the  actual 
value  given  separately  as  0.25000  W  [see  Footnote  C.47],  In  the  procedure 
used  by  NHE  (ll)the  "lower  bound  heat  transfer  coefficient,  (kR')n"  is 
calculated  with  this  assumed  zero  enthalpy  input,  and  it  is  then  assumed 
that  the  magnitude  of  the  enthalpy  input,  can  be  recovered  together  with 
any  rate  of  excess  enthalpy  generation  by  using  this  derived  "lower  bound 
heat  transfer  coefficient"  together  with  the  "true  heat  transfer  coefficient, 
(kR’)32  and  fi(0).  Let  us  assume  first  of  all  that  such  a  procedure  is  correct. 
Then  we  can  see  an  immediate  disadvantage  as  compared  to  the  method 
outlined  for  the  ICARUS  Systems  in  that  we  are  unable  to  determine 
whether  (kR')n  during  the  period  of  the  application  of  the  calibration  pulse 
in  ti  <  t  <  t2  is  the  same  as  for  t  <  t|,  or  t  >  t2  [see  Footnote  C.48].  The 
data  derived,  e.g.  see  Fig.  A. 8  below,  are  certainly  further  degraded  by 
using  an  incorrect  value  of  CpM.  However,  in  actual  fact,  the  procedure 
used  by  NHE  is  invalid  as  has  been  pointed  out  in  a  Report  (4)  and  in 
subsequent  correspondence  [see  Footnote  C.49],  It  is  difficult  to  see  why 
the  straightforward  procedure  outlined  in  the  Handbook  for  the  ICARUS-1 
System (2)  was  not  followed. 

Notwithstanding  the  reservations  (a)-(e)  about  the  evaluations  carried  out 
by  NHE,  columns  7-10  of  Table  A. 2  summarize  these  evaluations 
[columns  12,  14,  18  and  19  of  the  (kR')n  spreadsheets  given  by  NHE] 
while  the  derived  data  in  columns  11-14  of  Table  A. 2  are  also  based  on 
the  NHE  evaluations.  Columns  15-17  of  Table  A.2  make  a  correction 
for  (e);  the  remaining  required  corrections  have  been  made  for  all  the 
derived  data  given  in  Tables  A.3-A.9,  Sections  A.4-A.7  and  A.  10. 

We  consider  next  the  values  of  the  "lower  bound  heat  transfer  coefficient, 
(kR')i  i,”  compared  to  the  "true  heat  transfer  coefficient"  as  given  by  the 
NHE  evaluation.  Columns  7  and  9  of  Table  A.2  and  Fig.  A. 8,  as  well  as 
the  values  of  the  rates  of  excess  enthalpy  generation,  columns  8  and  10 
and  Fig.  A. 9.  Table  A.2  also  includes  the  1 1 -point  averages  of  109  (kR')n 
and  Qcxccss  in  columns  13  and  1 1  and  the  further  6-point  averages  of  109 
(kp  ')n  and  Qcxccss  *n  columns  14  and  12  respectively  using  the  data 
evaluated  with  the  0  -  /  series  as  given  by  the  "short  thermistor".  We  can 
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only  conclude  from  these  data  that  the  evaluations  are  incorrect  based  on 
the  following  evidence: 


(f)  it  is  impossible  for  the  "true  heat  transfer  coefficient,  (kR')32"  to  be  smaller 
than  the  "lower  bound  heat  transfer  coefficient,  (kR')n’‘  because  the  lower 
bound  value  is  based  on  the  assumption  that  there  is  a  zero  rate  of  excess 
enthalpy  generation  in  the  cell.  The  type  of  difference  seen  in  Table  A.2 
and  Fig.  A. 8  could  only  arise  if  the  cell  was  endothermic,  and  the 
endothermicity  has  already  been  fully  taken  into  account  using  the 
thermoneutral  potential.  Any  additional  endothermicity.  therefore,  requires 
that  the  cell  operates  as  a  spontaneous  refrigerator,  and  this  violates  the 
Second  Law  of  Thermodynamics. 

(g)  the  pronounced  variation  of  the  "lower  bound  heat  transfer  coefficient, 
(kR’)n"  with  time  following  the  application  of  the  heater  calibration  pulse 
at  t  =  ti  and  its  cessation  at  t  =  t2  implies  at  the  very  least  that  the  "raw 
data"  have  been  evaluated  using  an  incorrect  value  of  the  water  equivalent, 
CpM,  of  the  cell. 

(h)  the  excess  enthalpy  given  by  the  NHE  evaluation  is  apparently  negative 
both  for  t  <  ti  and  t  >  t2  which  is  a  further  illustration  of  the  apparent 
violation  of  the  Second  Law  of  Thermodynamics,  cf.  (f)  above. 

(i)  it  has  been  maintained  (ll)  that  the  NHE  evaluation  recovers  the  magnitude 
of  the  heater  calibration  pulse,  AQ,  during  its  period  of  application,  ti  <  t  < 
t2,  together  with  any  rate  of  excess  enthalpy  generation.  Fig.  A.9  shows 
that  this  is  incorrect:  the  values  of  the  rates  of  excess  enthalpy  generation 
[which  here  include  the  enthalpy  input  to  the  calibration  heater]  are  less 
than  AQ  in  the  period  t|  <  t  <  t2  if  we  take  QCXcess  =  0  as  the  base  line.  If 
we  fix  the  base-line  at  the  level  of  the  negative  rate  of  excess  enthalpy 
generation  for  t  <  h,  then  Qexccss  >  AQ  during  the  period  of  the  calibration 
pulse,  ti  <  t  <  t2  [see  Footnote  C.50]. 

We  conclude  that  the  evaluation  given  by  NHE  is  invalid  and  that  it  is  likely  that  this 
evaluation  is  subject  to  several  distinct  errors. 

A.4  ICARUS  Type  Interpretation  of  the  Experiment 

As  a  first  step,  we  correct  the  (kR')n  spreadsheet.  Table  A.2,  by  including  the  magnitude 
of  the  calibration  pulse,  AQ,  in  the  definition  of  the  "lower  bound  heat  transfer 
coefficient".  The  modified  values  of  109  (kR')n  are  shown  in  column  15  while  the  values 
of  109  (k/)n  and  109  (k/l^  are  given  in  columns  16  and  17.  The  values  of  109 

(kR')i  i  in  the  region  ti  <  t  <  t2  can  now  be  shown  together  with  those  for  t  <  fi  and  t  >  t2 
on  a  graph  using  a  single  scale  for  the  ordinate,  Fig.  A.  10.  While  we  cannot  be  certain 
whether  or  not  an  incorrect  choice  of  CpM  can  explain  the  fall  of  (kR')n  in  the  region  t>t| 
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[but  close  to  this  time]  or  the  rise  for  t>t2  [but  close  to  this  time],  it  is  clear  that  (kR')n 
drops  markedly  in  the  region  t|<t<t2  compared  to  the  values  for  t<t|  and  t>t2.  Such  a 
drop  in  (kR')n  can  only  be  due  to  the  neglect  of  the  build  up  of  the  rate  of  excess  enthalpy 
generation  during  1 1 <t<t2 .  It  follows  that  the  increase  in  temperature  due  to  the 
calibration  pulse  increases  the  rate  of  excess  enthalpy  generation.  In  fact,  the  experiment 
shows  a  very  early  establishment  of  "positive  feedback"  as  is  indeed  evident  from  the  plot 
of  the  "raw  data",  Fig.  A. 5  [see  Footnote  C.51],  It  is  very  important  that  the  presence  of 
"positive  feedback"  can  be  established  by  a  simple  examination  of  a  (kR')n  spreadsheet 
constructed  according  to  the  instructions  in  the  ICARUS- 1  Handbook  (2). 

It  should  be  noted  that  the  amplitude  of  the  calibration  pulse  would  have  had  to  be  AQ  = 
0.2763  W  in  order  to  bring  the  values  of  (kR')i ,  in  the  region  t|<t<t2  to  the  level  of  the 
regression  line  which  applies  to  the  data  for  t<t,  and  t>t2.  Such  a  change  in  AQ  is  beyond 
all  possibilities. 

The  next  step  is  to  prepare  a  modified  (kR')n  spreadsheet  where  we  correct  the  enthalpy 
inputs  [see  (c),  (d)  and  (e)  in  Section  A. 3]  and  present  the  data  in  a  form  suitable  for  the 
application  of  equation  (B.12)  columns  3  and  4  of  Table  A. 3,  while  column  5  gives  the 
result  of  the  application  of  equation  (B.  12).  In  view  of  the  early  intervention  of  "positive 
feedback",  we  would  only  expect  to  be  able  to  apply  equation  (B.12)  at  times  close  to  ti 
where  we  see  that  the  "true  heat  transfer  coefficient",  (kR')n  must  be  at  least  0.83808  x 
10  WIC4  while  the  water  equivalent,  CpM  is  of  the  order  of  454  JK'1  [in  agreement  with 
the  "guesstimate",  see  Section  A.3], 

The  influence  of  "positive  feedback"  on  the  failure  of  simple  methods  for  the  evaluation 
of  the  "lower  bound"  and  "true  heat  transfer  coefficients"  as  well  as  of  the  water 
equivalent  of  the  cell  is  also  shown  clearly  by  attempts  to  derive  (kR‘)  m  [which  rely  on 
the  combination  of  data  for  the  time  regions  t|<t<t2  and  t2<T,  (see  Section  B)].  This 
evaluation  has  been  found  to  be  especially  useful  in  the  analyses  of  data  sets  for  "blank" 
experiments  [e.g.  see  ^  ^].  Fig.  A.  1 1  illustrates  that  we  are  unable  to  obtain  a  satisfactory 
interpretation  of  such  data  for  the  experiment  as  is  also  shown  by  column  6  of  Table  A.3. 

Column  7  of  Table  A.3  gives  the  values  of  109  (kR')n  while  column  8  gives  the  running 
mean,  10  (kR  ')u  ,  of  column  7.  We  can  see  that  these  values  do  not  differ  greatly  from 
the  corresponding  values  in  columns  9  and  15  of  Table  A. 2  which,  are  based  on  the 
correction  of  (e)  alone  [see  Section  A.3]  of  the  NHE  analysis.  Fig.  A.  1 2  gives  the  plot 
of  the  data  in  column  7  versus  time  and  also  shows  the  variation  of  109  (kR')ii  with  time 
predicted  using  the  values  for  t<t j  and  the  known  behavior  established  with  appropriate 
"blank"  experiments  e.g.  (4)’(l3)'(l4).  As  jn  the  case  of  the  data  in  Fig.  A.  10,  we  can  see 
that  the  temperature  rise  induced  by  the  calibration  pulse  leads  to  a  decrease  in  (kR')n 
while  the  cooling  consequent  on  the  termination  of  the  pulse  leads  to  an  increase  in 
(kR')n.  These  changes  can  only  be  due  respectively  to  an  increase  and  decrease  in  the 
rate  of  excess  enthalpy  generation,  which  cannot  be  taken  into  account  in  derivim*  the 
values  of  (kR')i  i,  i.e.  the  effects  of  "positive  feedback". 


Fig.  A.  12  shows  that  we  still  observe  discontinuities  in  the  "lower  bound  heat  transfer 
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coefficient,  (kR')ii"  at  t]  and  t2.  However,  it  is  evident  that  there  can  be  no  mechanism, 
which  could  account  for  such  changes,  which  must  therefore  be  due  to  an  error  in  the 
analysis.  The  most  obvious  error  is  the  use  of  an  incorrect  value  of  CpM  [see  (b).  Section 
A.3],  The  analysis  of  the  time  dependence  according  to  equation  (B.12)  in  the  region  t  > 
ti  [but  adjacent  to  tj]  indicates  that  the  correct  value  is  ~  450  JK'1,  and  columns  9,  10 
and  11  of  Table  A.3  give  respectively  CpM  (dAO/dt),  109  (kR')i  ]  and  109  (kR)u  based  on 
this  value  of  CpM.  Fig.  A.  13  shows  a  plot  of  the  values  in  column  10  versus  time,  and  we 
can  see  that  the  discontinuities  in  the  heat  transfer  coefficient  at  t>ti  and  t>t2  [but  adjacent 
to  these  times]  are  now  eliminated.  However,  as  expected,  the  effects  due  to  "positive 
feedback"  are  maintained. 

Column  12  of  Table  A.3  also  gives  values  of  the  rates  of  excess  enthalpy  generation 
calculated  with  the  "true  heat  transfer  coefficient,  (kR')i2"  =  0.85065  x  10‘9  WK-4 
established  using  the  (kR')2i  spreadsheets  and  other  confirmatory  evidence  outlined  in 
Section  A. 5  below.  Fig.  A.  14  shows  that  there  is  indeed  only  a  small  rate  of  excess 
enthalpy  generation  for  t<ti  while  the  application  of  the  calibration  pulse  leads  to  a  build¬ 
up  of  this  rate  which  again  decreases  for  t  >  t2  [there  is  a  small  long-term  increase  in  the 
rate  of  excess  enthalpy  generation  for  t>t2].  Fig.  A.  15  shows  a  similar  calculation  but 
using  the  NHE  methodology  [note  the  difference  in  scales  of  the  y-axes  in  Figs.  A.  14  and 
A.  1 5],  We  again  see  a  near  zero  rate  of  excess  enthalpy  generation  for  t<ti,  while  for  t>ti 
but  adjacent  to  ti  we  now  see  the  step  due  to  the  calibration  pulse,  AQ  =  0.2500  W.  In  the 
region  t]<t<t2,  we  then  see  the  build-up  in  the  rate  of  excess  enthalpy  generation  due  to 
"positive  feedback".  At  t  =  t2  but  adjacent  to  t2  we  again  see  a  step  in  the  total  observed 
rate  of  excess  enthalpy  generation.  As  expected,  this  step  again  corresponds  to  the 
expected  value  AQ  =  0.2500  W;  at  longer  times  we  see  the  gradual  decrease  of  the  rate  of 
excess  enthalpy  generation  due  to  the  removal  of  the  effects  of  "positive  feedback". 

The  comparison  of  the  interpretation  of  the  (kR')n  spreadsheet  prepared  according  to  the 
instructions  for  the  ICARUS  Systems  <1)(2)  with  that  prepared  using  the  NHE 
methodology.  Section  A.3,  illustrates  the  importance  of  following  the  instructions  laid 
down  in  the  Handbooks  (,)  (2).  We  can  also  see  that  a  great  deal  of  information  about  the 
behavior  of  the  systems  can  be  derived  from  the  interpretation  of  the  correctly  evaluated 
(kR')n  spreadsheets  and  these  spreadsheets  serve  as  the  basis  for  the  next  stage  of  the 
analysis  of  the  measurement  cycles. 

This  next,  third,  step  in  the  data  analysis  is  the  preparation  of  the  (kR*)2i  spreadsheet  for 
the  whole  measurement  cycle  shown  in  Table  A. 4  which  can  be  used  equivalently  to 
produce  the  (kR’)3,  spreadsheet  [Table  A.4  is  an  abstract  of  the  two  spreadsheets]  [see 
Footnote  C.52],  In  evaluating  the  integrals  used  in  the  preparation  of  this  spreadsheet,  it 
has  been  noted  that  ti  and  t2  correspond  closely  to  the  relevant  measurement  intervals 
[this  is  not  illustrated  in  this  Report],  The  integrals  have  been  evaluated  using  the 
trapezium  rule  following  the  insertion  of  additional  data  points  t|  and  t2  [see  Footnote 
C  53] 

It  can  be  seen  that  Table  A.4  and  a  comparison  of  the  plots  of  (kR')2J  and  (kR')3i  versus 
time.  Fig.  A  16,  with  the  corresponding  plots  for  "blank  experiments",  e.g.  see  (4)(S)(I4), 
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shows  very  clearly  the  intervention  of  "positive  feedback"  due  to  the  superposition  of  the 
calibration  pulse.  If  we  focus  attention  first  of  all  on  the  behavior  of  (kR*)3i  for  t<fi,  then 
we  see  the  expected  small  decrease  with  increasing  time  [see  Footnote  C.54],  For  t>t|  we 
see  a  more  rapid  decrease  due  to  the  onset  of  "positive  feedback".  The  effects  of  this 
"positive  feedback"  decrease  for  t>t2,  so  that  we  observe  a  small  increase  of  (kR')3i  with 
increasing  time  in  this  region. 

The  variation  of  (kR')2i  with  time  can  be  interpreted  in  a  similar  way  provided  one  bears 
in  mind  that  there  is  now  no  region  in  lime  in  which  the  integrals  used  in  the  calculation 
of  the  heat  transfer  coefficient  are  independent  of  the  effects  of  "positive  feedback". 
This  influence  of  "positive  feedback"  on  the  integrals  used  in  the  evaluation  of  (kR')2i 
explains  why  we  cannot  obtain  a  satisfactory  evaluation  of  the  target  value  of  the  "lower 
bound  heat  transfer  coefficient,  (1cr')26i”  column  5  of  Table  A.4.  We  would  only  expect 
to  be  able  to  apply  the  ICARUS  methodology  in  a  region  of  time  where  the  influence  of 
"positive  feedback"  can  be  expected  to  be  adequately  small,  say  in  the  region  72,300  - 
75,300  s  of  the  measurement  cy  cle.  The  estimates  of  the  "lower  bound  heat  transfer 
coefficient,  (kR')26i"  and  of  CpM  are  shown  in  bold  type  in  column  5  of  Table  A. 4. 

Again,  if  we  bear  in  mind  the  influence  of  "positive  feedback"  on  the  integrals,  we  would 
only  expect  to  be  able  to  apply  the  ICARUS  methodology  to  the  initial  part  of  the 
calibration  period  to  give  estimates  of  the  target  value  (kR')36i  [say  between  32,400  and 
35,400  s  in  Table  A. 4].  These  estimates  of  the  "lower  bound  heat  transfer  coefficient, 
(kR')36i"  and  of  CpM  are  shown  in  bold  type  in  Column  10  of  Table  A. 4. 

It  is  important  here  to  draw  attention  to  a  restriction,  which  we  have  applied  to  the 
evaluation  of  the  integral  heat  transfer  coefficients.  Table  A.4  shows  that  abscissa  in 
columns  2  or  7  are  always  numerically  small  compared  to  the  ordinates  in  columns  3  or 
8.  In  consequence  the  heat  transfer  coefficients  (kR')2i  and  (kR')3i  in  columns  4  or  9  can 
be  evaluated  with  good  precision  notwithstanding  the  limited  accuracy  of  the  values  of 
the  water  equivalents,  CpM. 

This  effect  becomes  much  more  important  when  we  come  to  consider  the  evaluation  of 
the  "true  heat  transfer  coefficients,  (kR')262  and  (kR')362  [see  the  discussion  of  Table  A.5 
further  below].  For  the  particular  case  of  the  experiment  we  see  that  we  must  rely  on  the 
evaluation  of  the  "true  heat  transfer  coefficient,  (kR')362  in  view  of  the  intervention  of 
"positive  feedback"  even  though  this  heat  transfer  coefficient  is  not  the  target  value  of  the 
ICARUS  methodology  [see  Footnote  C.55],  Columns  1 1  and  12  of  Table  A.4  now  show 
that  the  abscissae  and  ordinates  have  comparable  magnitudes  in  the  time  range  in  which 
we  might  conceivably  be  able  to  neglect  the  effects  of  "positive  feedback"  and  the 
evaluation  of  (kR')362  therefore  fails 

The  comments  which  have  been  made  about  the  evaluation  of  the  integral  heat  transfer 
coefficients  using  the  whole  measurement  cycles  apply  equally  to  the  evaluations 
according  to  the  instructions  and  software  in  the  ICARUS  Systems  (l)<2).  Table  A.5.  The 
evaluation  of  (kR')26i  [not  show  n  in  Table  A.5]  fails  because  of  the  intervention  of 
"positive  feedback"  and  the  precision  of  (kR')3i  and  (kR')36i,  columns  4  and  5  of  Table  A.5 
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is  low  because  of  the  need  to  restrict  attention  to  the  region  t>ti  but  close  to  ti.  This  is 
equally  true  of  the  accuracy  of  (kR')32  and  (kR')362  shown  in  columns  8  and  9  of  Table 
A.5.  However,  the  evaluations  of  these  coefficients  are  instructive  because  it  is  virtually 
certain  that  the  value  of  the  "true  heat  transfer  coefficient"  quoted  by  NHE  is  either  the 
value  of  (kR')32  at  a  particular  time  or  else  (kR’)362  evaluated  over  a  particular  range  of 

time. 

Columns  10-13  of  Table  A. 5  illustrate  the  evaluations  of  (kR')22and  (kR')262  and  it  is 
clear  that  we  cannot  apply  the  ICARUS  Methodology  as  set  out  in  the  Handbooks  in 
view  of  the  intervention  of  "positive  feedback".  We  therefore  have  to  investigate  whether 
we  can  modify  the  approach  so  as  to  allow  the  determination  of  the  "true  heat  transfer 
coefficients".  We  have  to  note  that  it  is  unlikely  that  we  would  be  able  to  find  a  generally 
valid  procedure  because  it  is  in  general  not  possible  to  calibrate  closed  loop  systems 
subject  to  "positive  feedback".  However,  for  the  particular  example  of  Day  3  of  the 
experiment  we  can  see  that  the  effects  of  "positive  feedback  are  relatively  small  and, 
moreover,  confined  in  the  time-domain,  Fig.A.14.  We  can  therefore  include  the  observed 
values  of  the  rates  of  excess  enthalpy  generation  in  the  evaluation  of  the  integral  of  the 
enthalpy  input  and  use  this  modified  integral  to  re-evaluate  (kR')22and  (kR')262  columns  14 
-  16  of  Table  A. 5.  Fig.  A.  17  illustrates  this  evaluation.  It  can  be  seen  that  we  do  indeed 
now  obtain  a  satisfactory  fit  to  equation  (B.24),  which  explains  the  choice  of 

(kR’)262  =  0.85065  x  1  O'9  WK4  and  CpM  =450JK-' 
for  the  further  evaluation  of  the  data. 

In  view  of  the  fact  that  this  evaluation  of  the  "true  heat  transfer  coefficient, 
(kR')262  requires  the  development  of  a  special  approach,  it  is  necessary  [and  advisable]  to 
investigate  whether  the  value  obtained  can  be  confirmed  by  other  means  using  different 
parts  of  the  experiment  [i.e.  other  measurement  cycles].  Such  confirmations  can  be 
obtained  using  the  measurements  on  Day  61  and  the  first  57  hours  of  Days  1  and  2. 
These  confirmations  are  outlined  in  Sections  A. 5  and  A. 6  respectively. 

A.5  Application  of  the  ICARUS  Type  Interpretation  to  the  Data  for  Day  61 

It  has  been  shown  in  Sections  A.3  and  A.4  that  the  early  intervention  of  "positive 
feedback"  requires  us  to  modify  the  ICARUS  evaluation  strategies  in  order  to  achieve  the 
calibration  of  the  system  i.e.  to  determine  the  value  of  the  "true  heat  transfer  coefficient". 
It  is  therefore  important  to  find  confirmatory  evidence  that  this  heat  transfer  coefficient  is 
indeed  equal  to  0.85065  x  10'9  WK'4  as  given  at  the  end  of  the  previous  Section. 
Evidence  pertinent  to  this  conclusion  is  presented  in  the  present  Section  as  well  as 
Section  A. 6. 

We  note  in  the  first  place  the  values  of  the  total  excess  enthalpy  for  each  day  of  operation 
calculated  using  the  "true  heat  transfer  coefficient,  (kR')32"  =0.79350x10  WK  as 
given  by  the  NHE  evaluation  as  well  as  those  calculated  with  "true  heat  transfer 
coefficient  (kR')262,r  =  0.85065  x  10'9  WK'4  as  determined  in  Section  A.4  using  the 
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ICARUS  methodology.  Table  1.  These  values  are  plotted  in  Figs.  A.  18  and  A.  19 
respectively.  We  can  see  immediately  that  the  evaluation  given  by  NHE  must  be 
incorrect  because  we  obtain  negative  excess  enthalpies  for  some  of  these  days,  which 
contravenes  the  Second  Law  of  Thermodynamics  [c.f.  Section  A. 3].  On  the  other  hand 
the  evaluation  based  on  the  heat  transfer  coefficient  given  by  the  modified  ICARUS 
evaluation  scheme  only  gives  a  very  slightly  negative  excess  enthalpy  for  Day  61. 

It  is  therefore  reasonable  to  assume  that  the  rate  of  excess  enthalpy  generation  on  Day  61 
is  close  to  zero.  The  evaluation  of  the  "lower  bound  heat  transfer  coefficient,  (kR')n" 
from  the  relevant  (kR')n  spreadsheet.  Table  A. 6,  must  therefore  be  close  to  the  values  of 
the  "true  heat  transfer  coefficient,  (kR')i2".  Fig.  A. 20  gives  a  plot  of  the  relevant  data 
compared  to  the  plot  which  we  predict  using  the  value  (kRr)i2  =  0.85065  x  10'9  WK4  and 
the  variation  of  (kR')n  with  time  given  by  the  relevant  "blank"  experiments  (4)(13f(14)  u 
can  be  seen  that  the  observed  values  of  (kR')n  are  in  close  accord  with  those  which  we 
would  predict  on  the  assumption  that  there  is  only  a  low  rate  of  excess  enthalpy 
generation  on  that  day. 

It  can  be  seen  that  there  is  only  one  region  of  time  in  which  there  is  a  marked  deviation 
from  the  predicted  behavior  namely  for  0  <  t  <  ~  10,000  s  on  Day  61.  In  this  region 
(kR')n  is  markedly  larger  than  the  expected  value  and,  moreover,  decreases  rapidly  with 
time  to  these  predicted  values.  It  is  possible  therefore,  that  the  deviation  seen  in  this 
time  range  can  be  attributed  to  the  "overfilling"  of  the  cell  [see  Footnote  C.56],  Separate 
measurements  have  shown  that  the  pseudo  radiative  heat  transfer  coefficient  increases  by 
~  5%  over  the  expected  value  presumably  because  of  an  increase  in  the  conductive 
contribution  through  the  top  of  the  cell.  Another  explanation  already  noted  in  Section 
A.2  is  the  sudden  cooling  of  the  cell  by  the  large  addition  of  D2O  on  Day  61. 

A.6  A  Pre-ICARUS  Evaluation  of  the  True  Heat  Transfer  Coefficient 

It  is  possible  to  find  a  further  value  of  the  “true  heat  transfer  coefficient  (kR')i2  by 
applying  a  method  used  in  1992  (15),(16).  it  was  shown  at  that  time  that  the  "lower  bound 
heat  transfer  coefficient,  (kR')n”  decreases  markedly  from  the  expected  value  during  the 
initial  stages  of  the  measurement  cycles.  In  this  case  the  decrease  is  due  to  the 
exothermic  absorption  of  deuterium  in  the  lattice.  It  would  be  expected,  therefore,  that 
the  “lower  bound  heat  transfer  coefficient,  (kR')n”  would  rise  markedly  to  the  expected 
"true"  value  as  this  process  is  completed  with  the  proviso  that  we  can  observe  a  period  of 
operation  during  which  there  is  zero  excess  enthalpy  generation.  It  follows  that  we  can 
derive  a  value  of  the  "true  heat  transfer  coefficient,  (kR*)i2”  from  the  maximum  of  the 
"lower  bound  heat  transfer  coefficient,  (kR')n”  which  is  observed  w'ith  increasing  time. 

Fig.  A.21  shows  the  relevant  data  for  the  first  57  hours  of  operation  of  the  experiment 
[i.e.  up  to  the  time  of  application  of  the  heater  calibration  pulse  on  Day  3],  The  full  line 
at  the  top  shows  the  expected  variation  of  (kR’)n  with  time  based  on  the  value  of  (kRf)i2 
at  t  =  t2  on  Day  3  [i.e.  <kR')i2=  0.85065  x  10 9  WK4],  the  assumption  of  zero  excess 
enthalpy  generation  [i.e  (kR')n  =  (kR*)i 2]  and  the  known  variation  of  (kR')n  w'ith  time 
established  with  "blank"  experiments  (4)(1?)  (14)  Jt  can  be  seen  that  (kR’)n  does  indeed 
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rise  to  the  predicted  levels  as  the  charging  of  the  electrode  is  completed. 

Fig  A  21  also  shows  the  derived  rates  of  excess  enthalpy  generation  based  on  the 
experimental  values  of  (kR')n  and  the  assumption  that  the  "true  heat  transfer  coefficient" 
is  given  by  the  regression  line.  The  horizontal  line  here  gives  the  expected  value  of  the 
rate  of  excess  enthalpy  generation  based  on  the  assumption  that  the  current  efficiency  for 
the  charging  of  the  electrode  is  100%  and  that  the  enthalpy  of  absorption  of  deuterium  in 
the  lattice  is  40  kJ  Mole1.  It  can  be  seen  that  the  experimental  values  are  in  reasonable 
accord  with  this  prediction. 

Fie  A  21  furthermore  shows  that  there  is  a  small  build  up  of  excess  enthalpy  generation 
on  Day  3  following  the  completion  of  the  charging  process  [compare  J. 

The  data  shown  in  Fig.  A.21  also  lead  to  a  number  of  important  questions,  which  cannot 
be  resolved  at  the  present  time  [see  Footnote  C.57]. 

A.7  Day  68:  The  Period  0  <  t  <  21,300  s  During  Which  the  Cell  is  Driven  to 
Dryness 

We  consider  next  the  penultimate  day  of  the  investigation  of  the  experiment;  the  cell  is 
driven  to  dryness  during  the  first  part  of  this  measurement  cycle.  The  "raw  data"  for  this 
section  is  given  in  columns  1-5  of  Table  A.7  while  columns  6-9  summarize  the 
interpretation  in  the  relevant  NHE  spreadsheet  [the  part  using  the  temperature 
measurements  with  the  "short  thermistor"]. 

We  can  draw  a  number  of  important  conclusions  from  the  "raw  data"  alone.  We  note  in 
the  first  place  that  the  temperature  given  by  the  "long  thermistor"  is  now  slightly  higher 
than  that  given  by  the  "short  thermistor"  whereas  the  opposite  is  true  for  measurements 
made  at  low  temperatures.  At  first  sight  such  a  change  might  be  attributed  to  a  genuine 
effect  namely,  the  increase  in  the  enthalpy  input  in  the  bottom  part  of  the  cell  [containing 
the  Pd-B  cathode].  However,  such  an  interpretation  is  unlikely  because  the  temperature 
difference  between  the  two  thermistors  is  essentially  constant  for  say  20,000  s  even 
though  the  enthalpy  input  increases  by  a  factor  of  three.  It  is  more  likely  therefore  that 
this  particular  temperature  difference  is  a  further  manifestation  of  errors  in  the  calibration 
of  the  thermistors.  Experimentally,  furious  boiling  and  swirling  actions  were  centered 
around  the  cathode  during  this  period  which  may  have  affected  the  temperature  readings 
of  the  longer  thermistor.  There  was  also  about  2  cm  of  foam  on  the  liquid  interface. 

The  temperature  differences  between  the  two  thermistors  are  appreciably  larger  for  the 
last  four  data  acquisition  points  and  this  difference  is  especially  marked  for  the  last  point 
[0.590  K],  Such  a  difference  is  to  be  expected  if  the  "long  thermistor"  is  now  in  the 
relatively  concentrated  LiOD  solution  or  foam  residue  while  the  "short  thermistor"  is  m 
the  vapor  phase  However,  we  also  have  to  note  that  the  temperature  at  both  positions  is 
above  that  of  the  boiling  point  of  pure  D20  Evidently,  we  have  to  take  into  account  the 
increase  of  the  boiling  point  with  the  electrolyte  concentration  as  the  D20  is 
progressively  evaporated  [see  Section  B.2],  However,  we  also  have  to  take  note  of  the 
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fact  that  the  vapor  phase  can  be  superheated  [albeit  to  only  a  limited  extent,  see  Footnote 
C.58].  If  we  do  not  take  account  of  the  increase  of  the  boiling  point  with  concentration, 
we  arrive  at  the  impossible  result  of  negative  enthalpies  of  evaporation  with  increasing 
temperature  as  shown  by  the  NHE  evaluation.  We  also  have  to  use  the  correct 
atmospheric  pressure  in  the  calculation  of  the  rate  of  evaporative  cooling,  and  we  need  to 
change  the  thermoneutral  potential  and  the  water  equivalent  of  the  cell  in  the  NHE 
evaluation.  As  the  water  equivalent  of  the  cell  only  leads  to  a  significant  term  CpM 
(dA0/dt)  in  the  initial  stages  for  Day  68,  it  has  been  assumed  that  CpM  is  unchanged 

throughout  the  stage  leading  to  evaporation  to  dryness  [however,  see  further  comments  in 
Section  A. 8], 

The  modified  spreadsheet  is  shown  in  Table  A. 8,  which  is  still  an  intermediate  step  in  the 
evaluation  because  of  the  neglect  of  reflux  in  the  cell.  However,  this  intermediate 
spreadsheet  is  sufficient  for  the  next  stage  of  the  evaluation.  It  has  already  been  pointed 
out  in  Section  B.2  that  it  is  desirable  to  carry  out  such  a  first  stage  of  the  evaluation 
independently  of  the  complexities  of  the  evaluation  of  the  rate  of  evaporative  cooling  and 
reflux  in  the  cell.  We  can  carry  out  such  a  first  step  in  a  straightforward  way  by 
calculating  the  total  enthalpy  input  and  output  from  the  cell  using  columns  1-3  of  Table  9 
and  the  known  "true  heat  transfer  coefficient".  We  obtain: 

enthalpy  input  =378,105  J;  CPMA6  =  2,012  J 

enthalpy  output  =  212,686  J  for  (kR’)12  =  0.85065  x  10'9  WK"4 

=  223,320  J  for  (kR’)i2  =  0.893 18  x  10 9  WK’4 

The  enthalpy  available  for  evaporation  of  the  D20  is  therefore  163,407  J  if  the  lower 
value  of  the  heat  transfer  coefficient  applies  or  152,773  J  if  the  higher  value  applies. 
However,  we  require  207,031  J  for  the  evaporation  of  the  4.968  Moles  of  D20  in  the  cell. 
There  is  therefore  a  short  fall  of  either  43,624  or  54,258  J,  a  short  fall,  which  can  only  be 
supplied  by  excess  enthalpy  generation  in  the  cell.  If  this  excess  enthalpy  is  supplied  at  a 
uniform  rate  throughout  the  period  leading  up  to  evaporation  to  dryness,  then  this  would 
require  a  mean  excess  rate  of  2.048  or  2.547  W  depending  on  whether  the  "true  heat 
transfer  coefficient"  has  the  lower  or  higher  value. 

This  calculation  is  similar  to  one,  which  has  been  described  previously  (8)  except  that  the 
published  version  included  comments  on  the  time  dependence  of  the  rate  of  excess 
enthalpy  generation.  It  is  quite  obvious  that  the  rate  of  excess  enthalpy  generation  must 
increase  with  time  because  the  initial  rate  on  Day  68  is  less  than  1  W.  It  is  important 
therefore  to  try  to  establish  the  variation  of  the  rate  of  excess  enthalpy  generation  with 
time,  if  only  to  make  a  connection  with  the  initial  rate  of  "Heat  after  Death"  observed 
after  the  cell  has  reached  dryness  [see  Section  A. 8],  In  order  to  derive  this  variation,  we 
have  to  include  an  estimate  of  the  rate  of  reflux  in  the  cell,  and  this  part  of  the  calculation 
will  follow  the  scheme  outlined  in  Section  B.2  [see  Footnote  C.59],  The  details  of  the 
evaluation  are  given  in  Table  A. 9.  We  can  see  that  the  negative  values  of  the  enthalpies 
seen  in  Tables  A  7  and  A.8  are  now  eliminated  as  the  D20  in  the  cell  is  maintained  by  the 
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amount  of  reflux.  The  total  amount  evaporated  is  also  in  reasonable  accord  with  the 
amount  of  D20  initially  in  the  cell.  It  is  important  to  realize  that  we  have  assumed  that 
the  whole  of  the  heat  transfer  from  the  cell  in  the  region  filled  with  vapor  leads  to  re¬ 
condensation  [see  Section  B.2]  i.e.  we  have  overestimated  the  reflux  and  underestimated 
the  amount  evaporated.  We  should  also  note  that  the  calculation  is  improved  somewhat 
if  we  allow  for  the  fact  that  the  boiling  point  reaches  a  limit  due  to  the  limited  solubility 
of  LiOD  in  D20  at  the  boiling  point  [this  aspect  is  not  illustrated  in  this  Report]. 

Although  the  calculation  as  outlined  gives  a  reasonable  interpretation  of  the  behavior  of 
the  cell  as  the  contents  are  driven  to  dryness  [elimination  of  negative  enthalpies  of 
evaporation],  we  nevertheless  still  derive  negative  rates  of  excess  enthalpy  generation  at 
long  times.  This  is  undoubtedly  due  to  remaining  inaccuracies  in  the  calculation  of  the 
rates  of  evaporative  cooling.  At  the  present  time  it  is  best  to  restrict  attention  to  the 
earlier  part  of  the  period  leading  to  evaporation  to  dryness,  say  to  t  <  18,000  s.  We  see 
from  Table  A. 9  that  the  rate  of  excess  enthalpy  generation  reaches  ~  9.3  W  at  this  time, 
or  27  W  cm'3.  It  is  important  to  realize  that  similar  orders  of  magnitude  are  obtained 
even  with  the  interpretation  given  by  NHE,  i.e.  the  estimate  is  "robust"  [see  Footnotes 
C.60-62], 


A.8  Day  68:  The  Period  21,300  s  <  t  <  86,400  s  Following  Evaporation  to 
Dryness 

As  has  been  noted  in  Section  A.l  one  of  the  objectives  of  the  present  investigation  has 
been  the  search  for  the  presence  [or  absence]  of  the  effects  of  "Heat  after  Death".  The 
period  following  the  evaporation  to  dryness  on  day  68  is  an  example  of  the  protocol 
originally  described  as  Case  (5) (9),(l0)  i.e.,  cell  empty:  cell  allowed  to  boil  dry;  cell  then 
maintained  at  the  rail  voltage  of  the  galvanostat. 

The  original  investigation  was  divided  into  two  parts: 

(i)  the  investigation  and  interpretation  of  the  cooling  curves  following 
evaporation  to  dryness; 

(ii)  the  evaluation  of  thermal  balances  in  the  corresponding  period. 

Attention  here  will  be  confined  to  the  second  of  these  approaches. 

Table  A.  10  gives  the  "raw  data"  for  the  first  9000  s  of  operation  during  this  period 
[columns  1-4],  The  data  set  has  been  restricted  to  the  first  30  acquisition  intervals  so  as 
to  avoid  the  complications  due  to  the  heater  calibration  pulse,  which  was  still  applied  on 
Days  68  and  69.  Columns  5  and  6  gives  the  data  required  for  the  calculation  of  the  rates 
of  excess  enthalpy  generation  shown  in  columns  7  and  8.  It  has  been  assumed  that  the 
water  equivalent  of  the  empty  cell  is  28.3  JK  ',  the  value  for  the  cell  filled  with 
electrolyte  less  the  contribution  due  to  the  5  moles  of  D20  removed  during  evaporation  to 
dryness. 

The  values  of  the  rates  of  excess  enthalpy  generation  shown  in  column  7  have  been  based 
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on  the  “true  heat  transfer  coefficient,  (kR')|2”  observed  for  the  cell  filled  with  electrolyte, 
i.e.  0.85065  x  10'9  WK'4,  which  will  certainly  apply  to  initial  stage  of  the  observation  of 
"Heat  after  Death"  when  the  cell  is  filled  mainly  with  D20  vapor.  However,  calibrations 
of  cells  filled  with  air  (9)  (l0)  have  shown  that  the  heat  transfer  coefficient  falls  to  about 
0.75  of  the  value  for  the  cells  filled  with  electrolyte.  The  values  of  the  rates  of  excess 
enthalpy  generation  shown  in  column  8  have  therefore  been  calculated  using  (kR*)i  2  = 
0.65  x  10'9  WK'4. 

As  would  be  expected,  the  rates  of  excess  enthalpy  generation  shown  in  column  8  are 
lower  than  those  in  column  7,  but  the  decreases  in  both  cases  are  broadly  in  line  with 
those  expected  for  a  process  controlled  by  a  diffusional  relaxation  time.  The  initial  rate 
of  excess  enthalpy  generation  given  by  column  7  is  approximately  the  same  as  the  rate 
reached  during  the  period  0  <  t  <  21,300  s  as  the  cell  is  being  driven  to  dryness,  Fig. 
A. 22.  Such  a  correspondence  would,  of  course,  be  expected  if  excess  enthalpy 
generation  takes  place  in  the  bulk  of  the  metal  phase. 

We  note  also  that  the  rate  of  excess  enthalpy  generation  is  about  10  times  that  of  the  rate 
of  enthalpy  input  during  this  period  of  "Heat  after  Death". 

A.9  Day  69:  The  Period  2400  s  <  t  <  32,400  s 

This  period  is  of  special  interest  in  the  operation  of  the  cell  because  the  cell  was 
disconnected  from  the  galvanostat  at  2400  s  so  that  the  enthalpy  input  was  zero  during 
the  remaining  period  of  operation.  In  any  search  for  the  effects  of  "Heat  after  Death"  the 
protocol  there  should  be  described  as  Case  (6)  (9)’(10)  i.e.  cell  empty:  cell  allowed  to  boil 
dry;  cell  disconnected  from  the  galvanostat  and  that  the  application  of  Case  (6)  was 
preceded  by  a  period  covered  by  Case  (5)  as  described  in  Section  A. 8. 

The  examination  of  the  behavior  of  the  cell  has  been  restricted  here  to  the  time  t  <32,400 
s  as  the  usual  calibration  pulse  was  applied  at  t|  =  32,400  s.  The  Joule  heat  injected  by 
the  calibration  system  is  developed  in  a  small  volume  so  that  this  calibration  cannot  be 
used  to  derive  the  "true  heat  transfer  coefficient"  of  the  cell  for  the  particular  operating 
conditions  [see  Footnote  C.63], 

The  "raw  data"  for  this  period  of  operation  are  given  in  Table  A.l  1  and  a  "cooling  curve" 
is  plotted  in  Fig.  A. 23.  It  can  be  seen  that  although  the  temperature  differences  between 
the  cell  and  water  bath  are  small,  they  are  nevertheless  significant.  Inspection  of  Fig. 
A. 23  shows  that  there  must  be  a  source  of  enthalpy  in  the  system  firstly,  because  the  rate 
of  cooling  at  short  times  is  too  slow  to  be  accounted  for  by  the  cooling  of  a  calorimeter 
with  a  water  equivalent  of  28.3  JK'1  and  any  conceivable  value  of  the  heat  transfer 
coefficient;  secondly,  because  we  can  detect  at  least  one  period  during  which  the  cell 
contents  are  reheated! 

The  data  in  Table  A.l  1  are  suitable  for  the  analysis  of  the  cooling  curve  according  to  the 
method  originally  outlined  <9),(l0)  using  the  equation 
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ln[yo(l+y)/y(l+yo)]  +  tan"'(l+y)  -  tan''(l+y0) 

=  4(kR’)03ba,h  t/CpM 

(A.3) 

y  =  A0/0bath  andyo  =  A0o/0bath 

(A-4) 

and  A0o  is  the  initial  temperature  difference. 

Fig.  A. 24  shows  a  plot  of  the  experimental  data  according  to  equation  (A. 3);  the  full  line 
shows  the  predicted  behavior  using  CPM  =  28.3  JK'1  and  (1<r')  =  0.65  x  10  9  WK  4.  The 
deviations  from  this  plot  due  to  enthalpy  generation  are  similar  to  those  previously 
observed (9)/l0). 

We  can  also  make  thermal  balances  at  each  point  of  the  cooling  curves  using  particular 
values  of  the  water  equivalent  and  "true  heat  transfer  coefficient".  Those  based  on  CPM 
=  28.3  JK'1  and  (kR')i2  =  0.65  x  10'9  WK'4  give  initial  rates  of  enthalpy  generation  ~0.5 
W.  Unfortunately,  the  thermal  balances  in  the  period  preceding  the  disconnection  of  the 
cell  from  the  galvanostat  [i.e.  the  last  part  of  Case  5,  Section  A.8]  cannot  be  made  with 
sufficient  accuracy  to  allow  a  comparison  of  the  rates  of  enthalpy  generation  at  the  end  of 
the  period  following  Case  5  and  the  beginning  of  the  period  following  Case  6  [see 
comparison  of  the  rates  at  the  end  of  the  period  leading  to  evaporation  to  dryness  and  the 
beginning  of  the  period  following  Case  5,  Sections  A.7  and  A.8]. 

A.10  Days  25  and  26:  The  Period  Day  25  +  76,300s  <  t  <  Day  26  +  22,300s 

As  has  already  been  noted  in  Section  A.2,  there  were  frequent  changes  of  current  density 
in  the  experiment.  Consideration  of  Case  1  of  the  conditions  likely  to  give  demonstrations 
of  the  phenomenon  of  "Heat  after  Death" (9)<l0): 

Cell  full:  cell  operated  at  intermediate  temperatures;  cell  current  then 
reduced  in  stages  [Case  1  ] 

shows  that  the  change  of  current  close  to  the  start  of  Day  26  of  the  measurement  cycles  is 
likely  to  provide  the  best  example  of  this  particular  case,  see  Fig.  A.2.  There  are  two 
principal  reasons,  which  indicate  this.  In  the  first  place,  the  current  density  at  the  end  of 
Day  25  is  above  the  threshold  value  required  for  the  observation  of  the  phenomenon  (l2) 
while  on  Day  26  it  is  below  this  threshold  value.  Secondly,  the  cell  temperature  on  Day 
25  is  above  that  which  has  been  observed  to  be  important  for  the  onset  of  "positive 
feedback  (8)'(l5),(l6)  whereas  on  Day  26  it  drops  below  this  value.  We  would  therefore 
expect  to  see  a  marked  decrease  of  the  rate  of  excess  enthalpy  generation  at  the  start  of 
Day  26  from  the  value,  which  applies  at  the  end  of  Day  25  to  that  which  applies  towards 
the  end  of  Day  26  [see  Footnote  C.64], 

The  data  covering  measurements  in  the  last  stages  of  Day  25  and  the  beginning  of  Day  26 
are  given  in  Table  A.  12.  Here  columns  1-4  give  the  "raw  data",  columns  5  -  8  give  the 
derived  quantities  used  to  define  the  "lower  bound  heat  transfer  coefficient,  (kR')n" 
shown  in  column  9.  We  note  that  we  used  the  value  of  CPM  =  475  JK  1  in  view  of  the 


23 


overfilling  of  the  cell  on  Day  26  [see  Section  A. 2].  The  rates  of  excess  enthalpy 
generation  derived  are  plotted  in  Fig.  A. 25.  We  can  see  the  well-defined  fall  in  the  rate 
of  excess  enthalpy  generation,  which,  as  in  the  other  examples  of  Heat  after  Death 
discussed  in  this  Report,  is  consistent  with  a  diffusional  relaxation  time. 

Columns  1 1  and  12  of  Table  A. 12  also  give  the  "lower  bound  heat  transfer  coefficient, 
(kR')n"  and  rate  of  excess  enthalpy  generation  taken  from  the  spreadsheets  for  the 
evaluations  carried  out  by  NHE.  We  can  see  both  from  column  13  of  this  Table  and  from 
the  plot  in  Fig.  A. 26,  that  this  evaluation  predicts  a  negative  rate  of  excess  enthalpy 
generation  on  Day  25.  As  we  have  noted  elsewhere  in  this  Report,  such  negative  rates 
violate  the  Second  Law  of  Thermodynamics  and  are  certainly  due  to  the  use  of  the 
incorrect  value  of  the  "true  heat  transfer  coefficient,  (kR'))2"  given  by  the  NHE  analysis. 
Nevertheless,  we  can  see  from  Table  13,  or  Fig.  A.26,  that  we  can  detect  the  effects  of 
"Heat  after  Death"  on  Day  26  even  when  using  this  faulty  analysis.  Furthermore,  the 
increasing  values  of  the  "lower  bound  heat  transfer  coefficient  (kR')n"  on  that  day 
demonstrate  the  presence  of  a  rate  of  excess  enthalpy  generation,  which  decreases  with 
time. 

It  is  important,  however,  to  draw  attention  to  a  remaining  difficulty  in  the  interpretation, 
namely,  that  the  initial  rate  of  excess  enthalpy  generation  on  Day  26  is  larger  than  the 
final  rate  on  Day  25.  The  discrepancy  would  be  diminished  if  the  water  equivalent  were 
even  higher  than  475  JK'1  or  if  we  increased  (kR’)i2  in  view  of  the  evident  increase  in  the 
D20  content  of  the  cell.  Fig.  A.4.  A  higher  water  equivalent  is  not  likely  based  on  the 
experimental  measurements  of  the  electrolyte  volume.  It  does  not  seem  possible, 
therefore,  to  eliminate  the  effect  by  any  sensible  choice  of  the  values  of  CPM  and  (kR')i2 
thus  the  effect  is  likely  to  be  real.  If  this  is  so,  then  the  establishment  of  "Heat  after 
Death"  and/or  "positive  feedback"  would  be  more  complicated  than  is  indicated  by  the 
state  variables  alone.  For  example,  the  time  derivatives  may  also  be  involved  It  is 
evident  that  much  further  work  is  required  on  these  particular  aspects.  This  work  would 
be  justified  not  only  from  the  objective  of  clarifying  the  science  involved  but,  also, 
because  the  judicious  use  of  "positive  feedback"  and  "Heat  after  Death"  offers  us  the 
prospect  both  of  increasing  the  power  density  and  at  the  same  time,  of  increasing  the 
energy  efficiency.  It  should  be  noted  that  if  we  exclude  the  enthalpy  input  due  to  the 
cooling  of  the  cell,  the  rate  of  excess  enthalpy  generation  in  the  initial  stages  of  Day  26  is 
approximately  equal  to  the  enthalpy  input,  i.e.  a  power  gain  of  ~  100%  whereas  it 
approaches  ~  1000%  for  the  initial  stages  of  Heat  after  Death  according  to  Case  5  (9)’<l0). 
Section  A.8,  and  infinity  %  for  the  example  of  Case  6,  Section  A.9.  It  appears  that  if  the 
cooling  of  such  cells  is  prevented  [effectively  by  raising  the  temperature  of  the  heat  sink], 
then  enthalpy  generation  may  be  maintained  for  long  durations  [~1  week]  at  very  high 
energy  efficiencies  (26).  It  is  evident  that  this  aspect  of  the  work  requires  extensive  further 
investigation  particularly  with  regard  to  attempts  to  raise  the  power  density  of  such 
devices  while  maintaining  the  high  energy  efficiency. 

A.l  1  Further  Comments  and  Conclusions 

This  Pd-B  experiment  exhibits  all  the  key  features,  which  have  been  found  in  earlier 
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work.  These  are  in  the  main: 


(i)  excess  enthalpy  generation  in  the  early  stages  [t  <  2  Days]  due  to 
absorption  of  deuterium  in  the  lattice  followed  by 

(ii)  a  build  up  of  the  rate  of  excess  enthalpy  generation 

(iii)  the  development  of  "positive  feedback"  i.e.  the  increase  in  the  rate  of 
excess  enthalpy  generation  with  increase  of  temperature 

(iv)  a  marked  increase  in  the  rate  of  excess  enthalpy  generation  at  temperatures 
close  to  the  boiling  point  of  the  electrolyte 

(v)  a  variety  of  examples  of  the  phenomenon  of  "Heat-after-Death"  i.e.  a 
maintenance  of  elevated  rates  of  excess  enthalpy  production  following 
reduction  of  the  current  density  which  may  be  accompanied  by  the 
complete  evaporation  of  the  electrolyte. 

At  the  same  time  there  are  some  marked  differences  between  this  experiment  using  Pd-B 
and  the  earlier  investigations:  the  effects  of  some  of  these  differences  can  be  explained  in 
terms  of  the  earlier  results  while  some  of  the  results  are  surprising.  The  major  difference 
is  that  the  measurement  cycles  had  to  be  carried  out  at  rather  low  current  densities  [low 
for  the  observation  of  the  phenomenon]  in  view  of  the  relatively  high  surface  area  of  the 
electrode  [it  is  necessary  to  limit  the  power  input  to  the  cell  to  satisfy  the  design  criteria 
of  the  calorimeter].  As  the  rate  of  excess  enthalpy  generation  increases  markedly  with 
the  current  density  (l2),  the  values  achieved  in  the  experiment  were  necessarily  limited 
[the  rates  increased  to  ~  27  W  cm'3  on  Day  68  prior  to  evaporation  to  dryness].  A 
secondary  consequence  of  the  low  current  densities  was  that  the  electrode  was  polarized 
in  the  vicinity  of  the  region  for  the  onset  of  "positive  feedback"  for  most  of  the 
experiment  duration  [see  Fig.  A.2].  The  use  of  such  conditions  is  known  to  limit  the  rates 
of  excess  enthalpy  generation,  and,  in  the  limit,  may  destroy  the  phenomenon  (6)  (7)  [see 
Footnote  C.65]. 

The  major  unexpected  difference  for  this  Pd-B  electrode  has  been  the  observation  of  the 
development  of  "positive  feedback"  at  a  very  early  stage  of  the  experiments  [Day  3],  at  a 
low  current  density  and  at  a  low  temperature.  It  is  obviously  very  important  to  establish 
whether  this  early  establishment  of  "positive  feedback"  is  a  property  of  Pd-B  alloys  [the 
one  used  in  this  study]  or  whether  some  other  factor  is  involved.  The  marked  increase  in 
the  rate  of  excess  enthalpy  generation  at  temperatures  approaching  the  boiling  point  is 
probably  due  to  non-uniform  charging  [see  Footnote  C.66] 

A  major  feature  of  the  investigation  of  "Heat-after-Death"  in  the  experiment  is  the 
demonstration  that  the  rates  of  excess  enthalpy  generation  before  and  after  the  onset  of 
the  phenomenon  are  probably  identical.  Such  an  identity  would  be  expected  if  excess 
enthalpy  generation  takes  place  in  the  bulk  of  the  electrode,  but  these  effects  clearly 
require  further  investigation.  It  is  also  apparent  that  these  processes  relax  with  a 
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diffusional  relaxation  time  and  prolonged  maintenance  of  the  effects  evidently  requires 
special  conditions  (26)  [increase  of  the  temperature  of  the  heat  sinks].  Results  for  Day  25 
and  Day  26  shown  in  Figure  A-25  also  indicate  new  complications  in  the  establishment 
of  "Heat-after-Death"  and/or  "Positive  Feedback". 

This  investigation  of  this  Pd-B  experiment  has  demonstrated  yet  again  that  certain 
methods  of  data  evaluation  are  unsound  and/or  inaccurate  or  imprecise  [compare  e.g. 
(4),(5).(6).(7).(i3),(i4)j  Furthermore,  it  is  essential  to  avoid  the  effects  of  "positive  feedback" 
as  it  is  impossible  in  general  to  calibrate  closed  loop  systems  subject  to  such  feedback. 
Calibrations  can  only  be  achieved  if  the  effects  are  not  too  marked,  as  has  been  the  case 
for  Day  3  of  this  study.  Unfortunately,  it  is  almost  certain  that  the  investigations  carried 
out  by  NHE  have  relied  precisely  on  such  unsound  and  inaccurate  methods  of  calibration 
and  the  effects  of  "positive  feedback"  have  been  ignored.  However,  this  neglect  is 
probably  quite  general  and,  no  doubt,  accounts  for  many  of  the  contradictory  results  in 
this  field  of  research.  It  should  be  noted  that  much  of  the  pointless  controversy  in  this 
field  could  have  been  avoided  if  it  had  been  possible  to  replace  the  ICARUS  1  -  3 
Calorimeters,  Fig.  A.l,  by  the  ICARUS-4  version  [later  reclassified  as  the  ICARUS- 14 
Calorimeter]  Fig.  A. 27.  While  it  is  not  certain  that  this  particular  redesign  would  have 
eliminated  the  weak  time-dependence  of  the  heat  transfer  coefficients  observed  with  the 
ICARUS- 1  Calorimeter,  it  is  likely  that  this  would  have  been  true  and  that  these  systems 
could  have  been  developed  so  that  all  measurements  could  have  been  evaluated  with  a 
single,  predetermined  value  of  the  "true  heat  transfer  coefficient".  Finally,  it  is  important 
to  note  that  it  has  been  possible  to  achieve: 

(vi)  a  satisfactory  interpretation  of  evaporation  to  dryness  [Day  68] 

This  interpretation  has  had  to  take  into  account:  the  actual  barometric  pressure,  the 
change  of  the  boiling  point  of  the  solution  with  increasing  electrolyte  concentration, 
[saturation  of  the  electrolyte  -  not  discussed  in  the  present  Report],  changes  in  the  reflux 
ratio  [see  Footnote  C.67].  However,  prolonged  investigations  of  boiling  conditions  (l7)’(27) 
will  clearly  require  the  design  and  application  of  dual  calorimeters  such  as  the  ICARUS-9 
version  (l7)  (27).  it  is  also  important  to  determine  whether  the  marked  increase  of  the  rates 
of  excess  enthalpy  generation  at  temperatures  near  the  boiling  point  are  dependent  on  the 
establishment  of  boiling  conditions  or  are  simply  due  to  the  increase  in  temperature. 
While  it  is  certainly  desirable  to  develop  pressurized  systems  to  increase  the  boiling 
point,  significant  increases  in  the  boiling  point  could  also  be  achieved  by  using 
concentrated  electrolyte  solutions.  The  use  of  such  electrolytes  would  allow  the  extension 
of  the  range  of  applicability  of  the  ICARUS- 1  Calorimeters. 

Finally,  we  can  note  that  the  interpretation  of  this  experiment  gives  a  good  illustration  of 
the  need  to  evaluate  all  such  measurements  as  individual  "Case  Histories":  the  state  of 
development  of  research  in  this  field  in  1993  [when  the  first  ICARUS  System  was 
constructed]  was  certainly  not  at  the  point  at  which  such  interpretations  could  be  carried 
out  as  a  matter  of  routine.  Furthermore,  the  instrumentation  also  required  a  number  of 
additional  developments  to  facilitate  any  such  attempts  at  routine  evaluations.  The 
ICARUS- 14  System  [then  described  with  the  label  ICARUS-4]  was  to  be  the  next  step 
but,  as  has  already  been  noted,  this  modification  could  not  be  accomplished. 
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Section  B 


B.l  Introduction 

This  Section  contains  comments  on  the  modeling  of  the  Isoperibolic  Calorimeters  used  in 
the  ICARUS- 1  and  -2  Systems  leading  to  the  definition  of  a  number  of  versions  of  the 
heat  transfer  coefficient,  which  we  have  found  useful  in  our  investigations  [these  heat 
transfer  coefficients  define  the  behavior  of  the  calorimetric  systems].  The  content  of  this 
Section  is  based  on  Items  16  -  18  of  an  earlier  Document  (l3)  [see  Footnote  C.6]  and 
which  should  be  read  in  conjunction  with  the  present  Report.  The  earlier  Document (l3) 
also  contained  examples  of  data  evaluations  using  both  simulated  data  and  results  for 
"blank"  experiments  [with  Pt-cathodes;  compare  (l4)].  These  illustrations  <l3) demonstrated 
the  reasons  and  validity  of  the  choice  of  the  restricted  set  of  heat  transfer  coefficients 
specified  for  the  application  of  the  ICARUS- 1  Methodology 

B.2  Modeling  of  the  Calorimeters 

It  has  been  established  that  at  low  to  intermediate  cell  temperatures  [say  30°  <0<  80°C] 
the  behavior  of  the  calorimeters  is  modeled  adequately  by  the  differential  equation: 

CpM  (dA0/dt)  —  [Eccll(t)-Et|iermoneutral,bath]l  QKO 

change  in  the  enthalpy  input  rate  of  excess 

enthalpy  content  due  to  enthalpy 

of  the  calorimeter  electrolysis  generation 

+  AQH(t-t,)-  AQH(t-t2)  -  (3 1/4F [P/ { P * -P} ] [Cp, D2o.g  -Cp,d2o,.)A0+L] 

calibration  pulse  rate  of  enthalpy  removal  by  the  gas  stream  with 

Ethermoneutra!  referred  to  the  bath  temperature 

-  (kRo)03ba,h  [1-yt]  {f.(0)/03ba.h  +  4(pA0} 

time  dependent  effect  of  effect  of 

heat  transfer  radiation  conduction 

coefficient.  (B.l) 

With  the  calorimeters  supplied  with  the  ICARUS  Systems,  the  conductive  contribution  to 
heat  transfer  was  very  small.  In  fact,  if  this  term  was  "lumped"  into  the  radiative  term  by 
allowing  for  a  small  increase  in  the  radiative  heat  transfer  coefficient: 

Radiative  heat  transfer  =  (Icr'o)[1  -yt][(0bath  +A0)4-04bath  ]  (B.2) 

then,  the  values  of  the  "pseudo-radiative"  heat  transfer  coefficient  derived,  (kR'o)  )[  I  -yt], 
were  close  to  those  calculated  from  the  Stefan-Boltzmann  coefficient  and  the  radiative 
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surface  area  [typical  values  0.72  x  10'9  WK'4<(kR,o)[l-Yt]<0.76  x  10'9  WK"4].  However, 
for  the  Cell  No.38  used  in  this  experiment  this  "pseudo-radiative"  heat  transfer  coefficient 
is  found  to  be  ~  0.85  x  10'9  WK'4  so  that  the  conductive  contribution  was  evidently 
increased.  We  have  to  assume  that  this  increase  in  the  heat  transfer  coefficient  must  have 
been  due  to  a  "softening"  of  the  vacuum  in  the  Dewar  calorimeters  while  sitting  in  the 
laboratory  for  several  years. 

An  increase  in  the  "pseudo-radiative"  heat  transfer  coefficient  can  normally  only  be 
observed  if  the  cells  are  "overfilled"  with  D2O  during  the  periodic  replenishment  of  the 
cells.  This  "overfilling"  of  the  cells  leads  to  the  approach  of  the  electrolyte  level  to  the 
base  of  the  Kel-F  plug  sealing  the  cells  thereby  increasing  the  conductive  losses  through 
the  top  of  the  cell.  This  effect  leads  to  a  4  -  5%  increase  in  the  "pseudo-radiative"  heat 
transfer  coefficient  as  possibly  observed  in  the  results  for  Day  61  of  this  experiment  [see 
Figure  A. 20,  Section  A.5  and  Footnote  D.4], 

If  the  time  dependence  of  the  heat  transfer  coefficient  is  not  included  explicitly  in 
equation  (B.2)  then: 

Radiative  heat  transfer  =  (kR')  [(0bath+A0)4-04bath  ]  (B.3) 

where  the  radiative  heat  transfer  coefficient  (kR')  now  shows  a  weak  time-dependence. 

In  calculating  the  rate  of  enthalpy  removal  by  the  gas  stream 

(3I/4F)[P/{P*-P}][(Cp?o2o,g "  Cp,d2oj)  AO  +  L]  (B.4) 

we  have  always  assumed  that  the  partial  pressure  of  D20  [or  H20]  in  this  gas  stream  can 
be  calculated  using  the  Clausius-Clapeyron  equation  with  the  latent  heat  of  evaporation, 
L,  being  that  at  the  boiling  point.  Evaporative  cooling  only  becomes  an  important  term 
at  temperatures  close  to  the  boiling  point  [say,  at  A0>  70°C]  where  these  two 
assumptions  are  justified.  At  low  to  intermediate  temperatures,  AHevap(t)  is  a  minor 
correction  term  so  that  errors  due  to  the  two  assumptions  introduce  second  order  small 
quantities.  In  particular  the  errors  introduced  by  neglect  of  the  temperature  dependence 
of  the  latent  enthalpy  of  evaporation  are  <  0.1%  under  all  conditions  of  operation  of  the 
cells.  It  is  also  important  to  bear  in  mind  that  such  errors  are  further  reduced  for  all 
evaluations  of  the  "true"  heat  transfer  coefficients  as  these  evaluations  are  based  on 
differences  in  temperature  induced  by  the  calibration  pulses  [or  on  differences  in 
temperature  induced  by  "topping"  up  of  the  cells  or  perturbations  of  the  current  density: 
such  methods  of  calibration  are  not  considered  in  this  Report], 

However,  equation  (B.4)  does  not  give  the  rates  of  enthalpy  loss  tabulated  in  the  (kR')u 
spreadsheets  provided  by  the  NHE  analyses.  The  parameters  required  for  this  calculation 
were  contained  in  data  files  of  the  ICARUS- 1  and  ICARUS-2  software  and  were 
identical  for  both  systems.  The  values  installed  in  the  programs  as  supplied  were  [see 
Footnote  C.8]: 
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Cp.d2o,i  =  84.349  J  Mol1  K-' 

CpiD;0j  =  34.27  J  Mol'1  K'1 

0b»,i,»E  =  374.570  K 

F  =  96484.56  C  Mole'1 

R  =  8.314410  J  Mole'1  K'1 

L  =41,672.600  J  Mole'1 

P*  =  1.003  Ats 

The  Handbooks  (1),(2)  contained  specific  instructions  that  some  of  these  parameters  would 
need  to  be  changed  [here  Ogling  and  P*,  see  below]  as  well  as  instructions  as  to  how  such 
changes  in  the  parameter  listing  were  to  be  carried  out.  However,  it  appears  that  values 
of  the  rates  of  evaporative  enthalpy  loss  close  to  those  given  in  the  NHE  analyses  may  be 
calculated  for  low  to  intermediate  temperatures  using  the  parameter  listing  supplied  with 
the  instruments,  i.e.  the  changes  required  were  not  made.  [It  also  appears  that  the  latent 
enthalpy  of  evaporation  was  not  corrected  for  changes  in  temperature].  The  consequent 
errors  are  sufficiently  small  that  they  do  not  invalidate  the  analyses.  However,  the  values 
of  the  rates  of  evaporative  enthalpy  loss  contained  in  the  (kR1)  1 1  spreadsheets  of  the  NHE 
analyses  for  temperatures  close  to  the  boiling  point  cannot  be  calculated  with  any  values 
of  the  parameters  close  to  those  contained  in  the  listing  supplied  with  the  instruments. 
This  matter  has  not  been  investigated  further  as  it  is  in  any  event  necessary  to  make  three 
further  changes  if  experiment  cycles  close  to  the  boiling  point  are  to  be  evaluated  [see 
Footnote  C.9]. 

The  first  of  these  changes  is  the  adjustment  of  the  boiling  point  to  the  value  which  applies 
to  the  ambient  atmospheric  pressure.  It  should  be  noted  that  the  ICARUS-2  System  was 
supplied  with  the  means  for  the  continuous  recording  of  the  barometric  pressure,  but  this 
facility  was  evidently  disabled  following  the  installation  of  the  system  at  the  NHE 
laboratory.  The  values  of  P*  which,  have  been  used  in  the  present  interpretation,  have 
been  obtained  from  the  Sapporo  Airport  [thanks  to  Professor  Akito  Takahashi  of  Osaka 
University],  and  it  has  been  assumed  that  the  same  values  apply  to  the  experiments 
carried  out  in  the  NHE  Laboratory.  Furthermore,  it  has  been  assumed  that  the  pressure 
in  the  cell  is  the  same  as  the  ambient  pressure  although  it  may  well  be  that  the  pressure  in 
the  cell  was  somewhat  higher  than  this  value. 

The  second  change  is  that  it  is  necessary  to  take  note  of  the  fact  that  the  boiling  point 
corresponds  to  that  of  the  electrolyte  solution  in  the  cell.  It  has  been  assumed  that  this 
correction  is  given  by  that  for  an  ideal  solution 

^Pboiling  point  =  (R/L)(0  boiling  point)2  In  X,  (B.5) 

where  X(  is  the  mole  fraction  of  the  D20  in  the  electrolyte.  It  will  be  evident  that  this 
correction  becomes  especially  important  on  Day  68  when  the  cell  contents  are  driven  to 
dryness.  In  that  case  the  boiling  point  must  be  adjusted  at  each  measurement  interval  as 
the  D20  content  of  the  cell  decreases.  The  values  of  the  boiling  points  appropriate  for 
the  interpretation  of  the  experimental  data  for  Day  68  are  discussed  in  Section  A.7. 
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The  third  change  again  applies  specifically  to  Day  68,  namely,  an  allowance  for  the  effect 
of  reflux  in  the  cell.  In  order  to  evaluate  the  effects  of  reflux  we  need  to  take  note  of  the 
fact  that  the  vapor  space  in  the  cell  is  filled  predominantly  by  D2O  as  the  cell  is  driven  to 
dryness.  Thus,  even  at  the  start  of  Day  68  the  mole  fraction  of  D2O  in  the  vapor  space 
was  ~  0.85  for  this  experiment.  In  consequence,  heat  transfer  from  the  vapor  phase  to 
the  walls  of  the  Dewar  [to  provide  the  enthalpy  input  required  by  radiation  across  the 
vacuum  gap]  was  dominantly  from  the  D2O  content  of  the  vapor.  We  also  need  to  take 
note  of  the  fact  that  the  contribution  to  the  heat  capacitance  of  the  vapor  phase  in  the 
vicinity  of  the  boiling  point  due  to  the  D2O  content  of  this  phase  is 

-d(LP/P*)/dA0  =L2/R(0bolImg  poin,)2  exp[-LA0/R(0bojiing  Poim)2]  (B.6) 

where  A0  is  the  temperature  displacement  front  the  boiling  point.  This  heat  capacitance  is 
~  67  times  larger  than  that  of  an  equivalent  gas  space  filled  with  oxygen  and  hydrogen 
and  therefore  ~  380  times  larger  than  the  heat  capacitance  due  to  these  gases  for  the 
actual  working  conditions  at  the  start  of  day  68. 

This  marked  increase  in  the  heat  capacitance  of  the  part  of  the  cell  filled  with  gas  and 
D2O  vapor  has  two  consequences.  In  the  first  place,  the  heat  transfer  across  the  vacuum 
gap  must  be  maintained  at  the  same  value  as  that  which  applies  to  the  liquid  phase  [see 
Footnote  C.10].  Secondly,  the  radiative  output  across  the  section  of  the  Dewar  cells 
filled  with  vapor  must  be  balanced  by  the  condensation  of  an  amount  of  vapor  sufficient 
to  supply  the  radiative  enthalpy.  A  first  approximation  for  the  reflux  in  the  cell  is 
therefore. 

Rate  of  reflux  =  (kR’)|2  fi(0)  EAM/LAtM0  (B.7) 

where  AM  is  amount  of  D2O  evaporated  in  each  measurement  interval  At. 

Equation  (B.7)  will  be  an  upper  limit  for  this  extent  of  reflux  since  we  are  neglecting  the 
heat  transfer  to  the  walls  by  the  deuterium  and  oxygen  in  the  gas  space  as  well  as  the 
effects  of  the  re-heating  of  this  gas  space  by  the  liquid  in  the  lower  section  of  the 
calorimeter  [see  Footnotes  C.l  1,  D.5,  E.2], 

It  will  be  evident  that  analyses  based  on  the  use  of  equations  (B.4)  -  (B.7)  can  only  be 
approximations  for  a  number  of  reasons  in  addition  to  those  which  have  been  delineated. 
Two  of  the  most  obvious  deficiencies  are  the  use  of  "dilute  solution  theory"  in  the 
interpretation  and  the  neglect  of  hydrostatic  pressure  on  the  boiling  points  used  in  the 
Clausius-Clapeyron  equation.  It  follows,  therefore,  that  a  part  of  the  analyses  of  the  "raw 
data"  for  the  episodes  of  cells  being  driven  to  dryness  should  be  based  on  assumptions 
which  are  independent  of  the  use  of  equations  (B.4)  -  (B.7).  These  matters  are  considered 
further  in  Section  A. 7. 
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B.3  Definition  of  the  Heat  Transfer  Coefficients  Used  in  the  ICARUS 
Data  Analyses. 

The  heat  transfer  coefficients  will  be  described  by  the  suffices  used  previously  (13)  [as 
well  as  in  earlier  Reports  and  Correspondence]  as  defined  in  Section  A.2 


(kR')ij,i 


(B.8) 


where  i  =  1,2,3  denotes  "differential",  "backward  integration"  and  "forward  integration",] 
is  defined  at  appropriate  points  below  [see  also  Section  A.2]  and  1  =  1 ,2  denotes  lower 
bound"  and  "true"  respectively  [see  Footnote  C.12]. 


The  simplest  starting  point  is  to  assume 
system  i.e. 

Qr(t)  =  0 


that  there  is  no  excess  enthalpy  generation  in  the 

(B9) 


in  equation  (B.l)  and  to  evaluate  a  "lower  bound  heat  transfer  coefficient"  [i.e.  a 
coefficient  which  assumes  that  the  rate  of  excess  enthalpy  generation  is  zero]  at  a  time 
just  before  the  end  of  the  calibration  pulse,  t  =  t2'. 

(kR')l  =  [(Ecell(t)  -E,hermoneutral,ba.h)I-AHevap(t)  -  CpM(dA0/dt)]/ f,(0)  (B.  1 0) 

This  was  the  first  heat  transfer  coefficient  used  in  our  investigations,  hence  the 
designation  (kR')i .  It  should  be  noted  that  this  designation  should  really  be  changed  so  as 
to  be  consistent  with  the  definition  (B.8),  but  this  will  not  be  done  principally  because  the 
definition  (B.10)  was  subsequently  extended  to  any  part  of  the  measurement  cycle,  the 
coefficient  being  designated  (kR')ii.  We  should  perhaps  change  this  designation  of 

(kR')ioi  to  denote 

i  =  l,  differential 

j  =  0,  any  part  of  the  measurement  cycle 
1=1,  lower  bound 

but  the  description  (kR')n  will  be  retained  as  it  has  been  used  extensively  in  earlier  reports 
and  papers. 

Having  obtained  (kR')n,  it  is  frequently  desirable  to  establish  the  11-point  averages 
(]cr  ')i:i  so  as  to  decrease  the  "noise"  [see  Footnote  C.13,  D.6  and  E.j].  Such  averaging 

gives  ~  26  independent  values  for  measurement  cycles  lasting  1  day,  or  better  ~  52  values 
for  the  recommended  2-day  cycles.  In  turn,  it  isuseful  to  evaluate  the  6-point  averages 

of  (k/)n  which  have  been  desi§nated  as  (5V)u  •  11  is  not  useful  to  extend  this 

averaging  beyond  6  points  because  any  such  extension  makes  the  systematic  errors  [due 
to  the  residual  decrease  of  (kR')n  with  time]  larger  than  the  random  errors. 


The  values  of  (kR  ')u  and  (kR  ')n  [and  the  relevant  statistics]  were  to  have  been 
included  in  the  "(kR')n  -  spreadsheets"  but  it  is  apparent  that  the  group  at  NHE  did  not 
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carry  out  these  averaging  procedures. 


It  is  apparent  from  equation  (B.10)  that  we  need  accurate  values  of  CpM  to  make  (kR')n 
generally  useful,  but  it  is  also  apparent  that  the  group  at  NHE  retained  the  value  of  CPM 
specified  in  the  parameter  listing  rather  than  to  determine  the  correct  value  and  to 
substitute  this  corrected  value  in  the  listing  [see  Footnote  C.14].  A  first  approach  to  the 
determination  of  the  value  of  CpM  for  any  given  cell  is  to  rearrange  (B.10)  to  the  straight 
line  form 

y  =  mx  +  c  (B.ll) 

i.e. 

[{E  ccll(t)“E(hcrmoneutral,ba(h}I— AHCVap(t)]/  f|(0)  — 

[CpM(dA0/ dt)]/  f,(0)  +  (kR')ijj  (B.  12) 

and  to  derive  then  approximate  values  of  CpM  from  the  slopes  of  the  plots  in  regions 
where  the  temperature  is  varying  relatively  rapidly  with  time.  We  can  distinguish  four 
such  plots  designated  by  the  relevant  derived  heat  transfer  coefficients  (kR'o)i5|,  (kR'o))61, 
(kR'o),7i  and  (kR'o)m  according  to  whether  the  fitting  of  equation  (B.12)  is  carried  out  at 
times  somewhat  above  the  origin,  at  times  somewhat  above  fi  [the  time  of  application  of 
the  calibration  pulse],  at  times  somewhat  above  t2  [the  time  of  cessation  of  the  calibration 
pulse]  or  by  the  combination  of  the  last  two  time  regions,  see  Fig.  A.  1 1 .  However,  there 
is  a  measure  of  ambiguity  about  the  interpretation  of  the  values  of  (kRfo) j  j  |  derived 
which  has  been  discussed  in  Items  20  and  21  of (l3). 

It  should  be  noted  that  (kR’o)l5|  cannot  be  evaluated  systematically  for  this  Pd-B 
experiment  because  of  the  irregular  schedule  of  the  addition  of  D20,  see  Section  A.2. 
Evaluations  of  (kR'o)16|,  (kR'o)m  and.(kR'o)m  for  the  important  data  set  for  Day  3  are 
markedly  degraded  due  to  the  early  onset  of  "positive  feedback",  see  Section  A. 3.  The 
procedure  based  on  equation  (B.12)  has  limited  precision  because  of  the  need  to 
differentiate  the  inherently  "noisy"  experimental  data.  It  is  therefore  necessary  to  carry 
out  the  fitting  procedures  over  extended  regions  of  the  abscissa,  (dA0/dt)/ft(0),  so  that  the 
data  are  inevitably  affected  by  the  onset  of  the  "positive  feedback"  detected  for  the 
operation  of  the  cell  on  that  day. 


In  this  connection  it  should  also  be  noted  that  separate  investigations  have  shown  that 
(dA0/dt)  is  best  estimated  by  using  the  second  order  central  differences  [i.e.  the  chords  of 
the  curves].  More  accurate  values  could  be  derived  in  principle  by  using  higher  order 
differences.  However,  in  practice,  the  repeated  differentiation  of  the  experimental  data 
[implicit  when  using  higher  order  differences]  leads  to  an  increase  in  "noise"  if  we  use 
differences  higher  than  the  second  order  [see  also  Footnote  C.  1 5], 

In  the  absence  of  sufficiently  precise  determinations  of  CPM,  the  evaluations  must 
necessarily  be  restricted  to  regions  of  time  where  the  contribution  of  the  term  CPM 
(dA0/dt)/f|(0)  is  adequately  small.  In  that  case  it  is  adequate  to  use  a  "guesstimate"  of 
CPM.  This  matter  [including  the  evaluation  of  a  "guesstimate"  of  CPM]  is  considered 
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further  in  Section  A. 3. 


It  is  next  necessary  to  evaluate  a  "true  differential  heat  transfer  coefficient".  The  simplest 
procedure,  giving  (1cr')2  near  the  end  of  the  calibration  period  at  time  t  =  t2,  is  obtained 
by  including  the  calibration  pulse,  AQ: 


(kR')2—  {AQ  +  [Eccll  (A02,t2)  -  Ecc||(A0|,t2)]I  -  AHevap(A02,t2) 

+AHcvap(A0,,t2)  -  CpM[(dA0/dtU,2-  (dA0/dt)AOl.,2]}/f2(0)  (B.13) 

where  we  now  have 

f2(0)  ~  [Qballi  +  (A02,t2)]4  -  [9bath  +  (A0|,t2)]4  (B.14) 

[(kR')2  was  the  second  heat  transfer  coefficient  used  in  our  investigations].  In  order  to 
carry  out  such  evaluations,  it  is  useful  to  construct  A.4  or  A. 3  sized  plots  of  the  "raw 
data"  and  then  to  obtain  appropriate  averages  by  using  a  transparent  ruler.  This  type  of 
analysis  used  to  be  a  generally  accepted  approach  but  then  fell  into  disrepute.  However, 
the  methodology  is  now  again  accepted  giving  so-called  "robust  estimates". 

It  may  be  noted  that  the  errors  in  (kR')2  are  measures  of  the  accuracy  of  the  "true  heat 
transfer  coefficient"  as  the  estimates  are  made  in  terms  of  the  known  Joule  enthalpy  input 
to  the  calibration  heater.  Errors  in  (kR')i  or  (kR')u  are  measures  of  the  precision  of  the 
"lower  bound  heat  transfer  coefficients"  as  there  is  no  independent  calibration  and  there 
may  be  excess  enthalpy  generation  in  the  system.  It  is  important  that  (kR')i  and  (kR')2 
are  the  least  precise  and  least  accurate  coefficients,  which  can  be  obtained  from  the  "raw 
data".  Statements  that  the  errors  are  larger  than  this  [e.g.  see  (ll)]  simply  show  that 
mistakes  have  been  made  in  the  data  analysis  procedures  and/or  the  execution  of  the 
experiments. 

We  have  always  insisted  that  the  construction  and  evaluation  of  plots  of  the  "raw  data"  is 
an  essential  prerequisite  of  the  more  elaborate  data  evaluation  procedures.  For  one  thing 
it  shows  whether  the  "noise  levels"  in  the  experiments  were  sufficiently  low  to  justify 
more  detailed  evaluations  and  also  points  to  malfunctions  in  the  experiments.  It  also 
shows  immediately  whether  the  0-t  and  Eccii-t  transients  have  relaxed  sufficiently  to 
permit  the  evaluation  of  (kR')i  and  (kR')2.  Furthermore,  it  gives  immediate  indications  of 
the  presence  [or  absence]  of  "positive  feedback".  As  has  been  pointed  out  repeatedly,  all 
calibration  procedures  require  that  the  rate  of  excess  enthalpy  generation,  Q/(t),  be 
constant  during  the  calibration  periods.  These  matters  are  considered  further  in  the  main 
text.  Section  A. 3  and  A.4  [see  Footnote  E.5], 

Having  obtained  the  "true  heat  transfer  coefficient"  at  a  single  point  [usually  near  the  end 
of  the  calibration  pulse,  t  =  t2]  it  is  important  to  ask:  "what  is  the  true  heat  transfer 
coefficient,  (kR')i2,  at  any  other  time?"  We  can  make  such  an  evaluation  within  the 
duration  t|  <  t  <  t2  of  the  calibration  pulse  simply  by  using  equation  (B.13)  giving  (kR')i2 
rather  than  (kR')2.  Note  also  that  equation  (B.13)  can  be  rearranged  to  the  straight  line 
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form: 


{AQ  +[  Eccll(Ae2,t)  -  E°cell(A0i,t)]I  -  AHevap(Ae2,t)  +  AHeVap( A6 ,,t)}/f2(0) 

-  {CPM[  (dA0/dt)A92 ,  -  (dA0/dt)  A0, ,]}  /  f2(0)  +  (kR'o)162  (B.  1 5) 

which  is  applicable  at  times  close  to  and  above  fi.  It  is  evident,  therefore,  that  such  plots 
can  also  be  used  to  obtain  estimates  of  CPM,  but  the  accuracy  of  such  values  is  inevitably 
much  lower  than  the  precision  of  those  obtained  by  the  application  of  the  corresponding 
expression  for  the  "lower  bound  heat  transfer  coefficient,  (kR'o)i6i"  equation  (B.12). 
Nevertheless,  equation  (B.15)  is  useful  because  it  allows  the  removal  of  the  effects  of  the 
water  equivalent,  CPM,  on  the  "true  heat  transfer  coefficient,  (kR'o)i62",  simply  by 
extrapolating  to  zero  value  of  the  abscissa.  However,  the  time  corresponding  to  this  point 
will  not  be  accessible  experimentally  for  calibrations  carried  out  with  a  calibration  pulse 
of  6-hour  duration  for  polarizations  carried  out  at  low  cell  currents  [although  this  time  is 
probably  close  to  t  =  t2,  see  Footnote  C.  1 6]. 

In  the  regions  in  which  there  is  no  application  of  a  heater  pulse  i.e.  for  0  <  t  <  ti  and  t2  <  t 
<  T,  the  "true  heat  transfer  coefficient"  can  only  be  obtained  from  the  "heating"  and 
"cooling  curves"  i.e.  the  "driving  force"  is  the  change  in  the  enthalpy  content  of  the 
calorimeters  rather  that  AQ.  It  is  now  sensible  to  cast  equation  (B.  1 3)  in  the  form: 

{CpM[(dA0/dt)A02 1  -  (dA0/dt)A0(  t]}  /  f2(0)  =  -  (kR'o)i52  or  (kR'o)i72  + 


{ [Ecei,(A02,t)-Ecc„(A0 ,  ,t)]I-AHevap(A02,t)+  AHevap(A0 ,  ,t)  }/f2(0)  (B.  1 6) 

If  the  system  is  functioning  correctly,  then  it  will  be  found  that  the  L.H.S.  of  equation 
(B.16)  is  essentially  constant  [although  this  constancy  can  only  be  probed  over  a  short 
time  range].  The  second  term  on  the  R.H.S.  of  equation  (B.16)  will  be  much  smaller 
than  the  term  on  the  L.H.S.,  i.e.  it  is  in  the  nature  of  a  correction  term  to  give  "point-by¬ 
point"  values  of  (kR'o)i52  or  (kR'o)i72.  It  will  be  evident  that  the  accuracy  of  these 
versions  of  the  "true  heat  transfer  coefficient"  is  limited  by  the  accuracy  of  the  estimates 
of  CpM.  This  particular  part  of  the  methodology  is  therefore  only  useful  to  serve  as  a 
check  on  the  operation  of  the  cells  and  methods  of  data  evaluation.  Furthermore,  it  is  not 
possible  to  apply  equation  (B.16)  systematically  to  the  time  region  0  <  t  <  t,  for  this 
experiment  in  view  of  the  irregular  schedule  of  addition  of  D20  to  the  cell. 

As  has  been  noted  previously  (l3),  we  have  been  unable  to  combine  data  in  the  regions 
just  above  fi  and  t2  to  give  a  simple  equation  leading  to  (kR'o)ig2. 

The  assumption  underlying  this  part  of  the  account  presented  in  this  Report  is  that  we 
can  only  determine  (kR')i2  within  the  duration  of  the  calibration  pulse  t|  <  t  <  t2,  and, 
at  a  lower  accuracy,  (kR'ojia  and  (kifojra  in  regions  adjacent  to  the  origin  and  for  times 
adjacent  and  above  t2  respectively.  However,  this  conclusion  is  incorrect.  We  need  to 
make  the  additional  assumption  that  the  rate  of  any  excess  enthalpy  generation  is  constant 
during  any  particular  calibration  period  in  order  to  determine  (kR')|2.  This  means  that  we 
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can  only  obtain  a  single  value  of  this  heat  transfer  coefficient  per  calibration  period  and, 
consequently,  a  single  value  of  (kR')|2  -  (kR')n.  Two  important  points  follow  from  this 
conclusion.  In  the  first  place,  the  precision  of  (kR')i2  must  be  very  nearly  equal  to  the 
precision  of  (kR')i  i-  Secondly,  and  related  to  the  first  point,  we  see  that  if  we  extend  the 
assumption  that  the  rate  of  excess  enthalpy  production  is  constant  during  the  period  t|  <  t 
<  t2  to  saying  that  it  is  constant  for  the  whole  measurement  cycle,  0  <  t  <  T,  then  it  is 
immediately  possible  to  derive  (Icr’) i 2  over  the  whole  of  this  cycle.  Thus,  if  the 
difference  between  the  "true"  and  "lower  bound"  heat  transfer  coefficients  can  be 
established  at  any  one  time  [say  A(kR'),  at  time  t2]  then  [kR'(t)]i2  at  any  other  time  will  be 
given  by 

[(kR'(t)]i2  =  [(kR'(t)]n  +  A(kR*)t2f|(0)t2  /  fi(0)i  (B-17) 

The  ratio  f(0),2  /  fi(0)t  is  of  order  unity  which  implies  that  the  shift  (kR')i 2  -  (kR')n  is 

always  close  to  that  at  the  calibration  point.  Equation  (B.17)  also  points  to  a  further 
important  conclusion.  At  the  time  of  writing  of  the  Handbook  for  the  ICARUS-1 
Systems  (2)we  believed  that  the  precision  of  (kR')12  [and  of  other  versions  of  the  "true 
heat  transfer  coefficient"]  would  always  be  given  by  the  accuracy  of  that  coefficient 
which  is  certainly  lower  than  the  precision  of  (kR')n  [see  Footnote  C.17].  Equation 
(B.17)  shows  that  is  incorrect:  the  precision  of  (kR')i2  is  very  nearly  equal  to  the 
precision  of  (kR’)i  1  -  It  follows  that  changes  in  the  rates  of  excess  enthalpy  production 
can  be  established  at  the  same  level  of  precision  as  that  of  (kR')n.  The  same  comments 
apply  to  the  precision  of  the  "true  heat  transfer  coefficient,  (kR')22"  relative  to  that  of  the 
"lower  bound  heat  transfer  coefficient,  (kR’)2i”  which  is  discussed  below.  In  consequence, 
the  changes  in  the  rates  of  excess  enthalpy  production  can  be  established  with  relative 
errors  <  0.01%,  and  these  errors  determine  the  level  of  significance  with  which  such 
changes  can  be  discussed.  Of  course,  the  accuracy  of  the  true  heat  transfer  coefficients 
remains  determined  by  the  errors  of  differences  such  as  that  of  (kR')i2  -  (kR’)n- 

It  is  important  here  to  stress  once  again  that  any  attempt  to  calculate  the  variation  of  rates 
of  excess  enthalpy  generation  within  the  measurement  cycles  must  also  pay  due  regard  to 
the  fact  that  it  is  not  possible  to  calibrate  the  systems  if  the  rate  of  excess  enthalpy 
generation  varies  with  time.  It  is  also  important  that  this  comment  applies  equally  to  any 
calorimetric  system  which  we  might  wish  to  use.  If  the  rate  of  excess  enthalpy  generation 
does,  in  fact,  vary  with  time,  then  A(kR’)  must  be  derived  from  separate  experiments.  This 
is  the  situation,  which  applies  to  the  experiment  as  is  discussed  in  Section  A. 3  and  A. 4. 
The  comments  made  in  this  part  of  Section  B.3  should  be  read  in  conjunction  with  ^  \ 

The  discussion  of  the  accuracy  of  "true  heat  transfer  coefficients"  versus  the  precision  of 
the  "lower  bound  heat  transfer  coefficients"  prompted  our  search  tor  methods,  which 
would  increase  both  the  precision  and  accuracy.  The  reason  for  the  limited  precision  of 
(kR*)i  1  and  accuracy  of  (kR*)i2  is  mainly  due  to  the  need  to  differentiate  "noisy" 
experimental  data  sets  in  order  to  derive  CpM(dA0/dt)  [see  Footnotes  C.18  and  D.8]. 

If  we  wish  to  avoid  the  numerical  differentiation  of  the  experimental  data  sets,  then  we 
can  rely  instead  on  the  numerical  integrations  of  these  data  and  compare  these  to  the 
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integrals  of  the  differential  equation  representing  the  model  of  the  calorimeters.  For  the 
backward  integrals  starting  from  the  end  of  the  measurement  cycles  at  t  =  T  we  obtain 


(kR)2i  — 


Jnet  enthalpy  input(x)dx  CpM[A0(t)-  A0(T)] 


/  f|(0)dx 

T  T 


Jf,(0)dT 


(B-18) 


while  the  corresponding  equation  for  forward  integration  from  the  start  of  the 
measurement  cycle  is 


_  )netenthalpyinput(T)dT  CpM[A8(t)- A0(O)] 


|f,  (0)dx 


J  f.  (0)dx 


(B.19) 


Here,  the  suffices  21  and  31  denote  respectively  "backward  integration,  lower  bound"  and 
"forward  integration,  lower  bound".  (kR')2 1  and  (kR')3i  are  the  corresponding  integral 
heat  transfer  coefficients  defined  at  time  t. 


We  have  to  take  note  of  the  fact  that  care  is  needed  when  integrating  the  terms  fi(0)  and 
[net  enthalpy  input  (x)]  around  the  discontinuities  at  t  =  ti  and  t  =  t2  [also  the  times  t  =  0 
and  t  =  T  if  the  range  of  the  integrations  is  extended].  The  trapezium  rule,  Simpson's  rule 
or  the  mid-point  rule  has  been  used  at  various  times  to  carry  out  the  integrations.  Of  these 
rules,  only  the  mid-point  rule  is  strictly  speaking  correct  in  that  it  agrees  with  the 
mathematical  definition  of  an  integral.  It  is  quite  generally  assumed  that  integrations 
carried  out  using  the  trapezium  or  Simpson's  rule  will  converge  onto  the  "correct" 
algebraic  result  if  the  integration  interval  is  made  adequately  small,  but  this  does  not 
necessarily  follow.  This  is  a  matter  which  needs  to  be  investigated  for  each  particular 
case. 


It  may  be  noted  that  the  only  straightforward  way  in  which  we  can  integrate  around  the 
discontinuities  at  t  =  ti  and  t  =  t2  is  by  means  of  the  trapezium  rule,  and  this  is  the  method 
which  has  been  used  in  the  recalculations  presented  in  the  this  Report.  If  the  times  of 
application  and  cessation  of  the  heater  calibration  pulses  correspond  exactly  to  and  t2, 
respectively,  then  we  can  carry  out  the  integrations  around  the  discontinuities  by  inserting 
extra  data  points  at  these  times.  It  appears  that  the  data  sets  in  this  experiment  satisfy  this 
criterion,  although  this  is  not  generally  true  for  all  experiments  carried  out  with  the 
ICARUS-2  system;  lack  of  synchronization  of  the  calibration  pulses  with  t|  and  t2 
appeared  to  be  generally  true  for  measurements  with  the  ICARUS- 1  systems.  In  that 
case,  it  is  necessary  to  determine  these  times  separately  [this  can  be  done  adequately  from 
the  0  -  t  plots]  so  as  to  establish  the  integration  intervals,  and  it  is  then  necessary  to  insert 
4  additional  data  points. 

The  evaluations  of  (kR  )2|  and  (k  )22  [see  below]  and  of  (kR')3i  and  (kR')32  [see  also  below] 
were  to  have  been  carried  out  using  (kR')2|  and  (kR')3i  spreadsheets  produced  by  the 
software.  As  we  have  never  had  access  to  these  spreadsheets  [if,  in  fact,  they  were  ever 
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produced],  we  cannot  now  establish  whether  the  integrations  around  the  discontinuities 
were  carried  out  correctly,  although  we  believe  that  they  must  have  been  in  error.  In  any 
event,  all  the  integrations  used  in  the  evaluations  described  in  this  Report  have  been 
carried  out  using  the  "raw  data". 

The  adequacy  [or  inadequacy]  of  the  particular  integration  procedures  coupled  to  the 
adequacy  of  the  chosen  integration  interval  is  revealed  more  clearly  when  we  come  to  the 
use  of  equations  (B.18)  and  (B.19)  to  determine  CPM  and  to  carry  out  extrapolations  to 
remove  the  effects  of  the  second  term  on  the  R.H.S.  of  these  equations  on  the 
corresponding  heat  transfer  coefficients.  The  procedure  set  out  in  the  Handbook  tor  the 
ICARUS- 1  System  (2)  restricted  the  integrations  to  the  time  region  of  the  application  of 
the  heater  calibration  pulse.  For  backward  integration  we  obtain  (B.20) 

j  net  enthalpy  inpul(T)dt  c  M[Ae(t)  _A0(u)] 

U - T - =  - L+(k’*o)26l  (a20) 

jf1(0)dx  lf1(0)dx 

t  o  t  2 

while  for  forward  integration,  we  obtain 

j  net  enthalpy  input(T)dx  £  M[A0(t)_Ae(t|)] 

L - j - — E-f - +(kRO)36,  (B.21) 

Jf1(0)dx  |  fj(9)dx 

t,  C 

It  has  been  found  [e.g.  see  references  3,  4,  5  and  13]  that  equation  (B.20)  can  be  used  to 
derive  accurate  values  of  CpM  while  there  is  some  minor  degradation  when  using  forward 
integration,  equation  (B.21).  The  application  of  equation  (B.20)  to  the  data  sets  was  the 
"target  methodology"  of  the  ICARUS  Systems  [e.g.  see(2)]  and  the  derived  "lower  bound 
heat  transfer  coefficient,  (kR'o)26i"  was  described  as  (kR')2i  in  the  Handbook  and  the 
associated  correspondence.  We  have  since  then  used  the  more  extended  description, 
(kR'o)26i,  to  denote  the  fact  that  with  j  =  6,  we  are  carrying  out  the  evaluation  in  the  time 
region  t|  <  t  <  t2-  The  same  types  of  evaluation  may  be  used  to  derive  ^’0)251,  ^<’0)271 
and  (kR'o)28i  as  well  as  (kR’o)35i,  (kR'o)37i  and  (kR’o)38|.  It  is  only  necessary  to  start  the 
integrations  from  the  appropriate  times,  which  also  give  the  starting  values  of  0  for  the 
R.H.S.  of  the  relevant  equations.  Of  these  sets  of  estimates,  that  leading  to  (kR'o)28i  is 
especially  useful  and  this  particular  fit  also  gives  good  estimates  of  CPM.  However,  it 
should  be  noted  that  it  is  necessary  to  use  care  in  applying  these  procedures  to  the  data 
for  Day  3  of  this  experiment  because  of  the  early  onset  of  "positive  feedback"  see  Section 
A. 3. 


In  order  to  obtain  the  "true  heat  transfer  coefficients"  it  is  necessary  to  combine  the 


integrals  of  the  enthalpy  inputs  in  equations  (B.20)  and  (B.21)  with  thermal  balances 
made  at  one  or  a  series  of  points.  This  can  be  done  in  a  number  of  ways,  and  it  is 
important  that  this  part  of  the  evaluation  (2^  was  changed  in  the  summer  of  1994  following 
the  receipt  of  the  first  two  sets  of  data  collected  by  NHE (3),(4).  Attention  will  be  confined 
here  to  the  procedure  originally  suggested  in  the  Handbook  for  the  ICARUS- 1  System  (2) 
[see  Footnote  C.19].  We  make  a  thermal  balance  just  before  the  application  of  the 
calibration  pulse  and,  if  the  system  has  relaxed  adequately  and  d0/dt  =  0,  then  we 
consider  (kR')32  • 

0  =  [Net  enthalpy  input  (ti)][t-  t|]  +  Qf  [t-  ti]  - 

(kR')32{[  (6ba,h  +  A0  (t,)]4  -  04bath}  [t-t,]  (B.22) 

The  combination  of  this  equation  with  equation  (B.19)  eliminates  the  unknown  rate  of 
excess  enthalpy  generation.  We  obtain 


(kk)32  - 


j  net  enthalpy  input(x)dx-  net  enthalpy  input(t  j  )j[t  -  t,] 


t 

Jf2(0)dx 

t. 


CpM  A0(t)  -  A0(tj ) 


Jf2(0)dx 

ti 

(B.23) 


The  corresponding  equation  for  (kR’)22  follows  from  equation  (B.23)  on  replacing  ti  by 
t2  [see  equation  (B.24)  below]. 


We  note  again  that  the  group  at  NHE  did  not  follow  the  instruction  in  the  ICARUS- 1 
Handbook  (2)  to  use  measurement  cycles  of  2-day  duration  and,  for  the  reduced  time 
scales  of  1-day  cycles  in  particular,  it  is  necessary  to  include  the  term  CPM  (dA0/dt)  in  the 
thermal  balances,  equation  (B.22).  However,  the  team  at  NHE  continued  to  use  the 
original  form  of  the  equation  [see  Footnote  C.20],  It  also  appears  that  NHE  did  not 
follow  the  instruction  (2)  to  evaluate  (kR')32  at  times  close  to  t2.  This  matter  is  discussed 
further  in  Section  A.4. 


It  is  convenient  also  to  rewrite  the  derived  equation  for  (kR')22in  the  "straight  line  form". 

t 

J  net  enthalpy  input(x)dx  —  [net  enthalpy  input(t2)][t  —  t2] 

1 2 


j  f2(0)dx 
*2 

CpM[A0(t)-  A0(t2)] 

t  h  (  kR  °^262 

J  f9(0)dx 

t->  ^ 


(B.24) 
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(kR')22  and  (kR'o)262  were  the  versions  of  the  "true  heat  transfer  coefficient"  which  we  used 
in  our  investigations  prior  to  the  construction  of  the  ICARUS-1  System.  As  we  did  not 
wish  to  discuss  the  differences  between  (kR')32,  (kR'o)362,  (kR')22  and  (kR'o)262  and,  as  we 
expected  (kR’)32  to  converge  onto  (kR’)22  for  the  specified  2-day  measurement  cycles 
[within  the  error  limits  specified  for  the  ICARUS-1  System]  we  also  labeled  (kR')32as 
(kR')22  [see  Footnote  C.21]. 

It  should  be  noted  that  the  extrapolation  (B.24)  automatically  removes  the  effect  of  the 

t 

term  CpM[0(t)  —  0(^-2)]/  j*  f2(0)dx  on  the  true  heut  transfer  coeflicient  .  This  npplicntion 

t  t 

of  equation  (B.24)  [and  of  (kR'>22  evaluated  close  to  the  mid-point  t  =  t2]  was  one  of  the 
major  objectives  for  our  methodology  because  CPM  is  the  least  accurate  paiameter  in  the 
analysis. 

While  it  is  also  possible  to  write  equation  (B.23)  in  the  form  (B.24)  to  give  (kR'o)362,  this 
method  of  analysis  is  not  useful  as  the  range  of  the  extrapolation  required  is  too  long  ( 
[see  also  (l3)].  For  this  reason  it  was  recommended  in  the  ICARUS-1  Handbook  that 
(kR’)32  should  be  evaluated  at  times  close  to  t  =  t2  using  values  of  CPM  determined  from 
applications  of  equation  (B.21).  However,  in  view  of  the  errors  in  the  determination  of 
CPM,  these  values  of  (kR')32  are  inevitably  less  accurate  than  those  of  (kR')22  or  (kR'o)262 
[see  also  Section  A.4], 

We  note  also  that  care  must  be  taken  in  carrying  out  the  required  linear  regression  fitting 
procedures  as  is  illustrated  in  Section  A.4  [see  references  4,  5  and  13]. 

We  should  observe,  furthermore,  that  equation  (B.24)  is  soundly  based  [in  a 
mathematical  sense]  in  that  the  extrapolation  to  [A0(t)  -  A0  (t2)]  =  0  gives  the  value  of 
(kR'o)262  at  a  well  defined  time,  t  =  t2.  This  is  equally  true  of  all  of  the  coefficients  based 
on  "forward"  or  "backward"  integration;  however,  the  starting  points  for  these 
integrations  will  usually  be  chosen  to  be  t  =  0,  t  =  t|  or  t  =  T,  and  the  definition  of  the  heat 
transfer  coefficients  at  these  points  is  not  generally  useful.  The  exception  here  is  the 
"lower  bound  heat  transfer  coefficient,  (kR'o)26i"  which  is  also  defined  at  time  t  —  t2. 
We  observe  also  that  (kR’o)26i  and  (kR'o)262  are  the  most  precise  and  accurate  values  of 
the  "lower  bound"  and  "true  heat  transfer  coefficients"  which  can  be  derived  with  the 
methodology  as  presently  developed.  Furthermore, 

(kR'o)26i  =  (kR')n  at  t  =  t2  and  (kR'o)262=  (kR')i2  at  t  =  t2  (B.25) 

so  that  the  best  value  of  A(kR')t  that  can  be  obtained  for  use  in  equation  (B.  1 7)  is 

A(kR')t  =  (kR')i2  -  (kR')n  =  (kR'o)262  -  (kR'o)26i  (B.25A) 

This  "sound  basis"  of  the  heat  transfer  coefficients  derived  by  "forward"  and  "backward 
integration"  should  be  contrasted  with  the  corresponding  position  for  the  "differential 
heat  transfer  coefficients"  which  has  been  discussed  above. 
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In  the  final  part  of  this  Section,  we  need  to  consider  somewhat  further  the  time- 
dependence  of  the  various  forms  of  the  heat  transfer  coefficient  [compare 
We  observe  first  of  all  that  we  are  interpreting  here  the  systematic  variations  of  typically 
0.4%  of  the  differential  or  0.2%  of  the  integral  coefficients.  The  only  reason  why  we  are 
able  to  investigate  systematic  variations  of  such  small  quantities  is  the  very  high 
precision  of  the  methods  of  data  evaluation. 

We  observe  secondly,  that  as  the  differential  coefficients  are  evaluated  at  local  times, 
they  will  show  the  weak  time  dependence: 

(kR’)  =  (k’Ro)[l  -yt]  (B.26) 

[c.f.  equations  (B.2)  and  (B.3)].  In  the  definition  of  the  integral  heat  transfer  coefficients 
given  in  this  Section,  (kR)  has  been  regarded  as  being  constant  whereas  the  investigation 
of  the  differential  heat  transfer  coefficients  shows  that  we  should  really  include  the  time- 
dependence,  equation  (B.26),  i.e.  we  must  use  equation  (B.2)  in  the  integrations. 
Integration  of  this  equation  gives 

(kRo)  [/f](0)dx- 7tjf1(0>dx+  y}|f1(0)d'cdx]  (B.27) 

If  we  now  regard  fi(8)  as  being  constant  throughout  the  measurement  cycle  [which  is  a 
rough  approximation  for  the  case  of  the  "lower  bound  heat  transfer  coefficients"]  then  the 
integral  becomes 

(kRo)f,(0)t[l  —  yt/2]  (B.28) 

It  follows  that  the  heat  transfer  coefficients  given  by  equations  (B.18)  and  (B.19)  are 
given  by 

(kR)2i  =  (kRo)2I[l  + y(T-t)/2]  (B.29) 

and 

(kR)3I  =  (kRo)31[l  -yt/2]  (B.30) 


within  the  limits  of  the  approximation.  (kR’o)2i  and  (kR'o)3i  are  respectively  the  values  of 
(kR')2i  at  t  =  T  and  of  (kR')3|  at  t  =  0.  It  follows  that  the  slopes  of  the  plots  of  (kR')2i  and 
(kR')3i  versus  time  are  roughly  one  half  of  the  plot  of  (kR')n  versus  time. 

Equation  (B.27)  also  shows  the  way  in  which  we  can  test  whether  the  characteristics  of 
the  Dewar  cells  can  be  described  by  a  single,  time-independent  heat  transfer  coefficient. 
Thus  evaluation  of  (kR  )2i  according  to  equation  (B.18)  gives  us  the  heat  transfer 
coefficient 

(kR)2)  =  (k  Ko)2l  [1  -yt  +  yj  J  fl(0)dTdx/ifl(e)dx]  (B.3 1) 

I  T  T 
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so  that  the  "time-independent  heat  transfer  coefficient,  (kR'o)2i"  is  readily  determined. 
The  fact  that  heat  transfer  from  the  cells  can  be  represented  by  such  a  single  time- 
independent  heat  transfer  coefficient  has  been  demonstrated  several  times  [e.g.  see  Fig. 
51  of  reference  13].  Indeed,  such  representations  are  the  basis  of  our  statement  that  the 
"integral  lower  bound  heat  transfer  coefficients"  can  be  determined  with  a  precision  given 
by  relative  errors  of  less  than  0.01%  [see  Footnote  C.22], 

The  variations  of  (kR')iu  (kR')2i  and  (kR')3i  with  time  show  that  this  time  dependence  of 
the  heat  transfer  coefficients  must  be  taken  into  account  in  evaluations  of  the  rates  of 
excess  enthalpy  generation  aiming  at  the  highest  achievable  accuracy.  If  this  is  not  done, 
then  the  values  of  the  heat  transfer  coefficients  at  the  mid-points  t  =  t2,  should  be  used. 
In  that  case  the  values  of  the  rates  of  excess  enthalpy  generation  calculated  will  be 
slightly  too  small  for  t  <  t2  and  slightly  too  large  for  t  >  t2.  However,  the  total  excess 
enthalpy  calculated  for  a  complete  measurement  cycle  will  be  approximately  correct  [see 
Footnote  C.23].  We  must  also  note  that  the  differential  heat  transfer  coefficient,  (kR’)i2, 
must  be  used  in  the  evaluations  of  the  rates  of  excess  enthalpy  generation  and  the  integral 
heat  transfer  coefficients  in  the  evaluation  of  the  excess  enthalpy  [including  the  total 
excess  enthalpy  for  complete  measurement  cycles].  In  particular,  the  use  of  (kR')22  in  the 
evaluation  of  the  rates  of  excess  enthalpy  generation,  will  underestimate  these  quantities 
[see  Footnote  C.24], 

B.4  The  Specification  of  the  ICARUS-1  Data  Evaluation  Procedures  and 
Experimental  Protocols 

The  modeling  of  the  ICARUS- 1  Type  Calorimeters,  Fig.  A.l  has  been  investigated 
repeatedly  by  means  of  the  evaluation  of  data  sets  for  appropriate  "blank  experiments" 
[using  in  the  main  Pt-cathodes  polarized  in  D20-based  electrolytes].  The  objective  here 
has  been  the  definition  of  the  appropriate  "instrument  function".  It  has  been  found  that 
this  can  be  accurately  defined  by  equation  (B.l). 

The  next  step  in  this  initial  phase  of  the  work  has  been  the  definition  of  a  set  of  heat 
transfer  coefficients,  which  characterize  the  behavior  of  the  calorimeters  and  the 
investigation  of  their  precision  and  accuracy  leading  up  to  their  use  in  evaluating  the  "raw 
data  sets"  of  the  experimental  measurement  cycles.  The  "raw  data"  used  in  these 
investigations  have  been  both  those  for  the  appropriate  "blank  experiments"  and  those 
generated  by  simulations  of  the  cell  behavior.  An  illustration  of  this  phase  of  the 
investigation  has  been  given  (13)  [see  also  references  4.  5,  9,  10.  and  14], 

The  outcome  of  these  investigations  has  been  the  demonstration  that  it  is  useful  to 
determine  first  of  all  the  time-dependence  of  the  "differential  lower  bound  heat  transfer 
coefficient,  (kR')n",  as  well  as  of  the  derived  means,  (kR  ')51  and  (kR  ')„ .  However, 
these  coefficients  have  a  limited  precision  because  their  evaluation  requires  the 
differentiation  of  the  inherently  "noisy"  experimental  data.  Precise  and  accurate 
evaluations  are  therefore  best  based  on  the  "integral  lower  bound  heat  transfer  coefficient, 
(kR')2i"  and  the  "integral  true  heat  transfer  coefficient.  (kR')22"  as  well  as  on  the  values 
(kR'o)25i  and  (kR'o)2s2  derived  in  the  extrapolation  procedures.  These  extrapolation 
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procedures  lead  both  to  the  elimination  of  the  effects  of  the  water  equivalent,  CPM,  on 
their  values  as  well  as  to  reasonably  accurate  determinations  of  CPM.  The  differences 
between  the  "true"  and  "lower  bound  heat  transfer  coefficients",  (ki<’)22  -  (ki<”)n  or 
(kR'o)252  -  (ki^ofs i  can  then  be  used  to  define  the  "differential  true  heat  transfer 
coefficient,  (kR')i2". 

It  has  been  found  that  the  precision  and  accuracy  of  the  integral  heat  transfer  coefficients 
is  so  high,  that  it  is  possible  to  investigate  their  systematic  variations  with  time  [typically 
the  systematic  variations  of  just  ~  0.2%  of  their  numerical  values].  Furthermore,  it  is 
possible  to  reduce  such  data  to  a  single,  time  independent  heat  transfer  coefficient,  e.g.  of 
(kR'o)2i  with  relative  errors  below  0.01%.  This  result  is  hardly  surprising:  the  "physics" 
of  the  calorimeters  are  quite  simple  [they  are  "ideal  well-stirred  tanks"]  and  the  errors  are 
mainly  those  set  by  the  temperature  measurements.  It  is  also  relatively  straightforward  to 
specify  the  changes,  which  would  need  to  be  made  to  reduce  the  errors  further  -  say  to 
0.001%  -  if  that  should  ever  prove  to  be  necessary  or  desirable. 

Although  the  precision  and  accuracy  of  the  heat  transfer  coefficients  based  on  the 
"forward"  integration  procedures,  (kR'Xn  and  (ki<’)32,  was  known  to  be  lower  than  those 
based  on  the  "backward"  integrations,  (kR')2i  and  (kR')22,  the  ICARUS-1  Methodology 
was  nevertheless  based  on  such  "forward"  integrations  l2)  because  such  "forward" 
integration  was  easier  to  implement  and  to  combine  with  the  evaluations  of  the  data  sets. 
It  was  anticipated  that  the  extension  of  the  measurement  cycles  from  1  to  2  days  and,  in 
particular  of  the  calibration  periods  from  6  to  12  hours,  would  allow  the  determination  of 
(kR')3i  and  (kR'^with  the  required  and  specified  precisions  and  accuracies  (2).  These 
changes  in  the  measurement  cycles  were  also  expected  to  facilitate  other  parts  of  the 
investigation  such  as  the  determination  of  the  "true  heat  transfer  coefficient  (kRV'-  The 
production  of  plots  of  the  "raw  data"  and  the  inspection  of  these  plots  leading  to  the 
graphical  evaluation  of  (kR’)i  and  (kR')2  were  to  be  the  first  step  in  the  data  processing. 

Unfortunately,  the  protocols  laid  down  in  the  Handbook  for  the  ICARUS- 1  System  <2) 
were  not  followed  in  the  experiments  carried  out  by  the  team  at  NHE  [see  Footnote 
C.25],  Furthermore,  following  the  receipt  of  the  first  set  of  data  for  experiments  carried 
out  in  the  Sapporo  Laboratories,  it  became  clear  that  there  were  timing  errors  in  the 
ICARUS- 1  System  (:>).  These  timing  errors  did  not  affect  the  determination  of  (k|<’)2 1  and 
(kR')22  [see  Footnote  C.26].  It  was  therefore  recommended  (,)  that  the  protocol  set  down 
in  the  Handbook  (2>  be  strictly  adhered  to,  that  the  preliminary  evaluations  should  be 
based  on  (ku')i.  (kR1)’  and  (kR’)n  and  that  the  final  evaluation  should  be  based  on  (kn')2i. 
(kR')22,  (kR 0)251  and  (kR'o)2S2-  It  is  evident  that  these  instructions  were  ignored. 

The  development  of  the  various  aspects  of  the  data  analysis  described  in  the  sub-Section 
B.4  is  illustrated  in  part  by  the  analysis  of  the  experiment  described  in  Section  A. 4. 
Inevitably,  this  illustration  is  incomplete  because  of  the  very  early  development  of 
"positive  feedback"  in  this  experiment. 
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Section  C 


As  is  explained  in  the  main  text.  Section  A,  Section  C  contains  a  set  of  Footnotes  and 

Comments  which  should  be  read  in  association  with  the  main  text. 

C.l.  It  is  important  to  understand  that  the  ICARUS  1  System  (2)  [completed  in 
December  1993]  used  a  Keithley  199  scanning  DMM  unit  and  that  this  system 
introduced  timing  shifts  into  the  data  acquisition  processes  [the  well-known 
multiuser-multitasking  problem  associated  with  the  use  of  PC's].  This  applies 
equally  to  all  other  experiments  carried  out  before  ~  1st  January  1995.  These 
timing  errors  did  not  affect  evaluations  based  on  (kR')i,  (kR')2,  (kR')n  or  (kR')i2  but 
had  to  be  taken  into  account  in  some  evaluations  at  higher  levels  of  precision  and 
accuracy  based  on  (kR')2i  and  (kR')22;  they  were  especially  important  in  the 
evaluations  of  (kR')3i  and  (kR')32  as  was  explained  especially  in  three  reports 
(3),(4),(5)  |jor  the  definition  and  use  of  these  heat  transfer  coefficients,  see  Sections 
A  and  B]. 

There  is  also  extensive  further  correspondence  dealing  with  the  analyses  of  Data 
Sets  collected  with  the  ICARUS- 1  and  2  Systems.  However,  the  three  reports 
(3).(4).(5)  are  sufficient  to  explain  the  reasons  for  the  choices  of  the  methodologies. 

C.2  The  two  Reports  (3)(4)  were  sent  in  the  first  place  to  Technova  [see  also  Footnote 
D.2], 

C.3.  It  should  be  noted  that  it  is  possible  to  achieve  accurate  control  of  the  temperature 
of  the  heat  sink  by  using  this  simple  design.  In  turn  this  allowed  the  installation 
of  a  large  number  of  experiments  at  low  cost. 

There  are  misleading  statements  in  the  literature  about  this  aspect  of  the 
experiment  design.  In  fact,  the  design  follows  the  common  strategy  of  using  two 
thermal  impedances  in  series,  a  strategy,  which  is  required  for  experiments  aiming 
at  high  accuracy,  [however,  see  also  Section  A.2]. 

C.4.  The  major  reason  is  that  this  is  the  only  experiment  to  date  where  we  have  a 
listing  of  the  raw  data  together  with  the  evaluation  carried  out  according  to  the 
procedures  used  in  the  NHE  Laboratory.  This  is  evidently  of  key  importance  to 
(i)  and  (ii)  especially  in  view  of  the  assessment  of  the  ICARUS  methodology 
which  has  been  published  (ll)  [see  also  Footnote  D.3]. 

C.5.  (via)  several  "blank  experiments"  were  to  be  carried  out  [at  least  one  for  each 

cell  in  use;  the  use  of  Pt-cathodes  in  0. 1 M  LiOD/D20  was  recommended]. 
The  experiments  were  to  be  carried  out  at  a  single  current  density  using 
otherwise  the  protocol  described  in  (ia)-(va); 

(viia)  the  execution  of  the  "blank  experiments"  was  to  be  followed  by 
experiments  using  cathodes  made  of  Johnson  Matthey  Material,  Type  A. 
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The  protocols  (ia)  -  (va)  were  again  to  be  followed; 

(viiia)  the  first  step  in  the  data  evaluation  was  to  be  the  plotting  of  A3  or  A4 
sized  graphs  of  the  "raw  data".  The  heat  transfer  coefficients  (kR')i  and 
(kR')2  were  to  be  derived  for  each  measurement  cycle  [see  Section  B], 

(ixa)  the  next  step  was  to  be  the  construction  of  (kR')n  -  type  spreadsheets 
coupled  to  the  determination  and  interpretation  of  (kR')n,  (kR  ')„  and 

(kg'Jn  [see  Section  B.3],  The  further  evaluation  of  these  spreadsheets 
was  not  specified  in  1993;  this  was  a  matter  which  was  to  be  decided  by  a 
collaborative  program  between  NHE  and  I. M  R. A.  Europe. 

(xa)  after  the  execution  of  (viiia)  and  (ixa)  the  (kR')2i  -  type  spreadsheets  were 
to  be  prepared  and  values  of  (kR'o)26i  and  (k^o^i  and  the  associated 
values  of  CPM  were  to  be  determined.  These  values  of  CPM  were  to  be 
used  both  to  correct  the  evaluations  in  (ixa)  and  to  determine  the  "true  heat 
transfer  coefficients,  (kR')32"  at  times  close  to  the  end  of  the  calibration 
period.  Following  the  receipt  of  the  first  set  of  "raw  data"  collected  with 
the  ICARUS-1  system  (3)  this  objective  was  changed  to  the  evaluation  of 
(kR’)22  and  of  (kR'o)262  at  times  close  to  the  end  of  the  calibration  period. 

(xia)  it  was  envisaged  that,  following  the  completion  of  this  initial  stage  of  the 
investigation,  the  ICARUS  program  would  move  on  to  the  examination  of 
materials  variables  as  well  as  the  production  of  ICARUS-2.  This  second 
part  was  intended  to  deal  also  with  the  effects  of  increasing  the  current 
density,  the  analysis  of  data  close  to  the  boiling  point  [including  boiling 
episodes]  and  "Heat-after-Death".  The  effects  of  changes  of  current 
density  were  to  be  followed  in  line  with  earlier  investigations  (15),(16),  see 
(iiia),  using  an  initial  polarization  at  low  current  densities  for  various 
periods  of  time  followed  by  a  raising  of  the  current  density  sufficient  to 
drive  the  systems  through  the  region  of  "positive  feedback"  and, 
eventually,  towards  the  boiling  point.  The  analysis  of  data  close  to  the 
boiling  point  and  of  "Heat-after-Death"  was  to  be  covered  initially  by 
appropriate  "Case  Studies".  The  material  given  in  Section  A7  -  9  goes 
some  way  towards  fulfilling  the  needs  of  such  "Case  Studies". 

The  topic  of  "positive  feedback"  was  not  included  in  the  ICARUS-1  program  mainly 
because  these  effects  greatly  complicate  the  analyses  [see  Section  A. 3  and  A. 4]  but  also 
because  aspects  of  this  topic  were  still  being  actively  investigated  at  IMRA  Europe  [see 
Footnote  E.l  1],  However,  extensive  comments  on  this  topic  were  included  in  the  Report 
(4)  on  the  second  set  of  experiments  carried  out  by  NHE  as  the  systems  had  evidently  been 
driven  into  the  regions  of  "positive  feedback"  [see  also  <8)], 

C.6.  However,  the  material  outlined  in  <l3)  was  also  contained  in  earlier  letters  and 
reports  sent  to  Japan  and,  in  part,  in  papers,  which  have  been  published 
previously. 
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C.7.  Reference  13  deals  first  of  all  with  the  analysis  of  a  data  set  generated  by 
simulation  using  the  most  the  drastic  simplifications  of  the  differential  equation 
representing  these  Fleischniann-Pons  calorimeters.  The  second  part  deals  with  the 
analysis  of  a  blank  experiment  [Pt  cathode]  polarized  in  0.1  M  Li0D/D20.  This 
full  report  in  as  received  condition  [from  Fleischmann]  consists  of  190  pages  and 
will  be  published  in  another  report  (see  Ref.  1 3). 

C.8.  We  note  here  that  the  values  of  the  rates  of  evaporative  loss  given  in  the 
spreadsheets  contained  in  the  Report  have  been  calculated  with  an  earlier  version 
of  equation  (B.4)  using  the  values 

L  =  41,673.7  J  Mole'1 

F  =  96,484  C  Mole'1 

At  temperatures  close  to  the  boiling  point  the  values  of  P*  and  Oboiiing  were 
adjusted  as  described  in  Section  B.2. 

C.9.  Changes  in  the  parameter  listing  are  much  more  important  for  other  parts  of  the 
evaluation  e.g.  the  adjustment  of  the  heavy  water  equivalent,  CpM  of  the  cells  in 
the  evaluation  of  the  lower  bound  heat  transfer  coefficients,  (kR’)u  [see  Section 
A.3], 

C.10.  Independent  calibrations  show  that  the  heat  transfer  coefficient  for  cells  filled 
with  air  are  about  0.75  of  the  values  of  these  coefficients  for  cells  filled  with 
liquid  *9^l0l  It  follows  therefore  that  the  marked  increase  in  the  heat  capacitance 
of  the  cells  filled  with  D20  vapor  at  temperatures  close  to  the  boiling  point  must 
lead  to  the  maintenance  of  the  heat  transfer  coefficient  at  the  value  which  applies 
to  cells  filled  with  liquid. 

C.ll.  The  group  at  the  NHE  Laboratory  attempted  to  determine  the  values  of  AM 
directly  by  adding  a  condensation  section  to  the  cells  [see  Footnote  D.5]. 

C.12.  As  has  been  explained  previously  (l3),  we  had  hoped  to  circumvent  the  need  for 
such  an  extended  description  so  as  to  avoid  overburdening  our  accounts  with 
redundant  symbolism.  In  fact,  the  Handbooks  *l3(2)and  earlier  Reports* 
used  a  simplified  description.  Again,  as  has  been  described  previously  (l3),  it  is 
the  deviation  by  the  NHE  group  from  the  recommended  protocols  and  procedures 
and  their  use  of  inaccurate  and/or  invalid  methods  of  data  evaluation  (ll)  which 
makes  it  necessary  to  use  the  more  closely  defined  description,  (B.8). 

C.13.  Of  course,  other  averages  can  be  made  but  the  use  of  the  1 1 -point  average  has 
been  found  to  be  especially  useful. 

C.  14.  The  value  of  CPM  =  490  JK'1  was  appropriate  for  somewhat  larger  cells  used  in  an 
earlier  phase  of  our  investigations. 
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C.  1 5.  Objections  have  often  been  raised  to  the  procedures  which  we  have  adopted  based 
on  the  fact  that  it  is  alleged  that  we  have  not  "binned  the  data"  i.e.  we  have  not 
signal  averaged  before  the  data  analysis  [see  Footnote  D.7],  However,  "binning 
of  the  data"  must  always  be  approached  with  great  caution:  one  should  only  "bin 
data"  or  "bin  coefficients"  if  these  data  or  coefficients  can  be  expected  to  be 
constant  over  the  averaging  interval.  This  is  not  the  case  for  (kR')n  unless  the 
effects  of  the  term  CpM(dA0/dt)  have  been  taken  into  account.  Once  this  is  done 
we  can,  of  course,  bin  the  coefficients  as  we  have  done  in  deriving  (k/)^  and 

cv 

In  the  case  of  the  interpretation  of  data  obtained  with  calorimeters  relying  on 
radiative  cooling,  the  position  is  further  complicated  by  the  fact  that  the 
differential  equation  representing  the  calorimeters,  equation  (B.l),  is  both  non¬ 
linear  and  inhomogeneous.  It  does  not  follow  therefore  that  the  coefficients 
derived  by  averaging  the  data  will  be  the  same  as  the  averages  of  the  coefficients 
derived  by  using  the  "raw  data".  However,  it  was  confirmed  in  1992  [when  the 
approach  outlined  here  was  initiated]  that  one  does,  in  fact,  obtain  an  equivalence 
provided  attention  is  restricted  to  regions  where  (dA0/dt)  is  adequately  small.  It 
was  concluded  that  in  such  regions  the  differential  equation  (B.l)  could  be 
sufficiently  linearised  and  that  the  second  order  small  differences  were 
sufficiently  small  to  allow  such  averaging.  However,  as  the  values  of  (kR')i  j 
obtained  in  this  was  were  identical  to  those  obtained  using  the  procedure  outlined 
in  Section  B.3,  there  was  evidently  no  justification  in  pursuing  the  matter  further 
[nor  to  complicate  the  instrumentation!]. 

It  is  evident  therefore  that  these  particular  criticisms  were  invalid  [see  Footnote 
E.4], 

Averaging  procedures  are  considered  further  in  Section  B.3  in  the  discussion  of 
heat  transfer  coefficients  based  on  the  forward  and  backward  integration 
procedures. 

C.  16.  A  similar  comment  applies  to  the  determination  of  (kR'o)i6i;  the  time  at  which 
(dA0/dt)  =  0  will  usually  be  accessible  for  experiments  in  which  ft  =  9  hours. 
However,  no  such  point  can  be  defined  for  (kR'o)i7i  so  that  this  determination  is 
mathematically  questionable;  this  is  therefore  equally  true  for  (kR'o)i«i  although 
these  extrapolations  are  certainly  sound  from  an  operational  point  of  view. 

C.l  7.  The  validity  of  equation  (B.l  7)  was  established  at  the  time  of  construction  of  the 
ICARUS-2  System. 

C.l 8.  In  our  early  work  [prior  to  October  1989]  we  had  overcome  this  particular 
difficulty  by  using  (kR')2  as  starting  values  for  the  non-linear  regression  procedure 
leading  to  the  "true  heat  transfer  coefficient,  (kR‘)?".  Here  we  fitted  the  numerical 
integrals  of  the  differential  equation  (B.l)  to  the  data  sets  for  complete 
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measurement  cycles.  The  relative  errors  in  (kR')5  which  we  could  achieve  in  this 
way  were  below  0.1%  even  when  using  the  unsatisfactory  early  version  of  our 
calorimeters  [i.e.  those  not  silvered  in  their  top  portions]. 

This  use  of  the  non-linear  regression  procedure  had  the  further  distinct  advantage 
that  we  could  adjust  the  integration  interval  in  regions  where  the  temperature  was 
varying  rapidly  with  time  so  as  to  achieve  the  required  accuracy  in  the  integrals. 
This  was  not  possible  for  the  methods  leading  to  (kR')2i,  (kR')22,  (kR’)3i  and  (kR’)32 
discussed  in  Section  B.3  when  using  the  ICARUS-1  and  -2  Systems  because  the 
intervals  for  data  acquisition  were  fixed.  Indeed,  the  interval  300  s  was  chosen 
because  such  an  interval  does  not  degrade  the  evaluation  of  any  of  the  series 
(kR’)2i,  (kR')22,  (kR'o)25i,  (kR'o)252,  (kR'o)26i,  (kR'o)262,  (kR'o)27i  and  (kR'o)272. 
However,  it  does  degrade  the  evaluation  of  (kR*)3i,  (kR'o)35i,  (kR'o)36i  and  (kR'o)37i 
to  some  extent  and  leads  to  a  marked  degradation  of  (kR')32,  (kR'o)352,  (kR’o)362, 
and  (kR'o)372.  The  fact  that  the  data  acquisition  interval  was  too  long  for 
straightforward  estimates  of  the  (kR'),j  i  series  of  heat  transfer  coefficients  was 
already  pointed  out  to  NHE  in  the  ICARUS- 1  Handbook  ^  ^  [see  Footnote  E.6]. 

We  have  pointed  out  elsewhere  (15)  (16)that  a  major  reason  for  our  choice  of  non¬ 
linear  regression  was  that  the  pressure  of  events  in  1989  did  not  allow  us  to  go 
through  the  logical  sequence  of  using  linear  regression,  multi-linear  regression 
and  non-linear  regression  [in  fact,  we  had  to  opt  for  a  "catch-all"  methodology 
and  any  integrated  version  of  equation  (B.l)  was  certainly  expected  to  require  the 
use  of  non-linear  regression].  However,  our  use  of  this  methodology  was 
certainly  not  understood  [see  Footnote  E.7]  and  we  also  could  not  make 

non-linear  regression  "user  friendly"  with  the  computing  power  available  to  us  at 
that  time.  In  1991-92  we  therefore  investigated  the  application  of  linear  regression 
(i5).(i<>)  thjs  j,ecame  part  and  parcel  of  the  ICARUS-1  methodology  (2). 

C.19.  It  is  in  any  event  necessary  to  change  the  methodology  of  the  evaluation  in  view 
of  the  early  onset  of  "positive  feedback"  see  Section  A.4. 

C.20.  As  far  as  one  is  able  to  tell  from  the  limited  amount  of  information  available. 

C.21.  We  retained  the  designation  22  as  a  flag  to  indicate  that  the  "backward 
integration"  methodology  was  the  primary  objective  for  accurate  evaluations. 

C.22.  Reference  (13)  contains  extensive  discussions  of  the  errors  of  the  various  heat 
transfer  coefficients  and  the  causes  of  these  errors.  This  reference  will  be 
published  in  another  report. 

C.23.  This  will  explain  both  our  strategies  for  determining  the  heat  transfer  coefficients 
[which  give  the  values  at  t  =  t2]  as  well  as  giving  a  further  reason  for  choosing  2- 
day  measurement  cycles  with  t2  corresponding  to  the  end  of  Day  1. 

C.24.  The  need  for  posing  this  restriction  was  pointed  out  repeatedly.  However,  we 
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note  that  the  group  at  NHE  used  integral  heat  transfer  coefficients  to  evaluate  the 
rates  of  excess  enthalpy  generation  (ll)  and  this  applies  equally  to  their 
spreadsheets  for  the  experiment  considered  in  this  Report. 

C.25.  Or,  at  least,  they  were  not  followed  systematically. 

C.26.  The  points  at  which  data  were  acquired  were  always  correctly  timed.  This  is  true 
for  all  the  experiments,  which  we  have  carried  out  [see  Footnote  E.8], 

C.27.  and  as  has  been  pointed  out  elsewhere,  we  no  longer  have  the  "raw  data"  for  these 
experiments. 

C.28.  It  is  therefore  impossible  to  carry  out  the  tests  based  on  the  application  of  the 
protocols  specified  for  these  experiments,  especially  those  based  on  (iia)  see 
-  Section  A.2. 

C.29.  However,  we  should  note  that  as  far  as  the  results  for  this  experiment  are 
concerned,  the  early  development  of  "positive  feedback"  would  in  any  event  have 
precluded  the  evaluation  of  (kR')2  [see  Sections  A. 3  and  4], 

C.30.  We  should  make  a  special  note  here  of  the  fact  that  the  analyses  presented  in  this 
report  are  incomplete  in  that  the  effects  of  these  changes  in  current  density  have 
not  been  considered  except  for  that  on  Day  26  which  has  been  discussed  as  a 
special  case  of  "Heat-after-Death",  see  Section  A.  10.  In  fact,  it  is  possible  to 
develop  analyses  analogous  to  those  leading  to  (kR')22,  (kR')262,  (kR')32  and  (kR')362 
but  based  on  the  perturbations  due  to  changes  of  current  density.  This  is  a  matter, 
which  should  be  considered  further  in  an  addendum  to  this  Report. 

In  this  connection  it  is  of  some  interest  that  the  analysis  of  the  effects  of  changes 
of  current  density  in  one  of  the  data  sets  of  the  investigation  by  the  group  in 
Harwell (22)  gives  a  clear  demonstration  of  excess  enthalpy  generation.  A  paper 
dealing  with  this  topic  (23)  was  rejected  for  publication  [two  further  papers  dealing 
with  the  reanalysis  of  data  sets  generated  by  this  investigation  have  also  been 
rejected  for  publication]. 

C.31.  A  phenomenon  which  has  also  been  observed  and  commented  on  by  other 
authors. 

C.32.  It  has  been  argued  (6)(7)  that  this  phenomenon  is  linked  to  the  need  to  achieve  high 
levels  of  loading  of  the  cathodes  by  D+  which  is  probably  also  associated  with  a 
change  from  exo-  to  endothermic  absorption;  an  alternative  explanation  (24)  is  that 
these  phenomena  are  linked  to  the  formation  of  a  new  phase,  the  y-phase (25). 

It  is  not  clear  at  the  present  time  whether  [and  if  so  how]  the  onset  of  "positive 
feedback"  is  linked  to  the  increase  of  cell  temperature  and  the  results  of  the 
experiment  certainly  argue  against  such  an  interpretation  [see  Sections  A. 3  and 
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A. 4],  However,  it  seems  to  be  reasonably  clear  that  increases  of  temperature 
facilitate  the  phenomenon.  Furthermore,  maintenance  of  the  electrode  in  the 
region  of  the  onset  of  "positive  feedback"  [current  density,  temperature]  leads  to 
the  onset  of  oscillatory  phenomena (6)  (7)  (24)  which  reduce  and  may  even  eliminate 
the  effect.  It  is  the  sensitivity  of  "positive  feedback"  to  the  operating  conditions, 
which  is  a  major  factor  in  making  the  outcome  of  experiments  [i.e.  the  level  of  the 
rates  of  excess  enthalpy  generation]  dependent  on  the  experimental  protocols. 

C.33.  We  have  to  take  note  of  the  fact  that  the  usable  range  of  current  densities  was 
necessarily  restricted  by  the  maximum  output  of  the  galvanostat  [±  1A  at  ±  100 
V]  and  the  relatively  large  surface  area  of  the  electrode  [3.15  cm2]. 

C.34.  Furthermore,  the  current  densities  remained  in  the  low  range  for  the  observation 
of  excess  enthalpy  generation (12). 

C.35.  A  correction  to  Fig.  A.4  may  form  an  addendum  to  this  Report. 

C.36.  The  effects  of  such  changes  may  form  a  further  addendum  to  this  Report. 

C.37.  Which  furthermore  quoted  results  showing  that  the  volumes  of  the  electrolyzed 
gases  agreed  with  those  calculated  for  electrolysis  according  to  Faraday’s  Laws. 

C.38.  Nevertheless,  the  recombination  of  D2  and  02  formed  in  the  electrolysis  continues 
to  be  invoked  in  explanations  of  the  excess  enthalpy  generation  observed  e.g. 
and  the  numerous  comments  by  inter  alia  D.  Morrison  and  S.  Jones. 

C.39.  The  conditions  of  operation  above  Day  62  need  to  be  investigated  further. 

C.40.  Hypothetical  differences  in  temperature  due  to  inadequate  mixing  have  frequently 
been  invoked  in  hypothetical  arguments  about  the  performance  of  these 
calorimeters.  We  note,  however,  that  inadequate  mixing  would  not  give  rise  to  a 
systematic  and  time  invariant  difference  in  temperature  between  the  two 
positions.  Moreover,  such  differences  in  temperature  would  not  be  expected 
because  the  thermal  relaxation  time,  x  =  CpM/4(kR')e3ceii,  is  of  order  5000  s 
whereas  the  radial  and  axial  mixing  times  are  ~  3  and  ~  20  s  [as  revealed  by  tracer 
experiments].  Small  differences  in  the  cell  temperature  can  only  be  observed  in 
the  vicinity  of  the  electrodes  and  calibration  resistor  i.e.  within  the  Prandtl 
boundary  layers.  However,  the  volumes  of  liquid  within  these  boundary  layers 
are  negligibly  small  compared  to  the  volume  of  electrolyte  in  the  cells. 

C.41.  In  considering  the  precision  and  accuracy  of  the  measurements  of  the  "raw  data" 
we  should  note  that  the  cell  currents,  column  6  of  Table  A.2,  have  been  rounded 
to  5  significant  figures.  These  currents  should,  in  fact,  have  been  given  to  6 
significant  figures. 

C.42.  With  a  thermal  relaxation  time  of  5000  s,  the  temperature  perturbation  will  only 
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have  reached  98.67%  of  its  final  value  within  the  calibration  period.  The 
consequent  1.33%  error  is  outside  the  limits  specified  for  the  ICARUS  Systems 
[see  also  Section  A. 4], 

C.43.  Such  effects  can  be  seen  in  the  "raw  data"  of  some  of  the  experiments  carried  out 
by  the  group  at  Harwell  <22)as  well  as  in  some  of  the  early  measurements  carried 
out  by  the  group  at  NHE (4)  (5)  (24) 

C.44.  This  is  a  feature,  which  will  be  common  to  all  calorimetric  systems  used  to 
investigate  a  thermal  source  subject  to  "positive  feedback".  [Moreover  it  is  a 
feature  of  any  system  subject  to  "positive  feedback"].  It  is  likely  that  the  neglect 
of  this  fact  is  responsible  for  much  of  the  confusion  in  the  research  on  "Cold 
Fusion". 

C.45.  The  spreadsheets  contain  the  following  columns  of  entries: 

1 )  the  time  of  the  data  acquisition,  t/s; 

2)  Occii  short  thermistor,  the  temperature  measured  with  the  "short  thermistor"/K, 
see  Fig.  A.l ; 

2)  0Ccii.  long  thermistor,  the  temperature  measured  with  the  "long  thermistor"/K, 
see  Fig.  A.l; 

4)  Eeeu/V; 

5)  Iccii/A, 

6)  Obath,  the  temperature  of  the  water  bath,  /K; 

7)  (de  cell/ dt)short  thermistor?  /Ks’1; 

8)  the  enthalpy  input,  /W; 

9)  the  rate  of  evaporative  cooling,  AV; 

1 0)  the  Joule  heat  input  due  to  the  resistive  heater,  AV; 

11)  f]  (O)short  thcrmislorAC 

12)  the  rate  of  excess  enthalpy  generation,  AV; 

13)  109  (kR')n,  AVK"4; 

14)  an  entry  related  to  the  water  equivalent  of  the  calorimeter,  /JK'1; 

15)  (d0/dt)long  thermistor,  AC  S  , 

1 6)  the  rate  of  evaporative  cooling,  AV; 

17)  f]  (O)long  thermistor,  /K  ; 

1 8)  the  rate  of  excess  enthalpy  generation,  AV; 

19)  109  (kR*),,  AVIC4; 

20)  an  entry  related  to  the  water  equivalent  of  the  calorimeter,  /JK'1. 

It  is  evident  that  the  derived  data  in  columns  (7),  (9),  (1 1),  (12)  and  (13)  are  based 
on  the  temperature  measured  with  the  "short  thermistor",  column  (2),  while  those 
in  columns  (15),  (16),  (17),  (18)  and  (19)  are  based  on  that  measured  with  the 
"long  thermistor",  column  3.  The  basis  of  the  entries  in  columns  (14)  and  (20)  is 
not  clear,  nor  is  it  apparent  how  this  information  may  have  been  used. 

C.46.  The  water  thermostats  surrounding  the  cells  were  run  at  -  30°C  in  the  early  work 


50 


in  Salt  Lake  City.  In  1988  attempts  were  made  to  allow  for  this  shift  in  the 
reference  temperature  as  well  as  the  fact  that  electrolysis  takes  place  Irom  0.1  M 
LiOD  in  D20  and  not  D20  itself.  While  the  thermoneutral  potential  is  certainly 
not  1.527  V,  it  is  closer  to  this  value  than  to  1.54  V.  The  value  1.527  V  has 
therefore  been  used  in  the  calculations  summarized  in  Tables  A.3  -  A.9. 

C.47.  It  appears  that  the  facilities  for  monitoring  the  current  to  the  calibration  resistor 
and  the  voltage  across  this  resistor  were  disabled. 

C.48.  More  exactly,  whether  the  values  of  (kR')i  i  plotted  versus  time  fall  on  a  common 
straight  line  as  shown  in  e.g.  see  also  Fig.  A.20  in  Section  A. 5. 

C.49.  In  actual  fact,  the  matter  is  more  complicated  than  has  been  indicated  in  the  main 
text.  In  the  method  originally  proposed  by  NHE  the  "lower  bound  heat  transfer 
coefficient"  determined  before  the  application  of  the  calibration  pulse  was  used  in 
attempts  to  derive  the  rate  of  excess  enthalpy  generation  during  the  application  of 
this  pulse.  It  is  not  surprising  that  such  a  method  can  only  give  the  correct  result 
provided  there  is  a  zero  rate  of  excess  enthalpy  generation  for  the  period  t  <  ti 
before  the  application  of  the  pulse  [as  well  as  for  t  >  t2  following  the  termination 
of  the  pulse]  while  the  calibration  pulse  itself  [during  t,  <  t  <  t2]  leads  to  the 
generation  of  excess  enthalpy.  These  facts  were  illustrated  using  simulated  data 
but  the  Report  did  not  lead  to  any  further  discussion. 

It  is  of  interest  that  the  behavior  of  this  Pd-B  study  during  Day  3  to  some  extent 
satisfies  the  preconditions  for  the  application  of  the  NHE  methodology  [see  the 
main  text].  It  is  also  of  interest  that  the  method  proposed  by  NHE  bears  some 
resemblance  to  an  evaluation  used  previously^  ^  ^ except  that  the  eailier  method 
used  the  maximum  value  of  the  "lower  bound  heat  transfer  coefficient,  (kR  )n 
observed  during  the  early  stages  of  the  experiment.  The  use  of  such  a  value  of  the 
"lower  bound  heat  transfer  coefficient"  is  discussed  further  in  Section  A. 6. 

C.50.  It  so  happens  that  there  is  some  validity  to  this  conclusion  due  to  the  influence  of 
"positive  feedback".  This  matter  is  discussed  in  Section  A. 4. 

C.51.  It  is  evidently  very  important  to  establish  whether  this  early  onset  of  "positive 
feedback"  is  a  property  of  the  Pd-B  alloy  used  in  this  experiment.  It  should  be 
noted  that  it  was  proposed  some  years  ago  that  Pd-Ag  alloys  used  in  diffusion 
tube  separators  should  be  replaced  by  Pd-B  alloys.  Evidently  Pd-B  alloys  were 
expected  to  display  some  beneficial  properties,  and  the  literature  on  this  topic 
should  be  searched  and  the  main  manufacturers  should  be  consulted. 

C.52.  The  steps  outlined  in  Section  A. 4  do  not  correspond  to  the  sequence  in  the 
ICARUS-style  evaluation  (2).  In  that  evaluation  [as  modified  in  (3)j  the  first  step 
is  the  preparation  and  inspection  of  A-4  sized  plots  of  the  "raw  data"  leading  to 
the  evaluation  of  (kR’)i  and  (kR')2;  the  second  step  is  the  preparation  and 
interpretation  of  (kR')n  -  spreadsheets  while  the  third  step  is  the  preparation  and 
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evaluation  of  (kR,)2i/(kR,)3t  -  spreadsheets  leading  up  to  the  evaluation  of  (kR'o)26i, 
(kR'W,  (kR')36i  and  (kR1^.  It  should  also  be  noted  that  the  ICARUS  software 
prepares  the  more  restricted  version  of  the  (kR')2i/(kR')3i  -  spreadsheets  [restricted 
to  ti<  t  <t2]  as  shown  in  Table  A.5  and  not  the  version  for  the  whole  measurement 
cycle  0  <  t  <  T  shown  in  Table  A.4. 

It  should  be  noted  furthermore  that  the  tables  presented  in  this  Report  only 
contain  a  part  of  the  information  of  the  ICARUS-style  spreadsheets. 

C.53.  As  has  been  noted  in  Section  B.3,  it  is  possible  that  this  insertion  of  additional 
data  points  has  not  been  carried  out  in  the  systems  supplied  to  NHE.  This 
particular  question  could  only  be  resolved  by  re-implementing  the  relevant 
software. 

C.54.  The  values  for  the  first  20  -  30  points  must  be  excluded  as  the  benefits  of  using 
the  integral  coefficients  are  only  established  with  increasing  time.  Similarly,  the 
first  20  -  30  points  must  be  excluded  if  the  interpretation  is  based  on  backward 
integration  i.e.  if  we  consider  (kR')2i. 

C.55.  More  exactly  the  evaluation  of  (kR1^  was  restricted  to  times  close  to  t2  for 
calibration  periods  of  12  hours  i.e.  for  measurement  cycles  lasting  2  days  (2). 

C.56.  The  increase  of  the  pseudoradiative  heat  transfer  coefficient  due  to  the  "over 
filling"  of  the  cells  is  important  from  a  different  point  of  view.  It  has  been  shown 
in  Section  B.2  that  the  thermal  capacity  of  the  head  space  in  the  cells  [filled  with 
the  mixture  of  deuterium,  oxygen  and  D20  vapor]  must  increase  markedly  above 
that  predicted  for  the  head  space  filled  with  deuterium  and  oxygen  alone  at 
temperatures  close  to  the  boiling  point  of  the  electrolyte.  It  follows,  therefore, 
that  heat  transfer  from  the  cell  must  be  controlled  by  the  coefficient  for  cells 
"overfilled"  with  electrolyte.  This  matter  is  considered  further  in  Section  A. 6. 

C.57.  The  main  question  is  that  if  the  charging  of  the  electrode  is  indeed  100%  efficient 
and,  if  this  charging  is  completed  at  ~  130,500  s,  then  the  total  charge  taken  up 
by  the  electrode  is  in  excess  by  a  factor  of  five  of  that  predicted  for  a  charging 
ratio  D/Pd  =  1 .  Related  to  this  question  is  the  problem  that  the  total  excess 
enthalpy  during  the  charging  period  [~  3800  J]  is  in  excess  by  a  factor  of  two  of 
that  which  would  be  predicted  for  D/Pd  =  1  and  an  heat  of  absorption  of  40  kJ 
mole1.  This  suggests  anomalous  excess  enthalpy  production  during  the  charging 
period.  Experimentally,  vigorous  gas  evolution  was  observed  after  13,200  s 
[1=0. 1 50A],  This  corresponds  to  a  loading  of  D/Pd=0.5  at  that  time. 

C.58.  Evidently,  heat  transfer  to  the  walls  of  the  Dewar  cell  is  maintained  by  the  vapor 
phase  at  the  very  least  if  this  phase  is  filled  with  D20  vapor  at  temperatures  close 
to  the  boiling  point  of  the  electrolyte  [c.f.  the  increase  in  the  thermal  capacity  of 
the  vapor  phase.  Section  B.2],  It  follows  that  the  heat  transfer  coefficient  for  the 
cell  filled  with  vapor  is  the  same  as  the  heat  transfer  coefficient  for  the  cell  filled 
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with  electrolyte.  In  actual  fact,  the  heat  transfer  coefficient  for  the  cell  filled  with 
vapor  will  be  that  for  a  cell  "overfilled"  with  electrolyte  i.e.  some  5%  above  the 
value  0.85065  x  10'9  WK'4  [c.f.  the  discussion  of  Fig.  A.20], 

C.59.  It  is  important  to  realize  that  the  calculation  of  the  total  enthalpies  outlined  in  the 
main  text.  Section  A. 7,  is  independent  of  the  rate  of  reflux. 

C.60.  It  is  also  important  to  realize  that  the  interpretation  given  by  NHE  shows  the  other 
features  of  excess  enthalpy  generation  to  which  we  have  frequently  drawn 
attention  e  g.  the  fall  of  the  "lower  bound  heat  transfer  coefficient,  (kR')n",  with 
time  due  to  the  build  up  of  excess  enthalpy  production.  In  fact  excess  enthalpy 
generation  was  observed  throughout  the  time  range  except  for  the  last  four  data 
points  where  the  calculation  of  the  rates  of  evaporative  cooling  are  invalid. 
These  matters  passed  without  comment. 

C.61 .  The  group  at  NHE  attempted  to  derive  the  rate  of  evaporative  cooling  directly  by 
collecting  and  weighing  the  distillate.  However,  apart  from  the  inaccuracies 
introduced  by  the  unknown  hold-up  in  the  distillation  system,  the  data  derived 
cannot  be  interpreted  in  the  absence  of  information  about  the  rate  of  reflux.  It  is 
therefore  in  any  event  necessary  to  carry  out  the  calculation  outlined  in  Section 
A. 7  so  that  the  weighing  system  does  not  produce  any  additional  useful 
information. 

C.62.  The  operation  of  the  experiment  on  Day  68  should  be  contrasted  with  earlier 
measurements  rates  of  excess  enthalpy  generation  in  these  earlier 

experiments  reached  much  higher  values  than  those  reported  here  [up  to  two 
orders  higher].  Some  of  the  complications  in  the  evaluation  of  this  experiment 
could  therefore  be  avoided.  In  particular,  it  could  be  assumed  that  the  last  part  of 
the  earlier  experiments  were  controlled  by  a  process  akin  to  boiling.  This  was 
sufficiently  intense  that  the  effects  of  reflux  could  be  neglected. 

C.63.  As  has  been  noted  above  the  calibrations  used  in  an  earlier  investigation  were 
derived  by  using  a  heater  spiral  spanning  the  whole  volume  of  the  cell  i.e.  heat 
was  applied  uniformly  throughout  this  volume. 

C.64.  However,  if  we  describe  this  as  the  "conventional  wisdom"  of  the  research  in  this 
field,  then  we  can  see  that  the  results  for  this  experiment  fly  in  the  face  of  this 
"conventional  wisdom".  In  particular,  excess  enthalpy  generation  was  observed 
on  Day  3  of  the  measurement  cycles  at  a  current  density  below  the  threshold  value 
while  "positive  feedback"  was  established  at  a  temperature  below  this  further 
threshold.  We  can  therefore  only  regard  the  criteria  used  to  search  for  category 
#1  of  the  phenomenon  of  "Heat  after  Death"  as  rather  "broad  brush  indicators". 

C.65.  Possibly  because  of  the  cracking  of  the  electrodes  due  to  the  repeated  loading  and 
de-loading. 
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C.66.  The  non-uniform  charging  of  electrodes  is  probably  responsible  for  the  marked 
increase  in  the  rate  of  excess  enthalpy  generation  at  temperatures  approaching  the 
boiling  point.  The  electrodes  then  become  covered  by  a  porous  film  and  charging 
takes  place  at  very  high  current  densities  through  the  pores  in  the  film.  This 
process  is  favored  by  the  "tertiary  current  distribution"  in  the  film:  the  very  high 
current  densities  lead  to  increases  in  the  OD'  concentration  in  the  pores  which  in 
turn  enhances  the  non-uniform  current  distribution. 

A  program  of  work  on  the  non-uniform  charging  of  electrodes  [using  inter  alia 
electrodes  partially  covered  by  diffusion  barriers]  could  not  be  established  at 
IMRA  Europe  because  of  lack  of  resources.  This  is  probably  one  of  the  most 
important  neglected  aspects  of  the  work  because,  if  such  charging  should  prove  to 
be  successful,  it  would  allow  the  construction  of  energy  efficient  systems. 

C.67.  It  is  unlikely  that  the  variation  of  the  distillate  with  time  [as  determined  in  the 
NHE  investigations]  could  be  usefully  interpreted. 


Section  D 

As  is  explained  in  the  main  text,  Section  A,  Section  D  contains  a  set  of  Footnotes  and 

Comments,  which  should  not  necessarily  be  included  in  a  Report  given  extensive 

circulation. 

D.l.  [see  also  Footnote  C.l]  As  has  been  explained  in  earlier  correspondence,  M. 
Fleischmann  only  obtained  a  copy  of  the  Handbook  describing  the  ICARUS-2 
System  (I)  indirectly  during  1998  from  M.H.  Miles.  Prior  to  this,  M.  Fleischmann 
was  denied  access  to  this  Document. 

D.2.  [see  also  Footnote  C.2]  It  is  not  clear  whether  NHE  ever  took  any  account  of  these 
reports  [or  of  the  extensive  further  correspondence]  as  their  receipt  was  never 
acknowledged.  Therefore,  there  was  no  correspondence  concerning  their 
contents.  However,  we  note  that  the  whole  project  was  to  be  carried  out  in  the 
Public  Domain.  Presumably,  therefore,  these  Reports  should  be  made  available 
to  "third  parties"  if  this  should  be  required. 

D.3.  [see  also  Footnote  C.4]  At  the  7th  International  Conference  on  Cold  Fusion,  one 
of  us  [M.  Fleischmann]  was  also  given  a  CD  containing  the  raw  data  for  a  series 
of  experiments  carried  out  in  the  NHE  Laboratories.  This  CD  may  also  contain 
information  about  the  evaluations,  which  were  carried  out.  However,  none  of 
these  experiments  correspond  to  the  ones,  which  had  originally  been  requested 
[see  (5)and  the  preceding  reports  and  correspondence].  These  experiments  were  all 
carried  out  prior  to  June  1994  [although  some  of  the  experiments  were  continued 
to  later  dates].  It  had  originally  been  intended  to  use  this  date  as  a  "cut-off'  point 
[mainly  because  we  do  not  believe  that  there  has  been  any  significant  advance 
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since  that  time,  see  also  Footnote  E.1J. 


D.4.  The  possibility  of  the  "overfilling"  of  the  cells  was  to  have  been  avoided  by 
adding  an  audio  alarm  system  to  the  "switching  boxes'  which  were  part  of  the 
ICARUS-2  System.  This  modification  was  never  included  and  as  such  it  was 
never  used. 

D.5.  It  was  difficult  to  see  how  anybody  could  convince  themselves  that  such 
measurements  could  give  any  meaningful  results.  The  amount  distilled  was 
determined  by  weighing  the  distillate  and,  in  consequence,  was  reduced  by  the 
hold-up  in  the  vapor  filled  section  of  the  cell  as  well  as  the  condensation  unit. 
Furthermore,  the  rate  of  evaporation  had  to  be  derived  by  differentiating  the 
"noisy"  experimental  data.  However,  let  us  suppose  that  these  defects  could  have 
been  overcome.  In  that  case  we  would  at  best  have  derived  information  about  the 
reflux  ratio,  a  quantity  which  does  not  give  any  useful  information  about  the  rate 
of  excess  enthalpy  generation.  The  only  directly  useable  information  is  the 
detection  of  the  time  at  which  the  cells  are  driven  to  dryness.  However,  this  time 
can  be  determined  directly  from  the  "raw  data"  by  noting  the  fall  in  the  cell 
current  or  by  direct  visual  observation  or,  better,  by  using  video  recordings  [as  we 
had  done  in  1992 (15)(I6)], 

This  part  of  the  instrumentation  was  also  constructed  at  IMRA  Europe  and  the 
effort  involved  caused  considerable  disruption  of  the  research  program  [see 
Footnote  E.2], 

D.6.  Furthermore,  it  is  possible  to  use  filter  functions  other  than  the  "square-box"  filter 
of  the  11-point  averaging  and  to  design  such  filters  to  facilitate  the 
evaluation/demonstration  of  specific  points  of  interest  [e.g.  the  effects  of  positive 
feedback"].  A  program  of  work  on  these  aspects  was  started  in  1992  but  had  to  be 
abandoned  because  of  lack  of  time  and  resources  [see  also  Footnote  E.3]. 

D.7.  See  especially  the  comments  of  Stephen  Jones. 

D.8.  It  is  important  here  to  draw  attention  to  a  shortcoming  of  the  modeling  of  the 
calorimeters,  equation  (B.l),  as  well  as  the  evaluation  of  the  heat  transfer 
coefficients  [especially  of  (kR')n  and  (kR')i 2]*  In  this  modeling  and  in  these 
evaluations  we  assume  that  the  fluctuations  in  the  temperature  about  the  relevant 
mean  in  the  time-domain  are  scaled  by  the  total  water  equivalent  of  the 
calorimeter,  CPM,  to  give  the  fluctuations  in  the  enthalpy  content  of  the  system. 
It  may  well  be  that  this  is  an  overestimate  and  that  the  temperature  fluctuations 
take  place  in  just  a  small  element  of  the  liquid  in  the  calorimeter. 

The  detailed  consideration  of  this  topic  is  beyond  the  scope  of  the  present  report. 
However,  we  should  note  that  the  correct  calculation  of  the  fluctuations  in  the 
enthalpy  content  requires  both  temperature  measurements  at  many  locations 
coupled  to  the  evaluation  of  proper  weighted  means  and  the  development  of  the 


55 


theory  of  the  power  spectral  densities  of  the  fluctuations.  These  topics  were 
investigated  during  1989  but  were  subsequently  abandoned  because  of  lack  of 
resources  and  because  it  became  evident  that  (kR')n  and  (kR')|2  could  be  derived 
within  the  target  error  limits  by  taking  appropriate  averages  [e.g.  (k^~i)11  or 

(kR  ')n  ];  this  was  equally  true  of  the  integral  heat  transfer  coefficients. 
Furthermore,  it  was  also  clear  that  it  would  be  possible  to  obtain  the  required 
spatial  averages  by  using  the  ICARUS- 14  calorimeters  [originally  described  as 
the  ICARUS-4  versions]. 

We  note  here  that  the  work  carried  out  in  1989  relied  on  the  calculation  of  the 
autocorrelation  function  of  the  fluctuations.  However,  it  became  clear  that  it  was 
necessary  to  model  the  cross-correlation  functions  and  that  the  problem  of  the 
generation  of  thermal  noise  in  calorimeters  was  formally  similar  to  the  problem  of 
noise  generation  by  aerofoils. 


Section  E 

For  reasons  which  will  become  apparent,  this  Section  contains  a  set  of  Footnotes  and 
Comments  which  perhaps  should  not  be  included  in  the  main  Report. 


E.  1 .  One  of  the  authors  [M.  Fleischmann]  has  explained  on  other  occasions  that  all  the 
data  sets  which  were  originally  in  his  possession  were  removed  from  the  material 
returned  to  him  from  France.  This  includes  all  the  important  data  sets  collected  in 
the  NHE  Laboratory  prior  to  June  1994,  details  of  a  very  important  experiment 
carried  out  in  another  laboratory  in  Japan  which  mimics  the  procedure  which  had 
been  used  in  Utah  prior  to  October  1989  [i.e.  the  data  analyzed  in  (l2)],  the 
"Harwell"  data  sets  and  all  information  about  the  experiments  carried  out  in  the 
National  Cold  Fusion  Center,  Salt  Lake  City.  The  only  information  now  in  our 
possession  consists  of  parts  of  data  sets,  which  had  been  evaluated  at  various 
times  during  M.  Fleischmann's  stays  in  the  U.K. 

It  might  well  prove  possible  to  reconstruct  at  least  parts  of  the  original  data  sets 
from  this  fragmentary  material.  However,  it  is  not  clear  at  this  time  whether  such 
data  sets  could  be  evaluated  without  access  to  the  relevant  laboratory  notebooks. 
We  note  in  this  connection  that  one  of  us  [M.  Fleischmann]  has  frequently 
requested  access  to  this  information  and  has  suggested  procedures  for  the  further 
evaluations.  These  requests  have  always  either  been  ignored  or  else  blocked. 

E.2.  We  expressed  our  reservations  repeatedly  about  this  part  of  the  research  program 
but  never  received  a  reply.  It  should  be  noted  that  this  development  took  place 
following  completion  [in  June  1994]  of  the  first  part  of  our  investigations  of  the 
long-term  maintenance  of  boiling  conditions  using  the  calorimeters  in  the 
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sequence  ICARUS-4  to  ICARUS-9  (l7).  It  was  this  part  of  the  program,  which 
needed  to  be  developed  further,  and  the  ICARUS  10-13  Calorimeters  were 
designed  and  constructed  with  this  end  in  view.  However,  these  calorimeters 
were  never  put  into  use. 

E.3.  Part  of  this  program  was  aimed  at  self-tuning  control  of  the  systems  with  a  view 
to  implementing  automatic  "hill-climbing"  to  increase  the  enthalpy  output. 

E.4.  This  particular  critique  is  a  minor  [but,  nevertheless  very  good]  example  of  much 
of  the  criticism  in  this  field.  A  scenario  is  set  up  [or  a  point  is  raised]  leading  to 
highly  critical  statements  [should  one  say,  vehemently  critical  statements?]. 
There  is  no  doubt  that  these  statements  would  have  been  justified  if  the  original 
scenario  [or  point]  had  been  valid.  However,  the  original  scenarios  [or  points] 
were  not  valid  and  the  critique  is  therefore  quite  incorrect. 

Other  good  examples  of  such  critiques  include  statements  about  the  inadequacy  of 
stirring,  temperature  inhomogeneities  in  the  cells,  the  definition  and  use  of  pseudo 
conductive  or  pseudoradiative  heat  transfer  coefficients,  the  use  of  a  single 
thermal  impedance  [as  against  two  thermal  impedances  in  series]  the 
methodologies  of  evaluating  the  results  etc.  etc.  All  of  these  initial  scenarios 
should  have  been  checked  against  the  published  papers  or,  perhaps,  even  been 
raised  in  correspondence  prior  to  the  decision  to  publish  the  critiques? 

E.5.  As  has  been  noted  previously  (l3),  it  is  important  to  establish  whether  the  group  at 
NHE  ever  followed  this  particular  instruction  and,  if  they  did  so,  what  conclusions 
they  may  have  drawn  from  such  plots.  However,  it  is  unlikely  that  the  group  at 
NHE  ever  constructed  such  plots  because,  if  they  had  done  so,  they  would  surely 
have  understood  one  reason  for  the  specification  in  the  ICARUS- 1  System  that 
the  duration  of  the  calibration  pulses  had  to  be  lengthened  to  12  hours  [and  the 
reasons  for  other  specifications  of  durations  of  the  experiment  cycles].  They 
would  also  have  detected  the  presence  [or  absence]  of  "positive  feedback".  One 
of  the  authors  [M.H.M.]  reports  that  he  saw  no  evidence  for  the  construction  of 
such  plots  while  he  was  working  at  NHE.  Inadequate  inspection  of  "raw  data"  is 
a  malaise,  which  afflicts  much  of  modern  science. 

E.6.  However,  this  did  not  lead  to  any  discussion  of  the  advantages/short  comings  of 
the  various  methodologies  of  estimating  the  "true  heat  transfer  coefficient".  We 
note  that  the  ICARUS-1  System  Handbook  was  described  as  "Version  1-  Low 
Power  Measuring  System  for  Three  Cells"  with  the  qualification  on  page  62:  "The 
experimental  equipment  and  hardware  are  similar  to  those  which  we  currently 
use.  It  is  envisaged  that  updates  of  the  software  will  be  dependent  on  the  needs 
of  the  laboratories  taking  part  in  this  research  program". 

We  also  note  that  we  never  received  any  communication  from  the  group  at  NHE 
concerning  the  methods  of  data  evaluation.  It  is  evident  that  the  group  continued 
to  base  their  evaluations  on  (k|f)32,  and  it  appears  that  the  value  of  this  coefficient 


57 


determined  from  Day  3  was  also  used  for  the  whole  of  the  evaluation  of  the 
experiment  MC2 1  [Pd-B], 

E.7.  Intentionally? 

E.8.  This  is  an  aspect  of  the  greatest  importance  because  it  shows  that  it  would  enable 
the  precise  and  accurate  evaluation  of  all  the  data  sets,  which  have  been 
determined  hitherto  [see  Footnote  E.9],  In  view  of  the  difficulties,  which  were 
evidently  being  experienced  in  the  data  evaluation,  it  has  been  urged  at  various 
times  that  a  group  should  have  been  set  up  charged  with  the  task  of  evaluating  the 
various  data  sets.  The  task  of  the  parent  groups  would  then  have  been  reduced  to 
the  validation  of  these  evaluations.  There  was  never  any  reply  to  these 
suggestions  [see  Footnote  E.10], 

Short  of  this  suggestion,  one  of  us  [M.  Fleischmann]  urged  that  the  raw  data  for 
various  selected  experiments  should  be  released  for  further  study.  However  with 
the  exception  of  this  experiment  [considered  in  this  Report]  and  the  data  sets 
provided  by  Dr.  Asami  at  ICCF  7,  there  has  again  been  no  response  to  these 
proposals.  We  note  that  the  data  provided  at  ICCF  7  did  not  include  any  of  the 
data  sets  which  had  been  previously  selected.  This  selection  was  based  on  the 
desire  to  restrict  data  evaluations  to  experiments  carried  out  prior  to  June  1994 
[with  the  important  exception  of  this  experiment]. 

E.9.  It  would  evidently  have  been  sensible  to  restrict  such  evaluations  to  just  a  few 
data  sets  because  further  work  should  have  been  carried  out  using  ICARUS- 14 
Calorimeters  [originally  classified  as  the  ICARUS-4  Systems]. 

E.10.  However,  we  note  that  neither  the  original  authors  nor  Harwell  commented  on  the 

reassessments  (20)'(2l)of  parts  of  the  original  investigation  by  the  group  at  Harwell 
(22) 

E.l  1.  Furthermore,  one  of  us  [M.  Fleischmann]  had  believed  for  some  time  that  lines  of 
research  based  on  this  work  could  point  the  way  towards  aspects  of  importance  to 
National  Security  and  should  therefore  be  excluded  from  work  in  the  Public 
Domain  [see  also  Footnote  E.l 2],  It  should  be  emphasized  here  that  the  topic  on 
its  own  does  not  point  to  matters  of  interest  to  National  Security. 

E.l  2.  The  text  of  this  Footnote  may  be  made  available  on  application. 
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cathode  connection 


Fig.  A.l  Schematic  diagram  of  an  ICARUS- 1  Type  Electrochemical  Cell  used  in 
the  Pd-B  experiment.  At  NHE  the  long  thermistor  was  positioned 
somewhat  above  the  anode  spiral.  The  NHE  cell  was  25  cm  in  length  (with 
the  top  8.0  cm  silvered)  and  2.5  cm  in  diameter. 
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.  A. 2  Variation  ot  the  cell  temperature  (A)  and  current  density  (B)  for  the  Pd-B  experiment.  The  dotted  lines  delineate 

regions  for  the  expected  onset  of  significant  levels  of  "positive  feedback"  (A)  and  excess  enthalpy  generation  (B). 


lelineate  the  changes  in  cell  volume  above  which  the  level  of  electrolyte  must  be  assumed  to  approach  the  base  ot  the  Kel-r 
ilug  sealing  the  top  of  the  cell  (Fig.  A.l)  and  below  which  the  level  reaches  the  base  of  the  top  silvered  portion  of  the  Dewar- 


ean  temperature  =  295.198  K 
(temperature  difference)  =  0.0088  K 


Fig.  A. 6  The  cell  temperature  at  the  two  measurement  locations  for  the  first  32,400  s  of  the  first  measurement  cycle  (i.e.  the  period 

0  <  t  <  t}).  Mean  temperature  difference  =  0.0045  K;  a  (temperature  difference)  =  0.0027  K 


Cel)  tempeuture/K 


Fig.  A. 7  The  "raw  data"  (the  cell  temperature  and  cell  potential  as  a  function  of 
time)  for  the  third  measurement  cycle  of  the  Pd-B  experiment. 
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Systems  procedure  with  the  inclusion  of  the  calibration  pulse,  AQ=0.2500  W,  in  the  calculation  of  (k  during  the  period 

R  I  !  °  r 

t  <t<  t  .  It  was  assumed  that  C„M  =  490  JK 

I  2  p 


shown  in  Column  6  of  Table  A. 3  according  to  the  range  of  the  ordinates  and  abscissae  used  to  derive  the  straight  line  fits. 


with  all  remaining  corrections  to  the  calculation  (see  the  main  text).  It  has  been  assumed  that  CDM  =  450  JK  . 


Fig.  A.  14  The  rate  of  excess  enthalpy  generation,  Q,  as  a  function  of  time  for  the 
third  measurement  cycle  as  determined  by  the  ICARUS-Systems 
procedure  using  the  "true  heat  transfer  coefficient,  (k  ’)i2"  ~  0.85065  x 

R 

-9  -4 

10  WK  and  the  values  for  the  "lower  bound  heat  transfer  coefficient, 
(k  shown  in  Fig.  A.  12. 

R 


Fig.  A.  15  The  rate  of  excess  enthalpy  generation,  Q,  as  a  function  of  time  for  the 
third  measurement  cycle  determined  using  the  "true  heat  transfer 

coefficient,  (k  ’)l2"  =  0.85065  x  lO^  WK  .  The  values  of  the  "lower 

bound  heat  transfer  coefficient,  (kR')n,"  have  been  calculated  as  in  the 

ICARUS-Systems  procedure  (see  Fig.  A.  12)  but  with  the  exclusion  of  the 
calibration  pulse,  AQ  =  0.2500W,  during  the  period  t|  <  t  <  t2  (c.f.  the 

NHE  methodology). 
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Fig.  A. 18 


35000 
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0.85065  x  10  WK  as  given  by  the  ICARUS-Systems  analysis  modified  to  take  account  of  "positive  feedback 


0  875 
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Fig.  A. 20  The  variation  with  time  of  the  11-point  average  of  the  "lower  bound  heat  transfer  coefficient,  109(k  for  Day  61  of  this 

R 

experiment.  The  full  line  gives  the  variation  with  time  of  10  (k^)n  determined  for  relevant  "blank"  experiments  4  0^  14). 


absorption  of  deuterium  in  the  lattice  is  40  kJ  Mole  . 


Fig.  A  . 22  Comparison  of  the  specific  rates  of  excess  enthalpy  generation/W  cm  on 
day  68  of  this  experiment  during  the  period  0  <  t  <  21,300  s  (where  the 
cell  is  driven  to  dryness)  and  the  initial  period  21,300  s  <  t  <  30,300  s  of 
observation  of  "Heat  after  Death".  [Due  to  correction  for  electrode 
dimensions,  the  values  for  W/cm3  should  be  multiplied  by  a  factor  0.541. 


■  00  10C 
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The  cooling  curve  observed  on  Day  69  ol  the  operatr 
galvanostat. 


♦  ♦ 
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Fig.  A. 24  The  analysis  of  the  initial  portion  of  the  "cooling  curve"  shown  in  Fig.  A.23  according  to  equation  (A. 3).  The  full  line  shows 
the  R.H.S.  of  this  equation  plotted  with  C,,M  =  28.3  JK  (k  ’).,=  0.65  x  10 9  WR4  and  0  =  295  204  K 

K  u  bath 


Fig.  A.25  The  specific  rate  of  excess  enthalpy  generation/W  cm  for  the  last  part  of 
operation  on  Day  25  and  the  first  part  of  operation  on  Day  26  of  this 
experiment.  Evaluation  according  to  the  ICARUS-Systems  methodology. 
[Due  to  correction  for  electrode  dimensions,  the  values  for  W/cm3  should 
be  multiplied  by  a  factor  0.54]. 


Fig.  A. 26  The  specific  rate  of  excess  enthalpy  generation  for  the  last  part  of 
operation  on  Day  25  and  the  first  part  of  operation  on  Day  26  of  this 
experiment.  Evaluation  given  by  the  group  at  the  NHE  Laboratories.  [Due 
to  correction  for  electrode  dimensions,  the  values  for  W/cm3  should  be 
multiplied  by  a  factor  0.541. 


Fig.  A. 27  Schematic  diagram  of  the  ICARUS-14  Calorimeter.  This  calorimeter  was 
originally  designated  as  the  ICARUS-4  Design. 
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Table  A.l 

Column  1 
Column  2 
Column  3 
Column  4 
Column  5 
Column  6 
Column  7 
Column  8 
Column  9 

Column  10 

Column  1 1 


The  general  operating  characteristics  of  Experiment  M7c2  (Pd-B) 

Time/day 

2 

Current  density/A  cm 

The  maximum  cell  temperature  for  each  day  of  operation/K 

The  volume  of  D2O  added/ml 

The  total  volume  of  D2O  added/ml 

The  volume  of  D2O  electrolysed/ml 

The  total  volume  of  D2O  electrolysed/ml 

The  change  in  the  volume  of  electrolyte  at  the  start  of  the  relevant  day 
The  change  in  the  volume  of  the  electrolyte  at  the  end  of  the  relevant 
day 

The  total  excess  enthalpy  generated  each  day/J  based  on  the  "true  heat 
transfer  coefficient,  (kR')]2  =  0.79350  x  10'9  WK 4” 

The  total  excess  enthalpy  generated  each  day/J  based  on  the  "true  heat 
transfer  coefficient,  (kR’)12  =  0.85065  x  10~9  WK'4” 
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Table  A.  1 
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Table  A. 2  The  "raw  data  for  the  third  measurement  cycle  of  experiment  M7c2 

and  the  derived  "lower  bound  heat  transfer  coefficient,  (kR’)  j  j "  and  rates  of  excess 
enthalpy  production,  Q,  as  given  by  the  analysis  provided  by  the  group  at  the  N.H.E. 
Laboratory.  Values  ot  Q  have  been  derived  using  the  '‘true  heat  transfer  coefficient, 
(kR’)l2"  (°-793504  x  10’9  WKT4)  quoted  by  this  group. 

The  Table  also  lists  the  averages  Q  and  Q  as  well  as  109(kj^’)j  j  and 
1 09(kj^')  j  j  (and  the  associated  statistics)  for  the  evaluations  based  on  measurements 
with  the  short  thermistor. 

Columns  15-17  give  the  evaluation  of  109(kj^')j  j,  109(kj^') j  j  and 
109(kj^')jj  but  with  the  inclusion  of  the  enthalpy  input  AQ  =  0.2500W  due  to  the 
calibration  heater. 


Column  1 
Column  2 
Column  3 
Column  4 
Column  5 
Column  6 
Column  7 
Column  8 
Column  9 
Column  10 
Column  1 1 
Column  12 
Column  13 
Column  14 


The  elapsed  time/s 

The  temperature  of  the  water  thermostat/K 
The  temperature  measured  with  the  "short  thermistor"/K 
The  temperature  measured  with  the  "long  thermistor"/K 
The  Cell  Voltage/V 
The  Cell  Current/A 

109(Lr')j  ,/WK-4  for  measurements  in  Column  3 

Q/W  for  measurements  in  Column  3 

109(kj^')j  j/WIC4  for  measurements  in  Column  4 

Q/W  for  measurements  in  Column  4 

Q/W  of  the  values  in  Column  8 

Q/W  of  the  values  in  Column  1 1 

109(kj^')j  j/WK"4  o/WK"4  and  a/(kj^')j  j/%  of  the  values  in  Column  9 

109(kj^')j  j/WK"4  a/ WK‘4  and  a/(kj^')]j/%  of  the  values  in  Column  13 
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Column  15 


Column  16 


Column  17 


1  °9(kR‘)  1  ]/WK‘4  for  measurements  in  Column  3  but  with  inclusion  of 
the  enthalpy  input  AQ  =  0.2500W  due  to  the  calibration  heater. 

1  °9(^r')]  ]AVK'4  o/WK'4  and  a/(kR')]j/%  of  the  values  in  Column 

15 

]  °9(kR')i  i  AVK-4  o/WK"4  and  a/(k^)n/%  of  the  values  in  Column 

16 
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Table  A.3  Abridged  (kR’)]j  spreadsheet  for  the  third  measurement  cycle  of 
experiment  M7c2  .  The  evaluation  shown  in  Columns  2  -  8  has  included  all  the 
corrections  to  that  given  by  the  group  in  the  N.H.E.  Laboratories  with  the  exception 
of  the  use  of  CpM  =  490  JK"^.  The  evaluation  in  Columns  9-12  includes  the 
correction  of  this  parameter  to  CpM  =  450  JK'* 


Column  1 
Column  2 
Column  3 
Column  4 
Column  5 


Column  6 


Column  7 
Column  8 
Column  9 
Column  10 
Column  1 1 
Column  12 


The  elapsed  time/s 
The  total  input  power/W 
Values  of  -CpM(d0/dt)/f(0)  AVK"4 
Values  of  input  power/(0)  AVK"4 

The  derived  values  of  109(kR')j6j,  1  09(kR’) j 7 j AVK'4,  CpM/JK'^ 
and  the  relevant  regression  coefficients.  The  arrows  indicate  the 
ranges  of  the  fitting  procedures 

The  derived  values  of  1  09(kR’)  j  g  j  AVK-4,  CpM/JK"4  and  the  relevant 
regression  coefficients.  The  arrows  indicate  the  ranges  of  the  fitting 
procedures;  the  dotted  line  indicates  the  linkage  of  the  data  for 
measurements  in  t  j  <  t  <  t2  with  those  for  t2  <  T. 

109(kR')l  1  derived  from  Columns  3  and  4 

The  1 1 -point  running  mean  of  109(kR’)j  j 

Values  of-450(d9/dt)/f(e)  AVK-4 

Values  of  109(kR’)j  j  AVK-4  based  on  Columns  4  and  9. 

The  1 1 -point  running  mean  AVK"4  of  109(kR')j  j  shown  in  Column  10 
The  rate  of  excess  enthalpy  generationAV  based  on  the  "true  heat 
transfer  coefficient,  (kR’)]2"  =  0.85065  x  10"9  WK"4 
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18600  0.45467  0.00752  0.83675  _  0.84127  0.83932  0.00691  0.84366  0.83934  -0.00024 

18900  0.45515  0.00452  0.8391  _ 0.84362  0.84059  0,00415  .  0,84325  0.84056  -0.00002 

192001  0.455181  -0.00151 1  0.838831  I  1  0.83732  0.84041  -0.00139  0.83744  0.84041  0.00313 
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Table  A.4  The  combined  (kR')2]  and  (kR')3]  spreadsheet  for  the  whole  of  the 
third  measurement  cycle  of  experiment  M7c2 


Column  ] 
Column  2 

Column  3 
Column  4 
Column  5 

Column  6 

Column  7 

Column  8 
Column  9 
Column  10 

Column  1 1 


Column  12 
Column  13 


The  elapsed  time/s 

109  CpM  (0  -  0O)/  / f(0)dx  AVK"4 

Here  0Q  =  300.3 145K  (the  average  of  the  last  1 1  measurements) 

10 9  f input  dx/  f f(0)dx  AVK'4 
lo9(kR’)21  AVK"4 

]09(kR')25i  CpM/JK'^,  correlation  coefficient. 

The  arrows  indicate  the  ranges  of  the  fitting  procedures. 

109(kR')2(3i  CpM/JK'*,  correlation  coefficient. 

The  arrows  indicate  the  ranges  of  the  fitting  procedures. 

109  CpM(6-60)fe)dz  avk-4 

Here  0Q  =  300.353K  (the  initial  value  of  the  cell  temperature) 
lO^inputdx/  / f(0)dx  AVK"4 
]°9  (kR’)3  ]  AVK"4 

109(kR')35]  CpM/JK"*,  correlation  coefficient. 

The  arrows  indicate  the  ranges  of  the  fitting  procedures. 

109  CpM(0-0Q)A  ft f(0)dx  AVK-4  (evaluation  of  (kR^). 

Here  0O  —  300.3 175K,  the  average  of  the  11  measurements  preceding 
the  application  of  the  calibration  pulse. 

109  A  f  input  dz/A  f({Q)dx  AVK"4  (evaluation  of  (kR’)32) 
lo9(kR’)32  AVK"4 
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50400  0,06785  0.75805  0,8259 _ 0.04231  0.S6802  0.825721  0.306571  1.06197  0,7554 

50700  0.06901  0.75733  0,82634 _ 0.04214  0,86759  0.82544 _ 0,30102  1.05545  0.75444 

510001  0.069891  0.756591  0.82648 _ 0.04179  0.86716  0.82537  0.29442  1.04925  0.75482 
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Table  A.5  The  combined  abridged  (kR')2]  and  (kR’)31 -spreadsheets  prepared 
according  to  the  instructions  in  the  ICARUS-Systems  Handbooks^1^2)  (restriction  of 
the  range  of  the  integrations  to  the  region  of  application  of  the  calibration  pulse  t  j  <  t 
<  t2).  The  third  measurement  cycle  of  experiment  M7c2  .  The  evaluation  of 
(kR')3  ] ,  (kR')351>  (kR')32*  (kR')362>  (kR'>22  and  (kR’>262-  Modification  of  the 
procedure  for  the  evaluation  of  (kR')22  and  (kR‘)262  to  take  account  of  the  effects  of 
"positive  feedback"  and  evaluation  of  these  coefficients. 

Column  1  The  elapsed  time/s  (from  the  start  of  the  measurement  cycle) 

Column  2  1 09  CpM  (0-0o)/  f  *0)dx/  WK-4 

Here  0Q  =  300.3 175K,  the  average  of  the  1 1  measurements  preceding 
the  application  of  the  calibration  pulse 
Column  3  1 09^*input  dx/^f(0)dx  /WK-4 

Column  4  109(kR')31  AVK-4 

Column  5  1 09(kR')36  j  /WK-4  CpM/JK' 1 ,  correlation  coefficient. 

The  arrows  indicate  the  range  of  the  fitting  procedure. 

Column  6  1 09  CpM  (0-0Q)/A  J f(0)dx  /WK-4  (evaluation  of  (kR’)32  and 

(kR')362) 

Column  7  1 09  A ^input  dx/A >Ji f0)dx  /WK'4(evaluation  of  (kR’)32  and  (kR')3(52) 

Column  8  1 09  (kR')32  /WK-4 

Column  9  1 09  (kR')362  /WK-4,  CpM/JK-1,  correlation  coefficient. 

The  arrows  indicate  the  ranges  of  the  fitting  procedures 
Column  10  109  CpM  (0-0o)/A  J  f(0)dx  /WK-4  (evaluation  of  (kR')22  and 

(kR’)262) 

Here  0O  =  303.0747K,  the  average  of  the  last  1 1  measurements  during 
the  application  of  the  calibration  pulse. 

Column  11  1 09  A J\ input  dx/A  J  f(0)dx  AVK-4 


Column  12 


Column  13 


Column  14 


Column  15 


Column  16 


109(kR')22/WK-4 

109(kR')2£2  /WK-4  CpM/TK1,  correlation  coefficient. 

The  arrows  indicate  the  range  of  the  fitting  procedures. 

Modification  of  column  1 1  to  take  into  account  the  effects  of  "positive 
feedback" 

Values  of  lO^k^'^  taking  into  account  the  effects  of  "positive 
feedback" 

lo9(kR')262  /WK"4>  CpM/JK'1,  correlation  coefficient  taking  into 
account  the  effects  of  "positive  feedback". 
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Table  A.  5 
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Table  A.6 

Column  1 
Column  2 
Column  3 
Column  4 
Column  5 
Column  6 
Column  7 
Column  8 
Column  9 

Column  10 


Abridged  (kR')j  j -spreadsheet  for  the  operation  of  experiment  M7c2 
on  Day  61. 

The  elapsed  time/s 

The  cell  temperature/K 

The  cell  voltage/V 

The  cell  current/A 

The  rate  of  evaporative  coolingAV 

f(0)  /K4 

The  "lower  bound  heat  transfer  coefficient,  1  09(kR') j  j "  AVK"4 
109(k^T')i  j  AVK"4,  oAVK"4  and  a/t^n  /% 

The  "lower  bound  heat  transfer  coefficient,  109(kR')jj"  AVK-4  as 
given  by  the  evaluation  due  to  the  group  at  the  N.H.E.  Laboratories 
109(kJ)n  AVK'4,  aAVK"4  and  a/(kj)n  /%  of  the  values  of 
109(kR’)j  j  as  given  by  the  evaluation  due  to  the  group  at  the  N.H.E. 
Laboratories 
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Table  A.6 
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Table  A.7  The  operation  of  experiment  M7c2  on  Day  68  during  the  period  0  <  t 
<  21,300s  when  the  cell  is  driven  to  dryness.  The  evaluation  given  by  the  group  at  the 
N.H.E.  Laboratories. 


Column  1 
Column  2 
Column  3 
Column  4 
Column  5 
Column  6 


Column  7 


Column  8 


Column  9 


The  elapsed  time/s 

The  temperature  measured  by  the  "short  thermistor’VK 
The  temperature  measured  by  the  "long  thermistor’VK 
The  Cell  Voltage/V 
The  Cell  Current/A 

The  "lower  bound  heat  transfer  coefficient,  lO^kj^jj"  AVK-^  given 
by  the  evaluation  by  the  group  at  the  N.H.E.  Laboratories  and  based  on 
the  temperature  measurements  with  the  "short  thermistor" 

The  rate  of  excess  enthalpy  generation,  Q,  AV  given  by  the  evaluation 
by  the  group  at  the  N.H.E.  Laboratories  and  based  on  the  temperature 
measurements  with  the  "short  thermistor" 

The  "lower  bound  heat  transfer  coefficient,  1 j  j "  AVK'^  given 
by  the  evaluation  by  the  group  at  the  N.H.E.  Laboratories  and  based  on 
the  temperature  measurements  with  the  "long  thermistor" 

The  rate  of  excess  enthalpy  generation,  Q,  AV  given  by  the  evaluation 
by  the  group  at  the  N.H.E.  Laboratories  and  based  on  the  temperature 
measurements  with  the  "long  thermistor". 
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133 


18600  374.08  374.148  23.8283  0.9857  0.43198  4.3332  0.209  7.0142 

18900  374,158  374.184  25.8979  0.98207  0.33529 _ 5.5  0.22555  6.8202 

19200  374.181  374.213  28.8106  0.9772  0.47236  3.8563  0.31437  5.7566 


Table  A.8  The  operation  of  experiment  M7c2  on  Day  68  during  the  period  0  <  t 
<  21,300s  when  the  cell  is  driven  to  dryness.  Evaluation  by  the  ICARUS-Systems 
methodology  but  neglecting  the  effects  of  reflux  D2O  in  the  cell  and  with  P*  —  1 
atmosphere. 


Column  1 
Column  2 
Column  3 
Column  4 
Column  5 

Column  6 


Column  7 


Column  8 
Column  9 

Column  10 


Column  1 1 


The  elapsed  time/s 

The  enthalpy  input  AV 

The  enthalpy  change  of  the  calorimeter  AV 

f(0)  /K4  as  given  by  measurements  with  the  "short  thermistor" 

The  rate  of  evaporative  cooling  based  on  temperature  measurements 
with  the  "short  thermistor" 

The  "lower  bound  heat  transfer  coefficient,  10^(kj^')|  j "  AVK”4  given 
by  the  ICARUS-Systems  evaluation  and  based  on  measurements  with 
the  "short  thermistor" 

The  rate  of  excess  enthalpy  generation,  Q,  AV  given  by  the  ICARUS- 
Systems  evaluation  and  based  on  measurements  with  the  "short 
thermistor" 

f(0)  /K4  as  given  by  measurements  with  the  "long  thermistor" 

The  rate  of  evaporative  cooling  based  on  temperature  measurements 
with  the  "long  thermistor" 

The  "lower  bound  heat  transfer  coefficient,  1 0^(kj^')  j  j "  AVK-4  given 
by  the  ICARUS-Systems  evaluation  and  based  on  temperature 
measurements  with  the  "long  thermistor" 

The  rate  of  excess  enthalpy  generation,  Q,  AV  given  by  the  ICARUS- 
Systems  evaluation  and  based  on  measurements  with  the  "long 
thermistor" 
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18900  23.93393  0.07575  12.00282  21.82352  0.33314  6.21156  12.00827  22.88811  0.08078  9.24481 

19200  26,66153  0.07876  12.00785  22.59757  0.33188  6.22928  12.01455  24.63937  0,16175  8.27677 

19500  30.44557  0.13875  12.02483  28.46785  0.15293  8.38995  12.03238  32.27347  -0.16345  12.202 
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Table  A .9  The  operation  of  experiment  M7c2  on  day  68  during  the  period  0  <  t  < 
21,300s  when  the  cell  is  driven  to  dryness.  Evaluation  by  the  ICARUS-Systems 
methodology  with  the  inclusion  of  the  effects  of  reflux  of  D2O  in  the  cell.  P* 
assumed  to  be  equal  to  the  pressure  measured  at  the  same  time  at  Sapporo  Airport. 


Column  1 
Column  2 

Column  3 

Column  4 

Column  5 

Column  6 

Column  7 

Column  8 

Column  9 

Column  10 

Column  1 1 


The  elapsed  time/s 

The  rate  of  evaporative  cooling  AV  based  on  temperature 
measurements  with  the  "short  thermistor" 

The  amount  of  D2O  evaporated/Mole  in  each  300s  interval 
corresponding  to  the  rate  of  evaporative  cooling  shown  in  Column  2 
The  amount  of  D20  returned  to  the  cell  by  reflux/Mole  corresponding 
to  the  rate  of  evaporative  cooling  shown  in  Column  2 
The  amount  of  D20  removed  from  the  cell/Mole  corresponding  to  the 
rate  of  evaporative  cooling  shown  in  Column  2 

The  rate  of  evaporative  cooling  AV  based  on  temperature 
measurements  with  the  "long  thermistor" 

The  amount  of  D2O  evaporated/Mole  in  each  300s  interval 
corresponding  to  the  rate  of  evaporative  cooling  shown  in  Column  6 
The  amount  of  D20  returned  to  the  cell  by  reflux/Mole  corresponding 
to  the  rate  of  evaporative  cooling  shown  in  Column  6 
The  amount  of  D20  removed  from  the  cell/Mole  corresponding  to  the 
rate  of  cooling  shown  in  Column  6 

The  "lower  bound  heat  transfer  coefficient,  10^(kj^')j  j"  /WK"^  given 
by  the  ICARUS-Systems  evaluation  and  based  on  measurements  with 
the  "long  thermistor" 

The  rate  of  excess  enthalpy  generation,  Q,  /W  given  by  the  ICARUS- 
Systems  evaluation  and  based  on  measurements  with  the  "long 


thermistor" 
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141 


18900  27.9715  0.201362  0.040209  2.878859  27.9774  0.201404  0.043696  3.109763  0.16244  8.2642 

19200  28.3787  0.204293  0.042611  3.040541  28.5710  0.205677  0.046054  3.269386  0.36699  5.8109 

19500  35.2563  0.253803  0.045068  3.249276  35.4106  0.254914  0.04849  3.47581  0.13564  8.6033 


Table  A. 10  The  operation  of  experiment  M7c2  on  day  68  during  the  period 

21,300s  <  t  <  30,300s  following  evaporation  to  dryness. 


Column  1 

The  elapsed  time/s 

Column  2 

The  cell  temperature/K 

Column  3 

The  cell  voltage/W 

Column  4 

The  cell  current/A 

Column  5 

The  enthalpy  input  AV 

Column  6 

The  rate  of  change  of  the  enthalpy  of  the  calorimeter  AV 

Column  7 

The  rate  specific  of  excess  enthalpy  generation  AV  based  on 

the  "true  heat  transfer  coefficient,  =  0.85065  x  10'^  WK"^ 

Column  8 

The  specific  rate  of  excess  enthalpy  generation  AV  based  on  an 
assumed  "true  heat  transfer  coefficient,  (kj^,)]2"  =  0.65  x  10  ^  WK  ^ 
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Table  A.  10 


( 

I 


I 
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Table  A.  11 
day  69. 

Column  1 
Column  2 
Column  3 

Column  4 
Column  5 
Column  6 
Column  7 

Column  8 


The  operation  of  cell  M7c2  during  the  period  2400s  <  t  <  32,400s  of 


The  elapsed  time/s 

The  cell  temperature  measured  by  the  "long  thermistor"  /K 
The  first  term  on  the  L.H.S.  of  equation  (A.3). 

The  mean  temperature  of  the  water  bath  was  65^  =  295,204K 
The  second  term  on  the  L.H.S.  of  equation  (A.3) 

The  sum  of  Columns  3  and  4 

The  rate  of  change  of  the  enthalpy  of  the  calorimeter  AV 
The  rate  of  enthalpy  output  from  the  calorimeter  AV. 

The  "true  heat  transfer  coefficient,  has  been  assumed  to  be 

that  for  a  cell  filled  with  gas  alone  i.e.  0.85065  x  10“^  WK-^ 

The  rate  of  excess  enthalpy  generation  AV 
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Table  All 
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18900  295.25  5.61897  -0.04278  5.57619 

19200  295.263  5.37009  -0.04274  ~~1 32736 

195001  295.3951  4.19558  -0.04224  4.15329 


Table  A.12 

Column  1 
Column  2 
Column  3 
Column  4 
Column  5 
Column  6 
Column  7 
Column  8 
Column  9 
Column  10 
Column  1 1 

Column  12 


The  final  part  76,800s  <  I  <  86,400s  and  the  first  part  0  <  t  <  21,300s 
of  the  operation  of  experiment  M7c2  on  Days  25  and  26  respectively 

The  elapsed  time/day  +  5 

The  temperature  measured  by  the  "short  thermistor"  /K 

The  cell  voltage  /V 

The  cell  current/A 

The  rate  of  enthalpy  input  AV 

The  rate  of  evaporative  cooling  AV 

The  rate  of  change  of  the  enthalpy  content  of  the  calorimeter  AV 
f(0)  /K4 

The  "lower  bound  heat  transfer  coefficient,  10^(kj^')j  j"  AVK"4 
The  specific  rate  of  excess  enthalpy  generation  AV  cm“^ 

The  "lower  bound  heat  transfer  coefficient,  1 0^(kj^')  j  j "  AVK^4  given 
by  the  evaluation  of  the  group  at  the  N.H.E.  Laboratories 
The  specific  rate  of  excess  enthalpy  generation  AV  cm"'*  given  by  the 
evaluation  of  the  group  at  the  N.H.E.  Laboratories. 
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